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FOREWORD 

In  connection  with  the  widespread  use  of  supersonic  flying 
craft,  close  attention  has  been  given  to  problems  associated  with 
aerodynamic  and  radiant  heating  in  flight.  Of  particular  importance 
is  the  heat  shielding  of  the  crew,  of  the  structure  and  equipment, 
the  selection  of  structural  materials,  and  the  determination  of  the 
thermal  stresses  in  the  structure.  To  solve  these  problems,  it  is 
necessary  first  to  determine  the  temperature  of  the  exterior  shield¬ 
ing.  This  is  the  purpose  of  this  work. 

In  considering  aerodynamic  heat  transfer,  we  took  into  con¬ 
sideration  the  entire  range  of  velocities  and  altitudes  possible 
during  flight  in  the  atmosphere.  However,  not  all  the  problems 
dealt  with  in  this  connection  were  given  equal  attention.  This  is 
explained  by  the  limited  attention  devoted  to  certain  questions  in 
the  scientific  literature  and  by  the  fact  that  the  work  is  Intended, 
for  the  most  part,  for  educational  purposes. 

In  examining  the  problems  of  radiation  heat  transfer,  primary 
attention  is  given  to  emission,  since  it  is  of  great  significance 
from  the  standpoint  of  thermal  balance  in  aerodynamic  as  well  as  in 
radiant  heating. 

In  connection  with  the  Increasingly  wide  use  of  shielding 
materials  with  low  thermal  conductivity,  and  also  in  connection 
with  the  use  of  thick  shielding  in  the  case  of  intensive  heating, 
the  present  work  presents  a  method  for  determining  the  temperatures 
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of  thin  as  well  as  thick  shielding.  In  addition  the  criterion  for 
thin  shielding  is  presented,  as  is  a  method  for  determining  the  thick¬ 
ness  of  the  required  insulation. 

The  present  work  is  intended  for  preliminary  tentative  cal¬ 
culations;  consequently,  in  the  majority  of  cases  we  present  approxi¬ 
mation  methods  of  calculation  in  order  to  obtain  a  quick  result.  In 
connection  with  the  fact  that  we  are  pirmarlly  concerned  with  edu¬ 
cational  purposes,  we  consider  only  the  simplest  design  elements 
(plate,  shape,  delta  wing,  cone,  region  of  the  critical  point  of 
the  forward  blunt  end  of  the  body) . 

Although  this  work  takes  into  consideration  a  wide  range  of 
velocities,  it  should  be  borne  in  mind  that  the  accuracy  of  the  cal¬ 
culations  with  respect  to  heat  transfer  will  decrease  as  the  Mach 
numbers  Increase,  since  such  factors  as  dissociation,  recombina¬ 
tion,  and  the  development  of  the  thickness  of  the  boundary  layer 
can  be  calculated  only  with  a  relatively  low  degree  of  accuracy. 

The  radiation  of  the  air  behind  the  shock  wave  at  high  Mach  numbers 
may  lead  to  greater  heat  transfer  than  in  the  case  of  aerodynamic 
heat  transfer  and  establishes  the  applicability  limit  of  the  above 
methods  for  the  determination  of  the  temperature  at  the  critical 
point.  The  radiation  of  the  air  behind  the  shock  wave  is  given  only 
qualitative  consideration. 

In  this  work  it  has  not  been  our  Intention  to  examine  methods 
of  heat  shielding  and  consequently  certain  problems  of  heat  trans¬ 
fer  closely  related  to  measures  of  heat  protection  have  not  been 
dealt  with.  Heat  transfer  in  the  case  of  porous  cooling  and  sub¬ 
limation  Is  related  to  these  questions.  We  have  also  not  considered 
heat  transfer  through  multilayered  shielding. 

To  better  explain  the  complex  methods  of  thermal  calculations, 
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this  book  presents  numerical  examples  for  which  the  initial  data 
./ere  arbitrarily  adopted. 
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SYIVBOLS 


A  —  constant;  accomodation  coefficient; 

p 

a  —  coefficient  of  thermal  dlffuslvity,  In  m  /sec;  speed  of 
~  sound.  In  m/sec;  constant; 

c_  —  specific  heat  capacity,  in  kcal/kg-deg; 

c  —  specific  heat  capacity  of  gas  at  constant  pressure,  in 
P  kcal/kg-deg; 

c^  —  local  coefficient  of  aerodynamic  friction; 

E  —  energy  per  unit  time  in  kg-m/sec;  total  radiant  energy 
per  unit  surface,  per  unit  time,  in  kg/sec -m; 

2 

F  —  area,  in  m  ; 

G  —  weight,  in  kg; 

g  —  acceleration  of  gravity,  in  m/sec  ; 

H  —  flight  altitude  in  m;  ratio  of  displacement  thicknesses 
and  momentum  losses  for  an  incompressible  fluid; 

—  Mach  number  of  an  undisturbed  flow; 

Mg  —  local  Mach  number; 

Pr  —  Prandtl  number; 

2 

p^  —  pressure  in  undisturbed  flow,  in  kg/m  ; 

p.  —  relative  pressure  at  the  surface,  referred  to  the  impact 
pressure  of  undisturbed  flow; 

Pg  —  relative  pressure  at  the  surface,  referred  to  local  Impact 
pressure; 

p 

q  —  specific  heat  flow,  in  kcal/m  sec; 

q  -  specific  heat  flow,  developing  as  a  result  of  aerodynamic 
^  heating,  in  kcal/m^sec; 

q^  —  specific  heat  flow,  developing  as  a  result  of  aerodynamic 
^  heating  at  the  critical  point,  in  kcal/m^sec; 
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—  specific  heat  flow  of  radiant  energy  absorbed  by  a  body. 

—  In  kcal/m2sec; 

O 

—  specific  heat  flow  radiated  by  a  body,  in  kcal/m  sec; 

q  —  specific  heat  flow  involved  In  the  heating  of  a  body,  in 
®  kcal/m^sec; 

q,  —  flow  of  heat  removed  as  a  result  of  the  heat  conduction 

beyond  the  boundary  of  the  body  under  consideration,  referred 
to  1  m2  of  the  outer  surface  of  a  body,  in  kcal/m2sec; 

O 

'^pog  ~  specific  heat  flow  absorbed  by  a  body,  in  kcal/m  sec; 

O 

dokh  ~  specific  heat  flow  of  cooling,  in  kcal/m  sec; 

R  —  radius  of  surface  curvature,  in  m; 

Re  —  Reynolds  number; 

Re^  —  local  Reynolds  number; 

Rej^j,  —  critical  Reynolds  number; 

Re*  —  Reynolds  number  corresponding  to  the  determining  temperature; 
r  —  temperature  recovery  factor; 

p 

S  —  specific  heat  flow  of  solar  radiation,  in  kcal/m  sec; 

St  —  Stanton  number; 

£  —  proportionality  factor  for  the  Stanton  number; 

T  —  absolute  temperature,  in  °K; 

Tjjg  —  temperature  of  an  undisturbed  flow,  in  °K; 

Tp  —  temperature  of  the  adiabatic  wall,  in  °K; 

Ti  —  temperature  of  the  air  at  the  outer  limit  of  the  boundary 
layer  (local  temperature),  in  °K; 

Tj^  —  temperature  of  the  inner  structure,  in  °K; 

—  temperature  of  the  shielding,  in  °K; 

—  temperature  of  the  outer  surface  of  the  shielding,  in  °K; 

—  temperature  of  the  inner  surface  of  the  shielding,  in  °K; 

T*  —  determining  temperature,  in  °K; 

Ug^  —  rate  of  recombination  at  the  wall,  in  cm/sec; 

V  —  velocity,  in  m/sec; 
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—  local  velocity  at  the  limit  of  the  boundary  layer.  In  m/sec 

—  undisturbed  flow  velocity,  In  m/secj 

Vko3  —  first  cosmic  velocity  (circular  velocity).  In  m/secj 

X  —  distance  from  the  leading  edge  or  the  nose.  In  m; 

X  —  relative  coordinate  from  the  nose  or  the  leading  edge; 

y  —  coordinate  normal  to  the  surface.  In  m; 

z  —  coordinate  along  the  surface,  perpendicular  to  the  flow. 

In  m; 

p 

a  —  heat- transfer  coefficient,  in  kcal/m  sec -deg; 

—  heat- transfer  coefficient  at  the  critical  point  of  the 
blunt  nose.  In  kcal/m2sec -deg; 

a.  —  heat- transfer  coefficient  referred  to  Increment  in  heat 
content.  In  kg/m^sec; 

p 

a^k  “  equivalent  coefficient  of  heat  transfer.  In  kcal/m  sec -deg; 

P  —  coefficient  of  the  Integral  absorption  of  radiant  energy; 
the  Blot  criterion  of  similarity; 

Y  —  bulk  weight,  in  kg/m^; 

6  —  thickness  of  the  boundary  layer,  in  m; 

^ob  “  thickness  of  the  shielding,  in  m; 

6*  —  thickness  of  the  displacement  of  the  boundary  layer,  in  m; 

6  —  thickness  of  momentum  loss  of  the  boundary  layer.  In  m; 

e  —  coefficient  of  integral  radiation  from  a  smooth  surface; 

—  coefficient  of  monochromatic  radiation  of  wavelength  X; 

e  —  relative  teiqjerature  difference  along  the  thickness  of  the 
P  shielding; 

0  —  relative  temperature; 

i^k  —  half-angle  of  conic  flare,  in  radians; 

K.  —  adiabatic  exponent; 

X  —  coefficient  of  thermal  conductivity,  in  kcal/m*sec -deg; 
wavelength; 

Xjj^  —  wavelength  of  maximum  intensity,  in  M-j 

P 

—  coefficient  of  viscosity,  in  kg«sec/m  ; 


-  6  - 


4 

M-  —  coefficient  of  viscosity  in  the  case  of  the  determining 
temperature,  in  kg-sec/m^; 

p 

V  —  coefficient  of  kinematic  viscosity,  in  m  /sec; 

V  —  coefficient  of  kinematic  viscosity  in  the  case  of  the 

determining  temperature,  in  m^/sec; 

i  —  similarity  parameter  for  conic  flows;  relative  coordinate 

Pg  —  air  density  of  the  outer  limit  of  the  boundary  layer,  in 

kg •sec2/m; 

2  /  U 

—  air  density  in  undisturbed  flow,  in  kg'sec  /m  ; 

p*  —  air  density  in  the  case  of  the  determining  temperature, 
in  kg'secVm^; 

/  2  4 

0  —  Stefan- Boltzmann  constant,  in  kcal/m  sec -deg  ; 

T  —  time,  in  seconds; 

q)  —  Fourier  criterion;  the  angle  between  the  normal  to  the 
surface  of  the  body  and  the  vertical,  in  radians;  the 
angle  between  the  direction  of  flow  and  the  tangent  to 
the  surface  of  the  wing,  in  radians; 

X  —  coefficient  characterizing  the  distribution  of  velocities 
near  the  critical  point;  wing  sweepback  angle,  in  radians 

f  —  angle  between  the  direction  of  the  solar  rays  and  the 
normal  to  the  surface,  in  degrees; 

CD  —  exponent . 
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Chapter  1 

HEAT  TRANSFER  IN  AERODYNAMIC  HEATING 
§  l.I.  HEAT  PLOWS  DEVELOPING  IN  PLIGHT  AT  THE  SURFACE  OP  A  BODY 

In  the  case  of  constant  air  flow  about  any  type  of  body,  the 
flow  at  its  walls  will  have  a  velocity  equal  to  zero.  The  kinetic 
energy  of  the  flow  in  this  instance  is  converted  into  heat  and  the 
air  at  the  wall  is  heated. 

In  the  case  of  complete  conversion  of  the  kinetic  energy  of 
the  flow  into  thermal  energy,  the  decelerated  flow  temperature 
would  be 

However,  with  the  conversion  of  kinetic  energy  into  thermal 
energy,  a  part  of  this  energy  is  dissipated  and  the  temperature  of 
the  air  at  the  adiabatic  wall,  i.e.,  that  part  of  the  wall  neither 
absorbing  nor  radiating  heat,  will  be: 

where  r  is  the  temperature  recovery  factor 

The  temperature  T  is  called  the  recovery  temperature  or  the 

r 

temperature  of  the  adiabatic  wall,  since  the  temperature  of  the 
latter  will  actually  equal  T^. 

The  value  of  the  recovery  factor  depends  in  the  main  on  the 
structure  of  the  boundary  layer  and  on  the  physical  properties  of 
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(1.1) 

(1.2) 


the  air  which  are  determined  by  the  Prandtl  number 


where  a  is  the  coefficient  of  thermal  diffusivity. 

The  engineering  calculations  may  be  assumed  in  approximate  terms 
for  the  laminar  boundary  layer 

^  =  (Pr>=)’'’.  (1.4) 

for  the  turbulent  boundary  layer 

/•.-(IV”')''',  (1.5) 

where  the  Prandtl  number  should  correspond  to  the  determining  tem¬ 
perature  which  may  be  determined  from  the  formula  given  in  the  work 
of  Eckert  [35]: 

7-^ -  n  !-  0,5  ('/■„ -  7\)  H-  0,22  (7>  -  7,).  (1,6) 

The  Introduction  of  the  determining  temperature  produces  great 
differences  in  the  temperature  of  the  air  along  the  thickness  of 
the  boundary  layer.  It  leads  to  a  situation  in  which  the  physical 
properties  of  the  air  (viscosity,  density,  and  thermal  conductivity) 
along  the  thickness  of  the  boundary  layer  are  variable.  In  principle, 
it  is  possible  to  solve  problems  of  heat  transfer  without  the 
determining  temperature  by  means  of  the  corresponding  integration 
of  the  equations  of  the  boundary  layer.  However,  for  a  compressible 
gas  even  in  the  case  of  a  laminar  boundary  layer,  it  is  necessary 
to  carry  out  the  unwieldy  numerical  integration  of  two  second-order 
equations.  A.D.  Young  [30]  showed  that  the  integration  of  the 
equations  for  a  laminar  boundary  layer  of  a  compressible  gas  may 
be  substantially  simplified,  if  the  viscosity  is  represented  in 
the  following  form; 

In  this  case  we  will  obtain  results  even  in  the  form  of  formulas, 
although  admittedly  these  are  quite  cumbersome.  However,  use  of  the 
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above  formula  for  viscosity  Is  possible  only  within  a  comparatively 
narrow  temperature  range  in  the  boundary  layer,  since  cu  is  not  a 
constant.  Actually,  in  the  absolute  temperature  range  from  300  to 
500°K,  OD  =  0.72,  while  in  the  temperature  range  from  1000  to  2000°K, 

00  =  0.58. 

In  the  presence  of  a  turbulent  boundary  layer  it  is  not  possible 
to  obtain  simplified  solutions,  even  of  this  type.  In  connection 
with  the  difficulty  of  solving  boundary  layer  equations,  the  method 
of  the  determining  temperature  is  fruitful.  With  this  method  we  may 
use  the  results  of  the  solution  of  boundary  layer  equations  for  an 
incompressible  gas,  assuming  the  physical  properties  of  the  air  for 
the  determining  temperature. 

The  value  of  the  determining  temperature  depends  not  only  on 
the  temperatures  in  Formula  (1.6),  but  on  the  structure  of  the 
boundary  layer  and  on  the  Mach  number.  For  example,  G.  Young  and 
E.  Hansen  proposed  the  following  for  the  laminar  boundary  layer  in 
the  Mach-number  region  from  0  to  5  [35] 

'n-  +  o.rxs  (f.,  - 1\)  +  0. 1 9  (f,  -  7j). 

while  for  Mach  numbers  from  5  to  10 

T*  0,7/s  -I  O.SSr^  +  0,2:57sM7 

E.  Eckert  showed  that  it  is  possible  to  give  a  general  formula 
for  the  determining  temperature  (1.6)  for  the  entire  practical  range 
of  Mach  numbers.  Eckert's  verification  of  this  formula  for  a  laminar 
boundary  layer  by  means  of  the  results  based  on  an  exact  calculation 
without  the  determining  temperature  and  obtained  in  a  computer  by 
Young  and  Hansen  Indicated  that  in  the  range  of  Mach  numbers  from 
0.2  to  22  Formula  (1.6)  leads  to  an  error  not  exceeding  3/^  in  the 
determination  of  the  frictional  resistance,  and  consequently,  in  the 
determination  of  heat  transfer;  moreover,  in  the  majority  of  cases 
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(in  38  of  50)  the  error  was  less  than  Vfo  [35]. 

E.  Eckert  [35]  In  analyzing  the  Investigations  and  experimental 
data  produced  by  other  authors  comes  to  the  conclusion  that  Formula 
(1.6)  may  also  be  used  In  the  case  of  a  turbulent  boundary  layer. 

P.  Monaghan  [48]  came  to  the  same  conclusion.  Monaghan  found  on  the 
basis  of  theoretical  Investigations  and  an  analysis  of  the  experi¬ 
mental  data  obtained  with  Mach  numbers  up  to  M  =  8  that  the  deter¬ 
mining  temperature  for  a  turbulent  boundary  layer  may  be  the  same 
as  that  taken  for  a  laminar  boundary  layer. 

K.  Erike  [37]  shows  that  the  determining- temperature  method 
may  be  applied  not  only  to  laminar  and  turbulent  boundary  layers, 
but  also  to  dissociated  gas.  It  is  true  that  he  notes  that  the  ap¬ 
plication  of  the  determining-temperature  method  is  most  convenient 
for  preliminary  aerothermodynamic  calculations.  A  more  exact  solution 
could  be  obtained  with  the  aid  of  computers. 

As  can  be  seen  from  Formula  (1.6),  the  determining  temperature 
may  be  calculated  when  the  recovery  temperature  Tj,  and  the  tem¬ 
perature  of  the  wall  are  known,  but  these  temperatures  are  un¬ 
known.  Consequently,  the  aerodynamic  heating  in  the  case  of  steady- 
state  thermal  processes  should  be  calculated  by  the  method  of  succes¬ 
sive  approximations,  initially  determining  the  Prandtl  number  ap¬ 
proximately.  In  the  case  of  nonsteady- state  thermal  process,  as  is 
shown  in  §  4.1,  the  calculation  may  be  carried  out  without  succes¬ 
sive  approximations. 

The  dependences  of  the  Prandtl  number  on  the  temperature  accord¬ 
ing  to  N.B.  Vargaftlk  [ 25]  and  according  to  E.  Van  Dreist  [5]  are 
given  in  Fig.  1.1.  It  is  apparent  from  this  figure  that  in  the  tem¬ 
perature  range  of  from  0  to  l400°C,  the  Prandtl  number  changes  within  a 
relatively  narrow  range  from  0.68  to  0.72.  Taking  the  mean  value 
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0.70,  In  accordance  v;ith  Formulas  (l.4)  and  (1.5),  we  obtain  for  the 
calculation  of  the  first  approximation 

r ,  ---  0,83;  r,  -  •  0,89. 


The  specific  flow  of  heat  to  the  wall  from  the  air  heated  In 
the  boundary  layer  will  be  expressed  according  to  the  formula  of 
Newton 


(1.7) 


where  a  Is  the  heat-transfer  coefficient. 

The  values  of  the  specific  flows  of  heat  will  be  determined  In 

O 

kcal/m  sec.  Consequently,  the  units  of  the  heat-transfer  coefficient 

O 

will  be  expressed  In  kcal/m  sec 'deg. 

In  addition  to  the  Indicated  gas-kinetic  heat  flow  which  heats 
the  body,  the  outside  siorface  of  the  body  may  be  subject  in  flight 
to  the  influence  of  radiant  energy  emitted  by  the  body,  the  earth, 
the  moon,  comets,  and  by  the  galaxy;  only  flows  of  heat  coming  from 
the  sun  and  the  earth  are  of  practical  significance  In  flights  about 
the  earth.  The  remaining  heat  flows  are  significant  only  in  flights 
In  the  vicinity  of  corresponding  planets  or  in  galactic  flights. 


Fig.  1.1.  Dependence  of  the  Prandtl  number 
on  the  temperature  for  air.  l)  According  to 
N.B.  Vargaftilq  2)  According  to  E.  Van  Drelst. 

The  specific  heat  flow  of  direct  solar  radiation  absorbed  by 
the  irradiated  surface 
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'li'cos-V. 


(1.8) 

where  S  is  the  specific  heat  flow  of  solar  radiation,  normal  to  the 
surface;  f  is  the  angle  between  the  direction  of  the  solar  rays  and 
the  normal  to  the  surface  of  the  body;  ^  is  the  absorption  coeffi¬ 
cient  which  depends  on  the  material  of  the  surface,  its  structure, 
and  temperature. 

The  specific  heat  flow  of  solar  radiation  without  taking  into 
consideration  the  absorption  and  scattering  of  its  atmosphere  in  an 

p 

earth  orbit  (see  [3^])  is  S  =  0.332  kcal/m  sec.  At  altitudes  below 
40  km,  S  will  be  less  than  the  indicated  value  and  will  depend  on 
the  altitude  and  zenith  distance  of  the  sun,  and  also  on  the  meteoro¬ 
logical  conditions.  At  the  surface  when  the  sun  is  at  its  zenith, 

S  decreases  by  approximately  2  1/2. 

The  specific  heat  flow  reflected  from  the  surface  of  the  earth 
and  from  the  clouds  of  solar  rays  is  several  times  less  than  the 
specific  heat  flow  of  direct  solar  radiation.  At  an  altitude  of  500  km 
of  the  earth- sun  line,  in  approximate  terms 

'/or  =  0.016(l+9coso)j}.  (1.9) 

The  specific  heat  flow  of  the  earth's  radiation  will  be  even 
less:  at  a  height  of  500  km  it  is  approximately 

l-2cos?)3.  (1.10) 

The  heat  flow  of  the  reflected  solar  rays  of  the  earth's  radia¬ 
tion  decreases  with  altitude.  At  altitudes  of  several  thousand  kilo¬ 
meters  they  may  be  disregarded. 

The  total  specific  heat  flow  absorbed  by  the  body  in  flight  will 
be 

r  -<i,  I--/.,, 

where 


(1.11) 

(1.12) 
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The  heat  being  absorbed  by  the  body  will  Influence  the  radia¬ 
tion;  It  will  heat  the  body  and  will  be  drawn  beyond  the  limits  of 
the  section  under  consideration  as  a  result  of  heat  conduction. 

The  over-all  heat  balance  when  there  Is  no  cooling  will  be: 

*/  liar  “  "  113  Viiar  (l.l3) 

The  specific  radiation-heat  flow 

(1.14) 

where  a  is  the  Stefan-Boltzmann  constant  and  is  equal  to 

/  2  4 

<f-l3,G-10'=  cal/m  sec -deg  ; 

e  is  the  radiation  coefficient  which  depends  on  the  material  of  the 
surface,  its  structure  and  temperature.  Information  on  the  numerical 
values  of  e  is  given  in  Chapter  2. 

The  specific  heat  flows  involved  in  the  heating  of  the  body 
and  withdrawn  through  the  boundaries  on  the  side  depend  on  the  tem¬ 
perature  gradient  and  the  thickness  of  the  shielding  at  the  surface 
of  the  body  (^T/dy)^^  and  on  the  thermal  conductivity  X: 

where  y  is  the  coordinate  directed  from  the  surface  to  the  flow. 

The  thermal  conductivity  is  seldom  measured  in  kcal/m*hr -deg, 
however,  the  introduction  of  the  hours  as  the  unit  of  time  does 
does  not  correspond  to  the  technical  system  of  units  and  may  be 
Justified  only  in  the  case  of  constructions  in  which  the  thermal 
processes  last  for  hours  and  days.  In  the  case  of  flying  craft  sub¬ 
ject  to  substantial  aerodynamic  heating,  the  length  of  the  limiting 
thermial  processes  is  usually  measured  in  seconds  and  minutes.  Conse¬ 
quently,  the  coefficient  of  thermal  conductivity  in  the  given  work 
will  be  measured  in  kcal/m*sec -deg. 

The  temperature  gradient  at  the  surface  of  the  body  {^T/^y)g^ 
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in  the  general  case  should  be  determined  by  the  Fourier  equation  [g] 

where  £  is  the  specific  heat  capacity  of  the  shielding  material;  y 
is  its  bulk  weight. 

The  specific  heat  flow  involved  In  the  heating  of  the  shield¬ 
ing 


f  ‘1  -Z 


Ox 


should  also  be  found  in  the  general  case  with  Eq.  (l.l6). 

When  there  is  a  temperature  gradient  St/Sx  or  3t/3z,  there  will 
arise  in  the  shielding  heat  flows  perpendicular  to  the  heat  flows 
described  above.  They  will  be  proportional  to  X(3t/3x)  or  X(dT/^z) 
and  should  also  be  determined  in  the  general  case  with  Eq.  (I.l6). 

§  1.2.  HEAT- TRANSFER  COEFFICIENT  WITH  FLOW  ABOUT  A  FLAT  PLATE  AT 
ZERO  ANGLE  OF  ATTACK 


Determination  of  the  heat-transfer  coefficient  at  supersonic 
speeds  is  a  rather  complex  problem  and  at  the  present  time  we  have 
comparatively  exact  solutions  only  for  a  flat  plate  at  zero  angle 
of  attack,  and  that  without  taking  into  consideration  the  pressure 
produced  along  the  plate  by  a  boundary  layer  whose  thickness  is  in¬ 
creasing.  Problems  with  respect  to  determination  of  heat  transfer 
in  the  case  of  more  complex  flows  can  usually  be  reduced  to  those 
of  heat  transfer  for  flat  plates,  making  one  or  another  assumptions. 

The  heat-transfer  coefficient  may  be  expressed  in  terms  of  the 
dimensionless  Stanton  St  number.  The  Stanton  number  is  the  similarity 
criterion  for  heat  transfer  with  forced  motion.  Here 

(l-l?) 

Frequently,  in  addition  to  the  Stanton  number,  we  also  use  the 
dimensionless  Nusselt  number  which  is  equal  to 
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(1.18) 


moreover,  the  Nusselt  number  is  related  to  the  Stanton  number  in 
the  following  way: 


Ku  -  St  Re  Pr, 


(1.19) 


where  Re  is  the  Reynolds  number 


(1.20) 


The  Stanton  number  is  proportional  to  the  local  coefficient  of 
friction 

St-.s-f/.  (1.21) 

moverover,  the  proportionality  coefficient  £  is  a  dimensionless 
quantity.  At  low  velocities  and  small  temperature  differences  in 
the  boundary  layer,  £  in  the  laminar  boundary  as  well  as  in  the 
turbulent  boundary  layer  is  equal  to  (see  [35]) 

.5  .ip,-  ■/.  (1.22) 

2 

To  calculate  the  compressibility  of  the  flow  in  the  case  of 
laminar  boundary  layers,  the  Prandtl  number  in  Formula  (1.22)  should 
be  calculated  for  the  determining  temperature  [35]>  l.e., 

s,s:s;i{Pr*)  (1.23) 


In  the  case  of  a  turbulent  boundary  layer,  the  use  of  Formula 
(1.23)  is  less  well  founded',  however,  until  other  methods  are  developed 
we  may  assume  [35] 


(1.24) 

It  should  be  noted  that  as  a  result  of  the  fact  that  the  Prandtl 
number  varies  within  a  very  narrow  range,  the  range  of  the  possible 
changes  in  the  coefficient  £  will  also  be  narrow.  Adopting  as  the 
Prandtl  number  the  new  value  Pr  =  0.7,  we  will  obtain: 


6- '<.0,64. 
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This  value  may  be  used  for  the  first-approximation  calculation.  If  the 
coefficient  £  is  known,  the  calculation  of  the  Stanton  number  is 
reduced  to  the  determination  of  the  local  coefficient  of  friction. 

For  incompressible  flow  in  the  case  of  a  laminar  boundary  layer, 
the  local  coefficient  of  friction  may  be  expressed  ([35]^  [33]> 
p.  119)  in  the  following  way: 


(1.25) 


In  the  case  of  a  turbulent  boundary  layer,  we  may  use  the 
following  formula  up  to  Re  =  lo"^  (see  [48],  [33],  p.  425): 


0,1  592 


(1.26) 


The  formula  is  in  good  agreement  with  experiments  for  values 
10^  <  Re  <  10^  (see  [35],  [33],  p.  433) 

c,^-= (1.27) 

The  above  formulas  may  also  be  used  in  the  case  of  compressible 
flow,  if  the  Reynolds  number  is  calculated  from  the  determining  tem¬ 
perature,  using  Formula  (1.6).  Here  the  coefficient  c^  is  referred  to 
the  physical  parameters  of  the  air  for  the  determining  temperature, 
i.e.,  to  p*  and  p.*.  The  Stanton  number  in  Formula  (I.17)  is  referred  to 
Cpgj  for  compressible  flow  it  is  necessary  to  calculate  c^  or  the 
determining  temperature.  If  we  consider  the  above,  the 
friction  coefficient  referred  to  the  physical  properties  of  the  air 


at  the  limit  of  the  boundary  layer,  will  be: 


(•  -r’  <’* 


(1.28) 


while  the  Stanton  member 


Substituting  this  expression  in  Formula  (I.17),  we  obtain  for 


the  local  heat-transfer  coefficient: 
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a  ^  2 . 29) 

If  we  take  into  consideration  Formulas  (1.23),  (1.24),  (1.25), 
(1.26),  and  (1.27),  we  will  find  the  over-all  heat- transfer  coeffic¬ 
ient;  here  the  heat-transfer  coefficient  will  be: 

for  the  laminar  boundary  layer 

='..-=;i26(kc^^)  (1.30) 

for  the  turbulent  boundary  layer  and  Re  <  10^ 

a,-0,29(Re'^)-''’--  (Pr^)-  ( 1 . 31) 

6  Q 

for  the  turbulent  boundary  and  10  <  Re  <  10^ 

1 ,81  (ig  (1.32) 

The  value  of  the  density  for  the  detennining  temperature  may 
be  determined  from  the  characteristic  equation,  taking  into  con¬ 
sideration  the  fact  that  the  pressure  in  the  boundary  layer  does 
not  change: 

(1.33) 

The  Prandtl  number  in  these  formulas  may  be  determined  from  the 
graph  given  in  Fig.  1.1.  To  determine  the  Reynolds  number,  it  is 
necessary  to  know  the  relation  of  the  viscosity  to  the  temperature. 
For  this  we  may  use  the  formula  (see  [48]) 

P  -- 1 ,49 •  10"'  -  kg  s ec/ m^ .  ( 1  •  34 ) 

This  formula  is  thoroughly  validated  by  experiments  in  the 
temperature  range  of  100  to  2000°K.  The  coefficient  of  kinematic 
viscosity  in  the  Reynolds  number 

>=«*/?.  (1.35) 

The  specific  heat  capacity  of  the  air  may  be  determined  from 
the  graph  given  In  Fig.  1.2  This  graph  is  constructed  from  the  data 
of  N.V.  Vargaftlk  [25]  and  E.  Eckert  [35]. 

For  flying  craft  with  cosmic  and  near-cosmic  velocities  the 
thermal  regime  may  be  of  the  essence  at  altitudes  where  due  to  the 
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great  rarefaction  of  the  air,  boundary  layer  theory  is  not  applicable. 
This  is  the  region  of  free-molecular  flow;  It  Is  characterized 
by  the  fact  that  the  mean  molecular  free  path  (before  collision)  be¬ 
comes  greater  than  the  dimensions  of  the  body.  In  such  flow,  there 
Is  no  boundary  layer.  The  relation  of  the  mean  molecular  free  path 
to  the  linear  dimension  of  the  body  is  called  the  Knudsen  number. 

The  region  of  the  free-molecular  flow  Is  characterized  by  the  Knudsen 
number  (see  [50]) 

Kn= 


OfiP 
0^8 
0.2S 
0,U 

0  200  <t00  SOO  800  WOO  1200  IMOOtX 

Pig.  1.2.  Dependence  of  the  specific 
heat  capacity  of  the  air  on  the  tem¬ 
perature.  l)  According  to  E.  Eckert; 

2)  According  to  N.V.  Vargaftlk. 

The  heat  transfer  of  the  free-molecular  flow  depends  on  the 
relative  amount  of  energy  transmitted  by  the  molecules  to  the  air 
against  which  they  collide.  This  relative  energy  is  called  the  accom¬ 
odation  coefficient,  the  value  of  which  Is  Influenced  to  a  certain 
extent  by  the  material  and  the  structure  of  the  surface.  The  accom¬ 
odation  coefficient  values  are  close  to  unity;  Table  1.1  gives  the 
value  of  the  accomodation  factor  for  certain  materials  [39].  A  good 
agreement  between  the  theoretical  heat  transfer  and  the  experimental 
heat  transfer  Is  obtained  If  we  assume  A  =  0.9  [23]. 
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TABLE  1.1 

Accomodation  Coefficients  for  Certain  Materials 


MarepiMJi 

1 

UtieHT  aKKO- 
MOAaUHH 

2  A 

MarepHaa 

Koa^M- 
UHCNT  aKKO- 
MOAaitHH 

A 

A.IIOMHHHft  MaUIHIIUOtl 
Cpa60TKH  3 

0.95-i-0.97 

JlMTaa  CTaab  MauiHHHoit 
o6pa6oTKH  5 

0.87^0.88 

AaioMMiiHil  Tpa&aciiuil4 

0,89•^0.97 

Y  JliiTaa  crajib  TpaBjieHaa 

0,69-!-0.96 

AaiOMHHHt}  iiOAiipoaaii* 

0.87-^0.95 

JliiTaa  craab  n^iipo* 

BaHnaa  8 

« 

0»87*4.0,93 

l)  Material;  2)  accomodation  coefficient  A;  3)  machine- 
worked  aluminum;  4)  etched  aluminum;  5)  polished 
aluminum;  6)  machine- worked  lithium  steel;  7)  etched 
lithium  steel;  8)  polished  lithium  steel. 

The  heat  transfer  of  the  free-molecular  flow  may  also  be  deter¬ 
mined  from  Formulas  (l.l)  and  (1.7);  we  must,  however,  take  into 
consideration  the  fact  that  the  temperature  recovery  factor  in  free- 
molecular  flow  increases  and  becomes  greater  than  unity.  Theoreti¬ 
cally,  in  the  case  of  free-molecular  flow  past  a  flat  plate,  the 
recovery  factor  equals  (see  [23]) 

/•,.«  7/()  1,17. 

Experiments  confirm  the  theoretical  value  of  the  recovery 
factor  for  the  region  4  >  Kn  ^  Re  >  1.  For  Kn  ^  Re  >  4,  the  experi¬ 
mental  values  are  higher  than  the  theoretical  values  and  reach  r^^^  = 
=  1.4  (see  [23]). 

The  heat-transfer  coefficient  may  be  determined  from  Formula 
(1.17),  while  the  Stanton  nxmiber  for  free-molecular  flow  past  a  flat 
plate  at  zero  angle  of  attack  will  be  (see  [23]): 

St,„:-0,242A/:.  (1.37) 

where  A  is  the  accomodation  coefficient,  while  ^  is  the  ratio  of  the 
velocity  of  light  to  the  most  probable  molecular  velocity.  The  value 
of  C  may  be  expressed  in  terms  of  the  Mach  number  (see  [23]): 

(1.38) 
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Substituting  Formulas  ( 1.137)  and  (1.38)  in  Formula  (1.17),  we 
obtain: 

(1.39) 

It  Is  apparent  from  this  formula  that  the  velocity  of  the  flight 
has  no  Influence  on  the  heat-transfer  coefficient  In  free-molecular 
flow  while  If  we  take  Into  consideration  the  small  changes  in 
a,  ic,  c  ,  and  A,  the  density  of  the  air  Is  the  basic  factor  deter- 
mining  the  heat-transfer  coefficient. 

In  free-molecular  flow  there  will  be  no  shock  wave;  however, 
due  to  the  collisions  of  the  molecules  with  the  surface  these  mole¬ 
cules  may  be  dissociated  and  ionized.  Due  to  the  great  length  of 
the  molecular  free  path  in  free-molecular  flow,  at  great  hyper¬ 
sonic  velocities  an  extensive  ionization  zone  developes  in  front 
of  the  moving  body. 

For  the  range  of  Knudsen  numbers  [50] 

0,01<Kn<0,l 

the  flow  will  have  a  boundary  layer;  however,  the  velocity  at  the 
surface  will  not  equal  zero.  This  is  This  is  called  slippage  flow. 

The  sudden  increase  in  the  velocity  in  the  boundary  layer  with 
slippage  also  produces  a  sudden  increase  in  the  temperature.  The 
theoretical  determination  of  heat  transfer  in  the  case  of  flow  with 
slippage  is  quite  complex  and  extremely  unreliable  [50].  A  compari¬ 
son  of  experimental  data  on  heat  transfer  with  data  calculated  on 
the  basis  of  the  theory  of  constant  flow  (continuum)  with  a  laminar 
boundary  layer  indicates  that  the  heat  transfer  in  the  case  of  flow 
with  slippage  is  somewhat  lower  than  the  heat  transfer  obtained  from 
formulas  for  continuum  flow  [50]. 
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§  1.3.  THE  HEAT- TRANSFER  COEFFICIENT  IN  THE  CASE  OF  FLOW  PAST  A 
PLATE,  A  WING  PROFILE,  AND  A  BODY  OF  REVOLUTION 

At  small  angles  of  attack  for  a  flat  plate,  the  heat  transfer 
can  be  determined  from  formulas  corresponding  to  the  plate  with  a 
zero  angle  of  attack,  taking  Into  consideration  the  local  values 
of  the  velocity,  the  density,  and  the  temperature  of  the  air  at  the 
limit  of  the  boundary  layer.  It  Is  the  lower  surface  that  Is  of 
greatest  significance  from  the  standpoint  of  heat  shielding,  since 
the  highest  temperatures  are  attained  on  this  surface. 


Fig.  1.3.  Diagram  of  flow  about  the  nose 
of  the  profile,  l)  Shock  wave. 

In  the  case  of  flow  about  the  lower  surface  of  a  plate  with  a 
positive  angle  of  Incidence,  an  oblique  compression  wave  develops  on 
the  leading  edge.  Let  us  consider  the  more  general  case  of  a  wing 
profile,  when  the  tangent  to  the  profile  at  the  Investigated  point 
A  forms  the  angle  cp  with  the  direction  of  the  flow  on  the  lower 
surface  (see  Fig.  1.3).  In  this  case  the  local  pressure  p^  at  point 
A  on  the  surface  of  the  wing  may  be  represented  according  to  Buze- 
man  In  the  following  way  (see  [29]: 

.  .  .  (1-^0) 

We  are  frequently  restricted  to  two  tenns  of  this  expansion. 
This  Is  entirely  permissible.  If  the  Mach  numbers  are  not  substan¬ 
tially  greater  than  5-6.  At  great  hypersonic  velocities  the  neglect¬ 
ing  of  the  third  term  may  produce  a  substantial  error.  As  an  exartple. 
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Fig.  1.4  gives  the  error  as  a  function  of  M  when  cp  =  5*^.  As  can 
be  seen  from  the  curve,  when  >  20,  the  error  may  be  greater  than 
20^.  Consequently  we  Introduce  the  values  of  the  factors  for  a 
trinomial  expansion. 

If  we  assume  that  the  adiabatic  exponent  for  air  is  1.4,  for 
the  coefficients  c^^,  Cg,  and  c^  we  obtain  the  following  expression, 
if  cp  is  determined  in  radians  (see  [29]): 


I  [(Ml  -  2)H  1 .4Ml  I  (Ml  - 1  )-*: 
f,= (0.36MI  - 1 ,493Ml + 3,6M1  - 2MI  + 1 .33) (Ml  - 1)"’’*. 


(1.41) 

(1.42) 

(1.43) 


Here  the  coefficient  c^  reflects  the  Influence  of  the  shock 
wave.  When  the  flow  expands  (when  there  is  no  shock  wave),  the 
values  of  the  coefficients  c^^  and  Cg  will  remain  the  same,  while 
coefficient  c^  is  somewhat  changed  and  can  be  expressed  in  the  follow¬ 
ing  way; 

<?i=(0,4Ml-l.«13Ml  +  4Ml-2Ml  +  i.33)(Ml-l)-”.  (1 .44) 

In  the  case  of  flow  expansion,  cp  is  negative.  If  the  flow  ex¬ 
pansion  precedes  the  shock  wave,  as  for  exanple  in  the  case  of  flow 
about  the  rear  section  of  the  profile,  the  coefficient  c^  is  also 
calculated  from  Formula  (1.44),  but  a  term  expressing  the  influence 
of  the  shock  wave  is  added  to  the  expansion: 

where  9^^  is  the  angle  of  inclination  of  the  surface  in  the  forward 
section  In  front  of  the  shock  wave,  and 

fJ-Cj-c;-  ~(0,04Ml-0.32Ml+0.4Ml)(Ml-l)-*-*.  ( 1 .45) 

Figure  1.5  presents  a  graph  for  the  determination  of  the  coef¬ 
ficients  Cj^,  Cg,  c^,  and  c^  as  functions  of  the  number  In  undis¬ 
turbed  flow. 
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Having  determined  the  pressure  on  the  lower  surface  of  the  plate 
or  on  the  nose  of  the  wing  profile  by  means  of  Formula  (1.40),  we 
can  determine  the  air  density  at  this  same  surface  with  the  Hugonlot 
equation  (see  [29]): 


f  •  PtIPm  +  ® 


(1.46) 


Fig  1.4.  Relative  error  if  we  neglect 
the  third  term  of  the  expansion  in 
series  of  the  relative  pressure  values. 


If  there  is  no  shock  wave  (for  example,  for  the  upper  surface 
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of  a  plate  or,  In  certain  cases  for  the  profile),  the  density  must 
be  determined  from  the  equation  of  the  adiabatic  curve: 

p*/p. (ptipy' « (1.^7) 

If  the  flow  expands  behind  the  shock  wave,  we  Initially  deter¬ 
mine  the  density  behind  the  shock  wave  from  Formula  (1.46),  and 
then  by  means  of  the  equation  of  the  adiabatic  curve  we  determine 
the  density  at  the  surface  being  Investigated. 

The  air  temperature  at  the  limit  of  the  boundary  layer  will  be 
found  from  the  characteristic  equation: 

r./r.=(p,//»j(pjp.).  (1.48) 

Making  use  of  the  fact  that  the  flow-retardation  temperature 
does  not  change  even  In  the  presence  of  compression  waves,  we  ob¬ 
tain: 

Taking  into  consideration  that  k  =  1.4,  we  find: 

mI=5  l(r-/r»)(i  +o.2Mi)-ii.  (1*^9) 

Formula  (1.40)  and  the  values  of  the  coefficients  In  this  for¬ 
mula  correspond  to  a  wing  of  Infinite  span.  However,  this  formula 
can  also  be  used  in  the  determination  of  the  pressure  on  a  straight 
wing  of  finite  span  for  that  part  located  outside  the  Mach  end  cones. 

At  the  ends  of  the  wing  the  pressure  drops  linearly  along  the 
span;  moreover,  at  the  very  end,  the  excess  pressure  equals  zero. 
Hence,  It  Is  not  difficult  to  determine  the  pressure  at  any  point 
Inside  the  Mach  end  cones. 

In  the  case  of  aerodynamic  heating  the  case  in  which  the  lead¬ 
ing  edge  of  the  delta  wing  becomes  "supersonic”  is  of  particular  In¬ 
terest;  In  this  case  the  excess  pressure  is  greater  than  on  a  wing 
with  a  "subsonic"  leading  edge. 
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For  a  very  thin  wing  (a  plate) ,  the  excess  pressure  on  the  sec¬ 
tion  lying  outside  of  the  Mach  cone  emanating  from  the  peak  of  the 
triangle  forming  the  wing  will  be  (see  [26],  [2?]): 


where  x  is  the  sweepback  angle  of  the  leading  edge;  for  the  section 
of  the  wing  lying  inside  the  Mach  cone,  the  pressure  will  be  less 


and  expressed  as: 


P.  r  1  -1  arc  sin  .'n  ]/ 

KM.-I  — tg*x  "  /M»— IK  1— 


where  x  is  the  coordinate  along  the  flow,  while  z  is  the  coordinate 
perpendicular  to  the  plane  of  symmetry  (the  origin  of  the  coordi¬ 
nates  is  at  the  apex  of  the  triangle) . 

For  a  wing  of  infinite  span,  with  a  very  thin  profile  (a  plate), 
the  excess  pressure  according  to  the  linearized  theory  is 


The  greatest  heat  flows  are  developed  on  a  delta  wing  in  the 
section  lying  outside  of  the  Mach  cone.  The  pressure  on  this  sec¬ 
tion  of  the  wing  may  be  determined  from  the  formula 


;  m-i  - 

while  calculating  the  value  of  Pj^^  from  Formula  (1.40). 

It  should  be  noted  that  as  the  Mach  number  increases,  the  cor¬ 
rection  factor  approaches  unity. 

Thus,  when  M  =  5  and  x  =  60°,  this  factor  equals  1.07,  while 
when  M  =  10,  it  equals  I.015.  The  latter  formula  may  also  be  used 
for  the  determination  of  the  pressure  on  a  sweptback  wing,  l.e.,  on 
that  section  of  it  which  lies  outside  the  central  Mach  cone. 

The  formulas  presented  here  together  with  the  formulas  given 
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In  §  1.2  enable  us  to  find  the  heat-transfer  coefficient  on  a  flat 
plate  as  well  as  on  a  thin  wing,  If  the  angles  of  Incidence  are  rela¬ 
tively  small. 

In  §  3.1  we  give  an  example  of  the  determination  of  the  heat- 
transfer  coefficient  on  the  profile  of  a  wing. 

It  should  be  noted  that  the  presence  of  pressure  and  tempera¬ 
ture  gradients  on  the  surface  (for  example,  on  a  curvilinear  profile) 
increases  the  error  In  the  determination  of  the  heat-transf er  coef¬ 
ficient;  moreover,  this  error  Increases  as  the  Mach  number  Increases. 

The  heat  transfer  on  a  cylindrical  surface  with  an  axis  parallel 
to  the  flow  may  be  determined  from  the  formulas  for  a  flat  plate. 

If  the  thickness  of  the  boundary  layer  is  substantially  less  than 
the  curvature  radius  of  the  surface.  The  applicability  criterion  for 
the  theory  of  a  flat  plate  for  a  cylindrical  surface  may  be  expressed 
in  the  following  inequalities: 

for  a  laminar  boundary  layer  x/R  <  0.02  Re  ;  (l.50); 

for  a  turbulent  boundary  layer  x/R  <  0.3  Re^'^^,  (l.5l). 

The  heat -transfer  formulas  for  a  flat  plate  may  also  be  used 
for  a  cone;  however,  the  local  values  of  M^,  Pg,  and  T^  should  corre¬ 
spond  to  the  flow  past  the  cone  and  the  conical  nature  of  the  flow 
should  be  taken  into  consideration. 

Because  of  the  conical  surface  not  all  of  the  gas  streams 
at  the  surface  have  identical  lengths  of  contact  with  the  surface. 

The  streams  beginning  at  the  point  of  the  nose  have  the  greatest 
lengths  of  contact  with  the  surface.  Other  streams  are  shorter,  since 
they  approach  the  surface  at  some  distance  from  the  point  of  the 
nose.  Thus,  while  determining  the  heat  transfer  or  the  frictional 
coefficient  at  a  given  point  of  the  cone,  it  should  be  borne  in  mind 
that  the  mean  length  at  which  particles  of  gas  in  the  boundary  layer 
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are  subject  to  the  Influence  of  viscosity  forces  will  be  less  than 
the  length  of  the  generatrix. 

Consequently,  in  using  the  heat-transfer  formulas  for  a  flat 
plate  in  the  determination  of  the  heat  transfer  on  a  cone,  we  should 
find  the  Reynolds  number  from  the  given  length  which  depends  on 
the  length  of  the  generatrix  x  and  the  structure  of  the  boundary 
layer.  In  the  case  of  laminar  boundary  layers  Xp^^/x  =  l/3;  in  the 
case  of  turbulent  boundary  layers  x^^/x  =  l/2  (see  [5],  [30]). 

To  determine  the  parameters  of  the  flow  at  the  surface  of  the 
cone,  let  us  Introduce  the  similarity  parameter  for  conical  flows 

E  =  (1.52) 

where  is  in  radians . 

Using  the  denotations  of  Fig.  1.6,  according  to  N.P.  Krasnov 

([10],  p.  258),  we  will  have: 

2,09(1  +0.143/5’'’)»il.  (1.53) 

where 

-  Fi-P^ 

P*^-. - • 

(1.54) 

To  determine  the  remaining  local 
parameters  of  the  flow  at  the  limit 
of  the  boundary  layer,  it  is  neces¬ 
sary  to  take  into  consideration  the 
nonisotroplc  change  in  the  flow  with 
the  transition  through  the  compres¬ 
sion  wave;  these  changes  depend  on 

Pig.  1.6.  Diagram  of  the  sym-  the  angle  of  inclination  0  of  the 
bols  of  the  parameters  of  the  ° 

flow  about  the  cone,  l)  Com-  surface  of  the  compression  wave  and 
pression  wave;  2)  cone. 

The  angle  of  inclination  of  the 
compression  wave  according  to  N.P.  Krasnov  ([10],  p.  249),  is  deter- 


-  28  - 


mined  from  the  following  formula 

1.093  +  0.06  (1.55) 

Formula  (1.53)  and  (1.55)  are  approximations  of  the  exact  solu¬ 
tion  and  may  be  used  with  a  sufficiently  high  accuracy  for  0.5  <  4  <  “j 
>  2.5,  and  \  <  30°. 

If  we  know  the  angle  of  inclination  for  the  compression  wave, 
we  can  determine  the  relative  pressure  behind  the  discontinuity  ([10], 
P.  37): 


D 


(1.56) 


The  density  of  the  air  behind  the  compression  wave  is  found  from 
the  Hugonlot  equation 


_  =1+5---  *• 


(1.57) 


The  transition  to  the  density  of  the  air  at  the  surface  of  the 
cone  can  be  accomplished  by  means  of  the  equation  of  the  adiabatic 
curve : 


P,/Pc= (W7?c)''*=f(  J  (1.58) 

The  quantity  Pg/p^^  in  this  equation  can  be  determined  from 
Formula  (1.54),  using  elementary  transformations  of  this  value;  here 
we  will  obtain: 

(1.59) 


The  temperature  Tg  at  the  limit  of  the  boundary  layer  of  the 
cone  and  the  local  Mg  number  can  be  determined  from  Formulas  (1.48) 
and  (1.49),  respectively. 

In  the  case  of  a  pointed  body  of  revolution  with  a  curvilinear 
generatrix,  the  flow  parameters  at  the  point  may  be  determined  from 
the  formulas  for  a  cone  having  an  angle  (see  Pig.  1.7).  At  other 
points  on  the  nose  surface,  the  flow  parameters  can  also  be  found 
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approximately  from  the  formulas  for 
a  cone  with  an  angle  which  Is 
equal  to  the  angle  of  Inclination 
of  the  nose  surface  at  the  given 
point  (see  Fig  1.7). 

§  1.4.  HEAT- TRANSFER  COEFFICIENTS  IN 
THE  CASE  OF  FLOW  ABOUT  FRONTAL 
AREAS 

The  forward  noses  of  bodies  and 
the  leading  edges  of  the  lifting  sur¬ 
faces  are  related  to  the  frontal  sur¬ 
faces  .  Heat  transfer  In  the  case  of 
flow  about  frontal  surfaces  may  be 
determined  from  the  formulas  for  a  plate.  If  we  take  Into  considera¬ 
tion  the  local  characteristics  of  the  flow. 

With  flow  about  rounded  noses  or  leading  edges,  the  blunt  sec¬ 
tion  of  the  nose  in  the  region  of  the  critical  point  is  streamlined 
by  subsonic  flow  (Fig.  1.8).  On  the  basis  of  experiments  [5]  it  was 
determined  that  the  velocity  of  the  flow  along  the  meridional  cross 
section  is  proportional  to  the  relative  distance  from  the  critical 
point : 

V,=>xV^xlR.  (1.60) 

The  geometric  values  of  x  and  R  are  given  in  Fig.  1.8.  The  value 
of  the  coefficient  x  depends  on  the  shape  of  the  frontal  surface 
and  the  flight  Mach  number.  In  the  case  of  Incompressible  flow 
(Mjjp  =  O),  X  =  1«5  for  a  spherical  nose,  and  for  a  cylindrical  lead¬ 
ing  edge  X  =  2;  when  increases,  the  value  of  x  decreases  and 
may  be  determined  for  a  spherical  nose  from  the  formula  of  M. 

Romlg  [ 20 ] : 

x=o.8m:;®-“.  (i.6i) 


Fig.  1.7.  Diagram  of  denota¬ 
tions  for  the  geometric  char¬ 
acteristics  of  a  pointed  body 
of  revolution.  1)  Compression 
wave . 
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This  formula  Is  used  for  M  >1. 

00 

Because  of  the  large  negative  pressure  gradient  and  the  low 
Reynolds  numbers,  the  boundary  layer  must  be  laminar  In  the  region 
of  the  critical  frontal- surface  point.  The  Reynolds  number  for  the 
determining  temperature.  If  we  take  Into  account  Formula  (1.60), 
will  be 

Re*  =  JcKi/v* = X  V'-JC*/v*/?. 

Substituting  this  value  of  the  Reynolds  number  Into  Formula 
(1.30),  and  also  substituting  the  value  of  from  Formula  (I.60), 
we  will  obtain  the  expression  for  the  heat-transfer  coefficient  In 
the  region  of  the  critical  point  of  the  nose; 

«o=i4x''*v':^'/?"''*(i‘*p*)''*(Pr*r’'T;.  (1.62) 

Here  the  constant  A  Is  a  function  of  nose  shape.  Using  the  cor¬ 
responding  numerical  values  given  by  E.  Van  Drelst  [5]>  we  may  as¬ 


sume  the  following: 

for  a  spherical  nose  .  A  =  7.5 

for  a  circular  cylindrical  nose  . A  =  5.6 


The  determining  temperature  at  the  critical  point  may  be  deter¬ 
mined  from  the  formulas  for  the  deceleration  temperature  and  from 
Formula  (1.6);  taking  Into  consideration  the  fact  that  In  the  case 
under  consideration  T^  =  Tq,  where 

we  get; 

r*  -  r,-  0.5  (r,-  rj  -  0,22  (i  -  r)  (t*,-  t-.). 

To  determine  the  air  density  at  the  body  surface  from  the  deter¬ 
mining  temperature,  we  first  determine  the  pressure  behind  the  com¬ 
pression  wave  In  the  region  of  the  critical  point.  From  this  pres¬ 
sure,  we  find  the  density  of  the  air  behind  the  compression  wave 
and  then  from  the  adiabatic  equation  and  the  characteristic  equation 
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we  make  the  transition  to  the  density  of  the  air  In  the  case  of  the 
determining  temperature. 

For  a  spherical  or  straight  cylindrical  nose,  the  pressure  be¬ 
hind  the  compression  wave  may  be  found  from  Formula  (1.56),  If  we 
assume  that  =  90°.  The  density  behind  the  compression  wave  Is 
determined  from  the  Hugoniot  equation  (1.57) •  The  transition  from 
this  density  to  the  density  at  the  surface  can  be  made  by  means 
of  the  adiabatic  equation  (1.58).  The  sought  value  of  the  air  density 
for  the  determining  temperature  may  be  found  from  the  characteristic 
equation: 

p*/pc=(p«/p.)  (p*/p.)=(p»/p.)  (To/7*). 

In  determining  the  heat  transfer  on  the  leading  edge  of  the  wing, 
it  should  be  borne  In  mind  that  the  sweepback  of  the  wing  decreases 
heat  transfer.  In  approximate  terms,  this  decrease  In  the  heat  trans¬ 
fer  may  be  determined  from  Fig.  1.9  (see  [24]).  In  selecting  the 
sweepback,  it  should  be  borne  In  mind  that  a  very  large  sweepback 
(greater  than  70°)  decreases  the  critical  Reynolds  number. 


Fig.  1.8.  Diagram  of  flow  about 
a  round  nose.  1)  Compression 
wave;  2)  zone  of  subsonic  flow. 


2)  tfA 

Fig.  1.9*  Influence  of  the  sweep- 
back  of  the  wing  on  the  heat  trans¬ 
fer  at  the  critical  point  of  the 
leading  edge.  1)  (a  with  sweepback)/ 
/(a  without  sweepback) ;  2)  sweep- 
back,  deg. 
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If  a  rounded  nose  Is  streamlined  by  a  laminar  boundary  layer, 
the  heat  transfer  at  the  critical  point  will  be  at  Its  maximum.  How¬ 
ever,  If  the  boundary  layer  at  the  nose  becomes  turbulent,  the  maxi¬ 
mum  heat  transfer  will  be  In  the  vicinity  of  the  point  with  the  local 
number  =  1  [ 21 ] . 

§  1.5.  FACTORS  INFLUENCING  THE  TRANSITION  FROM  A  LAMINAR  BOUNDARY 
LAYER  TO  A  TURBULENT  BOUNDARY  LAYER 

The  structure  of  the  flow  In  the  boundary  layer  may  be  laminar 
or  turbulent.  In  the  case  of  a  laminar  boundary  layer,  the  frictional 
resistance  Is  substantially  less  than  In  the  case  of  a  turbulent 
boundary  layer.  The  decrease  In  the  friction  of  the  laminar  boundary 
layer  decreases  the  heat  transfer,  and  consequently,  the  temperature 
of  the  surface  also. 

As  an  example  Fig.  1.10  gives  steady-state  temperatures  of  the 
wing  surface  In  the  case  of  laminar  and  turbulent  boundary  layers 
when  Mj^  =  5.  It  Is  apparent  from  the  graph  that  In  the  case  of  the 
laminar  boundary  layer,  the  temperature  of  the  surface  Is  200  to 
235°  lower  than  In  the  case  of  a  turbulent  boundary  layer.  At  higher 
values  of  M^,  the  temperature  drop  across  the  surface  In  the  case 
of  a  laminar  boundary  layer  will  be  even  greater. 

Thermal  shielding  of  the  structure  Is  one  of  the  most  effective 
methods  of  ensuring  a  laminar  structure  for  the  boundary  layer.  The 
laminar Izat ion  of  the  boundary  layer  also  decreases  the  frictional 
resistance  and,  by  the  same  token,  decreases  the  weight  of  the  engine. 

If  there  are  no  factors  producing  turbulence  In  the  nose  sec¬ 
tion  of  a  streamlined  body,  the  nose  section  will  generally  be  stream¬ 
lined  by  a  laminar  boundary  layer.  The  transition  from  a  laminar 
boundary  layer  to  a  turbulent  boundary  layer  at  subsonic  flight  velo¬ 
cities  depends  on  the  Reynolds  number,  and  on  a  flat  plate  with  a 
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zero  angle  of  incidence,  this  transition  generally  takes  place  at 
Rej^r  =  5*10^  [33]. 

At  supersonic  flight  velocities,  as  experiment  shows,  the 
boundary  layer  becomes  turbulent  at  higher  Reynolds  numbers. 

The  cause  of  turbulence  in  a 
boundary  layer  is  the  instability  of 
the  laminar  boundary  layer  at  high 
Reynolds  numbers.  In  connection  with 
this,  the  random  velocity  pulsations 
In  the  boundary  layer  are  not  at¬ 
tenuated  but  are  developed  and  cause 
the  turbulence  of  the  boundary  layer. 
Time  is  required  to  develop  the  turbu¬ 
lence  and  the  disturbances,  and  con¬ 
sequently,  the  Reynolds  number  for 
the  initial  turbulence  is  greater  than  the  Reynolds  number  for  the 
initial  steady-state  motion. 

The  source  of  the  disturbing  pulses  is  usually  the  turbulence 
of  the  exterior  flow,  roughness,  undulation,  and  other  disruptions 
of  smooth  flow  about  the  surface.  By  decreasing  the  disturbing 
pulses,  we  Increase  the  time  required  for  the  development  of  the 
turbulence,  and  consequently.  Increase  the  critical  Reynolds  number 
for  the  transition  from  a  laminar  to  a  turbulent  boundary  layer. 

For  example,  by  decreasing  the  turbulence  of  the  flow  through  the 
use  of  damping  grids  we  succeeded  in  obtaining  a  critical  Reynolds 
number  of  RSj^j,  =  3*10^  in  addition  to  the  usual  value  5 ’10^  [33]  on 
a  flat  plate  at  subsonic  velocities. 

Although  it  might  have  been  assumed  that  the  viscosity  forces 
in  the  boundary  layer  would  be  a  stabilizing  factor,  in  reality  they 
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Pig  1.10.  Temperature  of  a 
wing  surface,  with  aerodynamic 
heating,  l)  Turbulent  boundary 
layer;  2)  laminar  boundary 
layer;  3}  wing  chord;  M  =  5> 

H  =  30  km,  a  =  10°,  e  =  0.8, 
c  =  0.03. 


reduce  the  stability  of  the  laminar  motion.  In  this  a  basic  role 
is  played  by  the  nature  of  the  change  in  the  velocity  through  the 
thickness  of  the  boundary  layer}  in  the  final  analysis  it  is  the 
nature  of  this  change  which  is  determined  by  the  viscosity  forces. 

If  the  wall  cools  the  boundary  layer,  for  example,  by  means  of 
its  heat  capacity  or  radiation,  the  density  of  the  air  at  the  wall 
is  Increased  and,  by  the  same  token,  the  kinetic  energy  referred  to 
the  volume  is  increased.  This  Increases  the  stability  of  the  boundary 
layer . 

Figure  1.11  illustrates  the  nature  of  the  change  in  the  tempera¬ 
ture  and  the  density  of  the  air  through  the  thickness  of  the  boundary 
layer.’  Curves  MN'P'  and  AB'C  correspond  to  the  adiabatic  wall, 
while  curves  MNP  and  ABC  correspond  to  the  wall  which  is  cooling 
the  boundary  layer.  In  the  section  BC,  the  density  is  Increased  and 
this  Increases  the  stability  of  the  boundary  layer. 


Fig.  1.11.  Temperature  distribution  and  air  den¬ 
sity  through  the  thickness  of  the  boundary  layer. 

With  supersonic  velocities  and  real  walls  there  is  always  heat 
transfer  from  the  boundary  layer,  although  this  is  the  result  of  the 
radiation  of  the  wall.  Consequently,  the  critical  Reynolds  numbers 
at  supersonic  velocities  are  substantially  greater  than  the  critical 
Reynolds  numbers  of  the  subsonic  section. 

L.  Lees  [47]  showed  theoretically  that  by  cooling  a  wall  stream- 
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lined  by  flow,  we  can  completely  stabilize  the  boundary  layer  of 
the  supersonic  flow,  i.e.,  the  critical  Reynolds  numbers  may  be 
made  higher  than  those  used  in  practice.  In  these  cases  the  practi¬ 
cal  criterion  of  stability  will  not  be  the  Reynolds  number,  but  the 
relative  wall  temperature  or 

E.  Van  Drelst  calculated  the  required  relative  temperatures  of 
the  wall  for  total  stabilization  of  the  boundary  layer  at  various 
Mach  numbers  [52]. 

Figure  1.12  shows  a  graph  of  the  stability  limits  according  to 
E.  Van  Drelst:  the  regions  within  the  curves  B,  C,  and  D  correspond 
to  the  laminar  boundary  layer  [51.  The  ratio  of  the  ordinates  of  the 
curves  B,  C,  and  D  tothe  ordinate  of  the  curve  A  characterizes  the 
degree  of  required  cooling  for  total  stabilization  of  the  boundary 
layer . 


Pig.  1.12.  Limits  of  the  total  stabili¬ 
zation  of  the  boundary  layer  according 
to  E.  Van  Drelst.  l)  Cone;  2)  plate; 

3)  in  flight;  4)  In  a  tube;  5)  laminar 
boundary  layer. 


The  region  above  the  stability  curves  do  not  necessarily  cor¬ 
respond  to  the  turbulent  boundary  layer.  For  these  regions  the  transi¬ 
tion  from  the  laminar  boundary  layer  to  a  turbulent  boundary  layer 
will  be  determined  by  the  Reynolds  number. 

Notwithstanding  the  great  number  of  theoretical  investigations 
of  boundary- layer  stability,  the  mechanism  of  the  transition  of  a 
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laminar  boundary  layer  Into  a  turbulent  boundary  layer  has  not  been 
completely  discovered  as  yet.  As  a  result,  the  determination  of  the 
critical  Reynolds  numbers  by  theoretical  means  Is  not  yet  possible. 

The  experimental  data  have  not  been 
adequately  systematized;  moreover, 
there  is  a  great  divergence  in  the 
result  of  various  experiments.  This 
is  explained  by  the  various  condi¬ 
tions  under  which  experiments  have 
been  conducted  and,  in  particular,  by 
the  fact  that  the  flow  turbulence, 
the  roughness  of  the  models,  and  the 
pressure  gradient  are  not  considered 
in  many  of  the  experiments. 

The  extent  of  turbulence  can  be  seen  from  the  experiments  on 
the  determination  of  the  critical  Reynolds  numbers  on  a  plate  when 
M  <  1  [12].  In  these  experiments  the  flow  turbulence  was  changed 
and  the  formation  of  a  temperature  gradient  was  excluded.  The  re¬ 
sults  of  these  experiments  are  given  in  Pig.  1.13>  whence  it  is 
apparent  that  depending  on  the  magnitude  of  the  turbulence,  the 
critical  Reynolds  nximber  may  change  within  a  wide  range.  We  should 
also  note  the  fact  that  with  a  turbulence  of  0.08^6,  =  3*10^  is 

almost  reached;  this  is  six  times  greater  than  the  critical  Reynolds 
number  ordinarily  assumed  for  a  plate  in  a  subsonic  region. 

We  should  also  devote  some  attention  to  the  experiments  carried 
out  In  order  to  determine  the  critical  Reynolds  numbers  on  the  cone, 
said  experiments  carried  out  by  E.  Van  Drelst  In  a  wind  tunnel  [53].  The 
cone  being  tested  was  cooled  from  within  by  means  of  liquid  and  gaseous 
nitrogen;  In  this  manner  It  was  possible  to  determine  the  critical 
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Pig.  1.13.  The  influence  of 
flow  turbulence  on  the  criti¬ 
cal  Reynolds  number  of  a  flat 
plate  at  subsonic  velocities. 
1)  Turbulent  region;  2)  tran¬ 
sition  region;  3)  laminar  re¬ 
gion;  4)  turbulence  in  %  of 
flow  velocity. 


Reynolds  numbers  for  various  degrees  of  cooling.  The  test  results  are 
presented  In  Fig.  1. 14.  Curves  A,  B,  and  C  correspond  to  the  limits  of 
total  stability  (see  Fig.  1.12).  Curves  A',  B',  and  C,  passed  through 
the  experimental  points,  represent  the  limits  above  which  the  boundary 
layer  becomes  turbulent. 


ReffO'* 


Pig.  1.14.  Limits  of  stability  for  the  boundary. 

In  accordance  with  experiments  carried  out  on  a 
cone  be  cooled  In  a  wind  tunnel,  l)  =  2.9; 

Re^/meter  «  2.13*10^;  2)  =  2.7;  Re^/meter  » 

2.4-lo'^;  3)  Mg  =  3.65;  Reg/meter  ==  1.97*lo'^; 

4)  Insulated  walls;  A)  total  stability. 

We  should  note  the  fact  that  the  Reynolds  numbers  for  the  experi¬ 
mental  points  shown  In  Pig.  1.14  correspond  to  the  end  of  the  transi¬ 
tion  (buffer)  zone  between  the  laminar  boundary  layer  and  the  turbulent 
layer.  In  accordance  with  measurements  of  wall  temperatures  In  the  des¬ 
cribed  experiment,  the  difference  In  Reynolds  numbers  between  the  begin¬ 
ning  and  the  end  of  this  transition  amounts  approximately  to  2* 10^. 
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For  a  comparison  against  other  experiments,  curves  A',  B',  and  C* 
have  been  constructed  again  In  Pig.  1.15,  where  has  been  selected 
as  the  argument,  and  the  degree  of  wall  cooling  -  characterized  by  the 
ratio  —  has  been  chosen  as  the  paramenter.  For  an  adiabatic  wall 

Tst/^r  =  Ij  In  the  case  of  wall  cooling,  <(  1.  The  cooling  of  the 

surface  Increases  the  critical  Reynolds  numbers,  particularly  for  M|5  <(  3* 
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Fig.  1.15.  Critical  Reynolds  numbers  for 
the  transition  from  the  laminar  to  the 
turbulent  boundary  layer,  according  to 
experiments  carried  out  both  In  a  wind 
tunnel  and  In  flight,  l)  Beginning  of 
transition;  2)  end  of  transition; 

3)  cones;  4)  V-2  rocket;  5)  Viking  10 
rocket;  6)  plate;  7)  solids  of  revolu¬ 
tions;  8)  In  flight;  9)  in  wind  tunnel. 

In  Fig.  1. 15  the  experimental  points  obtained  during  flight  tests 
of  cones  and  rockets  have  been  plotted,  as  have  the  experimental  points 
obtained  in  tests  of  plates  and  solids  of  revolution  in  wind  tunnels  [44], 
[22].  The  values  of  the  critical  Reynolds  numbers,  obtained  in  these  ex¬ 
periments,  show  significant  scattering,  but  on  the  average  (for  the  end 
of  the  transition)  are  close  to  the  critical  Reynolds  numbers  obtained 
by  E.  Van  Drelst  (see  the  curve  for  T^^/T^  =  !)• 

The  factors  responsible  for  this  divergence  between  various 
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experiments  cannot  be  ascertained  with  a  sufficient  degree  of  accuracy; 
however.  It  should  be  pointed  out  that  the  magnitude  of  the  critical 
Reynolds  number  determined  during  the  experiment  can,  in  addition  to  the 
factors  Indicated  above,  be  affected  by  the  method  employed  to  deter¬ 
mine  the  beginning  and  the  end  of  the  transition  from  the  laminar  to 
the  turbulent  boundary  layer.  As  was  pointed  out  above,  the  turbulence 
of  the  stream  has  a  significant  effect  on  the  critical  Reynolds  numbers. 
In  the  experiments  used  for  the  preparation  of  Pig.  1. 15t  the  turbu¬ 
lence  of  the  air  varied,  and  Is  completely  unknown  for  a  number  of  cases. 

For  Mg  ^  2  In  the  case  of  -  1,  It  can  be  assumed,  on  the 

average,  that  the  beginning  of  the  transition  occurs  at  Reg  -  1.5*10°> 
and  that  the  end  of  the  transition  takes  place  at  Reg  =  3» 5* 10°«  The 
last  value  corresponds  to  the  Van  Drelst  curve  for  3- 0«  The  dif¬ 
ference  between  the  above- Indicated  Reynolds  numbers  corresponds  to  the 
difference  obtained  In  Van  Drelst 's  experments  [53]- 

Note  should  be  taken  of  the  fact  that  the  critical  Reynolds  num¬ 
bers  for  plates  and  solids  of  revolution,  according  to  the  cited  ex¬ 
periments,  show  no  pronounced  divergence  and  can  therefore.  In  approx¬ 
imate  terms,  be  assumed  to  be  Identical. 

Very  few  experimental  data  v/lth  respect  to  the  critical  Reynolds 
numbers  are  presented  In  the  literature  for  the  case  of  M  >  4.  The  fol¬ 
lowing  experiments  might  be  cited,  in  the  GALCIT  (USA)  [Guggenheim  Alr- 
onautlcal  Laboratory,  California  Institute  of  Technology]  wind  tunnel. 

It  was  found  that  In  the  case  of  M  =  5- 8  the  laminar  boundary  layer  on 
a  plate  Is  preserved  to  Reynolds  numbers  of  5. 5’ 10^.  In  this  case,  the 
laminar  boundary  layer  was  extremely  stable  [8].  In  flight  tests  (No.  27) 
of  the  V-2  rocket  It  was  found  that  at  Mg  =  4. 2  the  beginning  of  the 
transition  from  the  laminar  to  the  turbulent  boundary  layer  takes  place 
at  Reg  =  2.5*10^  [44]. 
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M.  Bertman  found  in  experiments  with  cylindrical  bodies  In  a 
wind  tunnel  at  M  =  6.9  that  the  critical  Reynolds  number  may  range  from 
1.2* 10^  to  4.5*10^»  depending  on  the  thickness  of  the  leading  edgej  In 
this  case,  with  an  Increase  In  the  thickness  of  the  leading  edge  the 
critical  Reynolds  number  Increases.  The  above-mentioned  critical 
Reynolds  number  of  4.5*10^  corresponds  to  the  Reynolds  number  calcu- 

4 

lated  on  the  basis  of  the  leading-edge  thickness  and  Is  equal  to  3* 10 
[38]. 

The  above —  enumerated  experiments  at  M  )>  4  demonstrate  that  with 
great  M(ach)  numbers,  the  critical  Reynolds  numbers  do  not  diminish 
and  lie  approximately  within  the  same  range  as  In  the  case  of  M  <  4; 
therefore,  for  M  >  4  the  same  critical  Reynolds  numbers  can  be  employed 
as  In  the  case  of  M  ^  4. 

In  the  case  of  surface  cooling,  l.e.,  for  critical 

Reynolds  numbers  Re^  will  Increase  as  a  result  of  the  reduction  In  vis¬ 
cosity  and  the  Increase  In  density  at  the  surface  of  the  wall.  In  approx 
Imate  terms.  It  may  be  assumed  that  the  critical  Reynolds  numbers  are 
Inversely  proportional  to  Tg^/T^.  This  Is  completely  confirmed  by  the 
Van  Drelst  experiments  for  the  case  of  M  )>  3*5>  said  experiment  shown 
In  Fig.  1.15.  For  M  <  3* 5  the  effect  of  cooling  will  be  more  pronounced, 
but  the  cooling  Itself,  caused  by  thermal  radiation,  will  not  be  great. 

In  the  case  of  great  M(ach)  numbers,  the  quantity  T^^/T^  can  be 
substantially  reduced  and  In  this  connection  we  can  anticipate  a  pro¬ 
nounced  Increase  In  the  critical  Reynolds  number.  However,  In  actual 
fact,  with  the  cooling  of  the  boundary  layer  the  displacement  thickness 
of  the  boundary  layer  diminishes  and  at  a  certain  cooling  ratio  the  re¬ 
duced  displacement  thickness  may  become  close  to  the  height  of  the 
roughness  protuberances.  In  this  case,  any  further  cooling  of  the  bound¬ 
ary  layer  will  not  only  fall  to  Increase  the  critical  Reynolds  number 
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but  will  actually  cause  It  to  drop,  since  the  roughness  Is  one  of  the 
significant  factors  resulting  In  the  Instability  of  the  boundary  layer, 
particularly  If  the  roughness  exceeds  the  displacement  thickness.  This 
phenomenon  may  be  referred  to  as  the  reversal  of  the  cooling  effect. 

The  presence  of  a  pressure  gradient  along  the  flow  has  a  pro¬ 
nounced  effect  on  the  critical  Reynolds  number:  In  the  case  of  a  pos¬ 
itive  pressure  gradient,  the  critical  Reynolds  numbers  diminish, 
whereas  In  the  case  of  a  negative  pressure  gradient,  they  Increase. 

With  a  negative  pressure  gradient  the  following  pressure  Is  exerted  on 
an  elementary  portion  of  air  having  a  volume  dxdydz 

——  dxdydz, 
dx 

which  accelerates  the  motion  of  this  portion  of  air.  The  decelerating 
force  of  friction,  acting  on  this  same  portion  of  air,  will  be: 

The  negative  pressure  gradient  coun¬ 
teracts  the  decelerating  effect  of  the  vis¬ 
cosity  forces,  as  a  result  of  which  the  sta¬ 
bility  of  the  boundary  layer  Is  Increased. 

In  the  case  of  supersonic  velocities, 
significant  negative  pressure  gradients  are 
formed  on  certain  types  of  wing  profiles  and 
nose-parts  of  airframes,  and  these  gradients 
enhance  the  stability  of  the  boundary  layer. 
Figure  1.16  shows,  as  an  example,  the  pres¬ 
sure  distribution  over  the  lenticular  pro¬ 
file  of  a  wing  at  M  =  2.13  [31]*  We  can  see 
from  the  cited  graph  that  with  small  angles  of  attack  the  negative 


Pig.  1. l6.  Distribution 
of  pressure  over  the 
lenticular  profile  of  a 
wing.  M  =  2.13;  "C  =  1.1; 

Re  =  6.4*10^j  l)  Upper 
surface;  2)  lower  sur¬ 
face. 
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gradient  covers  almost  the  entire  cord  of  the  profile. 

The  Influence  that  the  pressure  gradient  exerts  on  the  critical 
Reynolds  number  can  be  determined  in  approximate  terms  on  the  basis  of 
the  kinetic-energy  balance  of  the  boundary  layer.  To  simplify  the  der¬ 
ivation  of  the  corresponding  relationships,  we  will  assume  that  the 
velocity  of  the  flow  and  the  pressure  gradient  along  the  surface  are 
constant,  which  is  close  to  the  truth  for  thin  profiles  in  the  case  of 
small  angles  of  attack,  as  well  as  in  the  case  of  elongated  nose  sec¬ 
tions. 

We  choose  the  displacement  thickness  5  as  the  parameter  which 
characterizes  the  thickness  of  the  boundary  layer  in  the  case  of  a  pres¬ 
sure  gradient.  The  energy  Introduced  into  the  boundary  layer  per  unit 
time  as  a  result  of  the  pressure  difference  along  the  flow  will  be 


For  a  compressible  fluid,  the  displacement  thickness  is  equal  to 
(see  [48]) 


On  a  flat  plate,  for  the  laminar  boundary  layer,  in  an  incompress¬ 


ible  fluid  (see  [33]) 


«• 


(1.63) 


In  a  compressible  fluid,  the  displacement  thickness  can  be  found 
by  Introducing  into  Formula  (I.63)  the  ratio  of  the  shape  parameters 
and  the  densities: 


C.-1.72 


Hum  f  - 

hhV^' 


(1.64) 


Here  p 


* 


and  Re 


* 


are  calculated  for  a  definite  temperature. 


The  shape  parameter 


(1.65) 


where  6  Is  the  thickness  of  momentum  loss.  For  the  laminar  boundary 
layer  In  an  Incon^jresslble  fluid  (see  [33]»[^8]) 

//=2.59. 

and  In  a  compressible  fluid  (see  [48]) 


(1.66) 


Expression  (l.64)  was  derived  on  the  basis  of  Formulas  (1.25)  and 
(1.28)  and  the  following  relationship  (see  [48]) 

H  \  »*  /  Wck/ 

Substituting  the  value  of  6  Into  the  expression  for  Ep  In  ac 


*  *41/ 


cordance  with  Formula  (1.64)  and  assuming  the  quantity  (P  IJ-  ) 

to  be  constant  along  the  wall,  we  will  obtain: 

_  1 .72  ^  (P V)’'*  V'i'  j 

Tlie  energy  of  the  viscosity  forces  absorbed  by  the  boundary  layer 
per  unit  time  will  be 

dx. 


Taking  Into  consideration  the  value  of  c^.  g^j^  In  accordance  with 
Formulas  (1.25)  and  (1.28),  we  will  obtain: 

=.0,332  (pV)’'*  W'’| 

The  energy  lost  In  the  boundary  layer  will  be: 

E. 


According  to  Formulas  (l. 25), (l. 28),  and  (1.64),  we  will  find; 

tftm 


-  44  - 


Let  us  make  the  following  denotation: 


(1.67) 


so  that 

Let  (Hei^p)o  critical  Reynolds  number  for  the  given  flow  par¬ 

ameters  at  the  wall,  without  any  effect  due  to  the  pressure  gradient, 
and  let  'the  critical  Reynolds  number,  with  the  pressure  gra¬ 

dient  having  been  taken  Into  consideration.  The  value  of  (Re.  „)  will 

Kj?  p 

be  set  equal  to  that  Reynolds  number  at  which  the  energy  losses  In  the 
boundary  layer  will  be  equal  to  the  energy  losses  in  the  boundary  layer 
for  which  the  effect  of  the  pressure  gradient  has  not  been  taken  into 
consideration.  Denoting  the  quantities  which  cori^espond  to  the  flow  with 
a  pressure  gradient  with  the  subscript  "p,"  and  the  quantities  corres¬ 
ponding  to  a  flow  for  which  the  effect  of  the  pressure  gradient  has  not 
been  taken  into  consideration  with  the  subscript  "Q"  we  will  have  the 
following  for  the  critical  numbers : 

Let  us  substitute  the  value  of  E  from  Formula  (1.68)  into  the 
last  expression; 

i/p 

Since  the  quantity  c^.  is  inversely  proportional  to  x  '  , 

From  this  expression  we  can  find  the  sought  value  of  x  and  the 


ratio  Xp/xQ,  which  correspond  to  the  critical  Reynolds  numbers,  with 

(R^)/  ^ 

(R«icp)d  -"fo  ’  (1*70) 

In  order  to  sin?)llfy  the  determination  of  Xp,  a  graph  constructed 
in  accordance  with  Formula  (I.69)  is  presented  in  Fig.  I.17. 


-  45  - 


Pig.  1.17.  Graph  for  the  determination  of  the 
critical  length  x  in  the  presence  of  a  pres- 

sure  gradient  as  a  function  of  the  critical 
length  X-,  for  a  flow  without  a  pressure  grad¬ 
ient.  ^ 

For  an  experimental  verification  of  Formula  (I.69)  we  will  present 
it  in  a  different  form  for  the  M  <  1  case.  Let  us  introduce  the  Pohlhau- 
sen  parameter 


dp  V 
dx  liK,  ■ 

The  thickness  of  the  boundary  layer  at  whose  boundary  the  velocity 
differs  from  the  streamline  velocity  by  0.1^  (see  [33]>  page  118),  will 


be  equal  to 


S= 


64C 

/T5  ’ 


?tnd  consequently 

X=-18  — JC. 
dx 

Taking  into  consideration  that  for  the  case  M  <  1  the  value  of 
H  =  2. 59»  Formulas  (I.69)  and  (I.70)  can  be  presented  in  the  following 
form;  —  (t-nrma)*. 

(1.71) 

Figure  1.18  shows  a  congjarlson  of  Formula  (I.71)  and  the  experi¬ 
mental  data  derived  in  a  subsonic  wind  tunnel  [12]. 
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Fig.  1.18.  Comparison  of  theo¬ 
retical  and  experimental  values 
of  relative  critical  Reynolds 
number  in  the  presence  of  a  pres¬ 
sure  gradient. 

§1.6.  FEATURES  OF  HEAT  TRANSFER  IN  THE  CASE  OF  HYPERSONIC  VELOCITIES 

In  the  case  of  hypersonic  velocities  (M  >  5)  certain  mechanical, 
physical,  and  physicochemical  phenomena  develop  at  the  surface  of  the 
body  past  which  the  flow  is  moving,  and  these  phenomena  change  the 
physical  parameters  of  the  air  and  may  either  be  endothermic  or  exo¬ 
thermic  reactions. 

Ihe  dissociation,  recombination,  and  ionization  of  the  air  can  be 
Included  among  such  phenomena.  Moreover,  as  a  result  of  the  increase 
In  the  displacement  thickness  of  the  boundary  layer  in  the  case  of  hj^jer- 
sonic  velocities,  the  increasing  pressure  of  the  flow  at  the  surface 
of  the  body  begins,  in  this  case,  to  exert  a  pronounced  effect  on  heat 
transfer. 

The  above-mentioned  phenomena  are  of  great  significance  for  heat 
transfer  in  the  case  of  M  >  10. 

In  the  case  of  air  dissociation  produced  by  high  air  temperature, 
molecules  of  oxygen  Og  and  nitrogen  Ng  deponpose  into  atoms,  and  in 
this  case  a  substantial  quantity  of  heat  is  absorbed.  For  exangjle,  at 
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M  =  20,  the  temperature  of  an  ideal  gas  not  subject  to  dissociation 
must  be  17>500°C  in  back  of  a  normal  shock  wave,  and  for  the  case  of 
a  real  gas,  taking  into  consideration  the  absorption  of  energy  for  dis¬ 
sociation,  at  an  altitude  of  50  km  the  temperature  will  be  only  6500°  K. 
For  oxygen,  total  dissociation  sets  in  at  6,000°K,  and  for  nitrogen,  at 
10,000°K  [7]. 


Pig.  1. 19.  Relative  gas-kinetic  specific 
heat  flow  as  a  function  of  the  rate  of  air 
recombination,  l)  Plight  velocity,  2.4km/ 
sec;  2)  "Pyrex";  3)  oxides,  metals;  4) 
u  tcni/sec. 


Since  the  wall- surface  temperature  in  the  case  of  hypersonic  vel¬ 
ocities  will  always  be  lower  than  the  air  temperature  behind  the  shock 
wave,  an  Inverse  process  —  the  recombination  of  the  atomic  gases  —  may 
take  place  at  the  wall.  In  this  case,  thermal  energy  will  be  liberated. 
Below,  we  present  some  data  on  the  effect  that  dissociation  has  on  heat 
transfer  at  the  critical  point  of  a  blunt  body. 

The  quantity  of  heat  liberated  in  the  case  of  recombination  is  a 
function  of  the  rate  of  recombination  which,  in  tuni,  is  determined  by 
the  concentration  of  atoms  at  the  surface  and  the  catalytic  properties 

-  48  - 


of  the  wall.  Figure  1. 19  shows  the  theoretical  relationship  between  the 
relative  specific  heat  flow  to  the  wall  ('^J  =  q^/q^  _  and  the  rate  of 
recombination  at  the  wall  [8]j  the  radius  of  nose  curvature  In  this 
case  Is  assumed  to  be  0. 5ni  and  the  wall  temperature  Is  taken  as  700°K. 

It  follows  from  the  cited  graph  that  with  recombination  rates  In 
excess  of  lO'^  cm/sec,  the  effect  of  recombination  on  heat  transfer  be¬ 
comes  close  to  the  maximum  and  Is  a  weak  function  of  any  further  In¬ 
crease  In  the  rate.  In  the  case  of  an  Infinite  recombination  rate,  heat 
transfer  will  be  almost  the  same  as  In  the  case  In  which  there  Is  no 
dissociation  [28],  l.e.,  the  recombination  almost  completely  offsets 
the  effect  of  dissociation. 

Conversely,  at  rates  of  recombination  below  10  cm/sec.  Its  effect 
Is  Insignificant  and  heat  transfer  as  a  result  of  dissociation  may  dim¬ 
inish  substantially  as,  for  example,  with  a  flight  velocity  of  7.8  km/se 
this  drop  will  be  greater  than  by  a  factor  of  three. 

The  rates  of  natural  recombination  (without  catalysts)  are  quite 
low.  If  atoms  Impinge  on  the  wall,  the  rate  of  recombination  may  In¬ 
crease  manyfold.  In  this  case,  the  wall  acts  as  a  catalyst. 

The  catalytic  properties  of  the  wall  are  a  function  of  the  material 
of  which  the  wall  Is  made.  Metals  yield  higher  rates  of  recombination, 
whereas  for  nonmetalllc  walls  the  rates  of  recombination  are  substan¬ 
tially  lower.  Figure  1.20  [8]  shows  the  rates  of  recombination  for  oxy¬ 
gen  and  nitrogen  at  walls  made  of  various  materials;  these  rates  were 
derived  experimentally.  There  Is  little  data  on  rates  of  recombination 
and  therefore  the  graph  shows  curves  as  well  as  Individual  points 
(circles)  and  approximate  areas  (large  circles  with  chemical  denotations 
of  elements. 

A  comparison  of  the  curves  shown  In  Figs.  1.19  and  1.20  shows 
clearly  that  a  metallic  wall  will  receive  only  slightly  less  heat  than 
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at  which  recombination  takes  place  at  an  Inflnately  great  rate.  It 
should  be  borne  In  mind  that  with  an  Increase  In  temperature  the  rate 
of  the  catalytic  reactions  also  Increases;  this  can  also  be  seen  from 
Pig.  1.20.  Therefore,  In  first  approximation,  dissociation  need  not  be 
taken  Into  consideration  In  the  case  of  a  metallic  surface. 


u  cujceic 


Pig.  1. 20.  Catalytic  properties 
of  various  materials  In  recombin¬ 
ation  of  air  atoms  at  the  wall, 
ij  KCl;  2)  Li  Cl;  3)  "Pyrex"  glass; 

A)''Pyrex''  glass. 

In  the  case  of  a  wall  made  of  a  glass-based  material,  the  reduc¬ 
tion  In  the  transfer  of  heat  may  be  significant,  since  the  rates  of 
recombination  In  this  case  will  not  be  great. 

Chlorides,  oxides  of  metals,  and  certain  other  nonmetalllc  mater¬ 
ials  will  yield  Intermediate  heat-transfer  values.  It  should  be  pointed 
out  that  since  with  an  Increase  In  wall  temperature  the  rate  of  recom¬ 
bination  Increases,  a  comparatively  slight  drop  In  heat  transfer  with 
respect  to  the  transfer  of  heat  In  the  case  of  an  Infinite  recombina¬ 
tion  rate  may  be  permitted  for  the  group  of  materials  under  consider¬ 
ation. 

In  the  case  of  flow  past  plane  and  cylindrical  surfaces,  the 
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phenomena  of  dissociation  and  recombination  take  place  under  somewhat 
different  conditions  than  in  the  case  of  flow  past  the  area  close  to 
the  critical  point.  The  theoretical  investigations  of  heat  transfer  in 
the  case  of  a  laminar  boundary  layer  show  that  the  effect  of  dissoci¬ 
ation  on  heat  transfer  in  this  case  is  not  too  great  [19]- 

Given  high  temperatures,  characteristic  in  the  case  of  hypersonic- 
velocity  flow  past  bodies,  the  atoms  and  molecules  of  the  air,  in 
addition  to  dissociation,  will  also  be  subject  to  ionization.  Given  the 
air  temperatures  that  prevail  about  a  body  flying  through  the  atmosphere 
at  hypersonic  velocities,  the  degree  of  ionization  will  not  be  great. 

For  example,  with  M  =  20  approximately  1%  of  the  air  will  become  ionized 
(see  [3])-  This  degree  of  ionization  has  little  effect  on  heat  trans¬ 
fer. 

Although  for  M  <  25  the  effect  of  ionization  on  heat  transfer  is 
Insignificant,  it  maybe  Important  in  the  case  of  radio  communications 
[21]. 

In  determlng  heat  transfer,  it  is  necessary  to  know  the  tempera¬ 
ture  difference  T^  -  T^^  [see  Formula  (I.?)].  On  the  other  hand,  the 
specific  heat  of  the  air  is  Included  in  Formula  (1.29)  in  order  to  de¬ 
termine  the  heat-transfer  coefficient.  The  product  of  these  two  param¬ 
eters  must  characterize  the  change  in  gas  energy  as  the  gas  is  cooled 
from  a  temperature  of  T  to  the  temperature  T  . . 

I*  St 

In  the  case  of  hypersonic  velocities,  it  frequently  becomes  neces¬ 
sary  to  deal  with  an  extremely  great  temperature  difference  T^  —  T^^  in 
which  range  the  specific  heat  of  the  air  may  undergo  pronounced  change. 
Therefore,  Instead  of  the  following  product 

(n-rjc; 

it  is  more  effective  to  introduce  the  enthalpy 

i^jc,dr  (1.72) 
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and  to  express  the  specific  flow  of  heat  during  aerodynamic  heating  in 
the  following  form: 


9'»  *“  *1  4t)» 


(1.73) 


where 


a,=gs*c)(.*V%. 


(1.74) 


She  value  of  the  heat-transfer  coefficient  may  be  determined  by 

means  of  Formulas  (l.  30),  (l.3l)j  (1.32),  (1.39),  and  (I.62),  if  the 
* 

quantity  c  in  this  formula  is  removed. 

She  enthalpy  value  for  various  temperatures  can  be  found  from  Fig. 
1.21  (see  [48]). 

Formula  (1.62)  was  derived  for  the  determination  of  heat  transfer 
at  the  critical  point  of  a  rounded  nose.  In  the  case  of  great  hypersonic 
velocities,  a  structurally  simpler  formula  may  be  employed;  this  form¬ 
ula  was  derived  on  the  basis  of  experiments  carried  out  in  a  wind  tunne] 
(see  [4l]): 


(1.75) 


where  Is  the  first  cosmic  velocity  at  the  ground  equal  to  =  7. 

km/ sec. 


Fig.  1.21.  Enthalpy  of  air,  re¬ 
ferred  to  temperature,  as  a 
function  of  air  temperature. 


Since  in  the  derivation  of  Formuls 
(1.75)  we  have  taken  into  consideration 
the  experiments  that  were  carried  out 
at  hypersonic  velocities  all  the  way  tc 
cosmic  velocities,  the  effect  of  the  dj 
soclatlon  and  recombination  of  the  air 
has,  therefore,  been  accounted  for  in 
this  formula.  Ihe  experiments  which 
served  as  the  basis  for  this  fonnula 


were  carried  out  with  glass  "pyrex”  models,  and  this  material  exhibits 
weak  catalytic  properties  (see  Pig.  1.20).  However,  the  transfer  of  hea 
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In  these  experiments  was  measured  by  means  of  calorimeters  which  con¬ 
sist  of  platinum  plates  placed  into  the  stream  (see  [49]).  We  can  see 
from  Fig.  1.20  that  platinum  exhibits  high  catalytic  properties  (Ug^  > 
10^  cm/sec^,  and  therefore  Formula  (1.75)  yields  the  transfer  of  heat 
for  the  surfaces  that  exhibit  the  greatest  recombination  rates. 

The  enthalpy  value  1q  at  the  critical  point  may  also  be  derived  on 
the  basis  of  the  fact  that  all  of  the  kinetic  energy  in  the  decelerated 


flow  is  converted  into  thermal  energy.  Then 


JL  yi 

2  “ 


427g 


(1.76) 


where  1^^  is  the  enthalpy  of  the  air  in  an  undisturbed  flow;  this  value 
of  the  enthalpy,  as  well  as  Ig^  and  I^qq  found  from  Fig.  1.21. 

Because  of  the  existence  of  the  boundary  layer,  the  main  stream 
is  deflected  from  the  wall  by  approximately  the  displacement  thickness 
(see[48]).  With  great  hypersonic  velocities,  the  displacement  thickness 
of  the  boundary  layer  shows  a  pronounced  increase,  and  this  results  In 
an  Increase  for  the  angle  of  inclination  of  the  limit  of  the  dlsplacet- 
ment  thickness,  and  consequently  It  leads  to  an  Increase  of  pressure  In 
the  stream  and  at  the  wall.  The  Increased  pressure  results  In  Increased 
air  density  which.  In  turn,  produces  an  Increased  heat-transfer  coeffic¬ 
ient. 

The  Increase  In  pressure  produced  by  the  boundary  layer  can  be 

* 

determined  by  assuming  the  boundary  layer,  having  a  thickness  o  to 

have  solidified.  In  this  case,  a  flat  plate  will  be  similar  to  a  profile 
exhibiting  maximum  thickness  at  the  trailing  edge. 

The  local  angle  of  Inclination  for  the  limit  of  the  boundary  layer. 


corresponding  to  the  displacement  thickness,  will  be 

di* 

^  dx 


(  l.  77) 


This  angle  Is  added  to  the  local  geometric  angle  for  the  flow;  In 
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the  case  of  a  cone,  the  half-angle  of  the  cone  flare  increases  by  the 
angle  A<p. 

The  value  of  the  derivative  in  Formula  (1.77)  is  easily  determined 
for  a  plate  with  a  laminar  boundary  layer  by  Fonnula  (1.64); 

(1.78) 

dx  2 

Hence  we  can  see  that  with  an  Increase  In  the  displacement  thickness 
there  Is  an  increase  in  the  local  angles  of  stream  inclination. 

Having  determined  the  value  of  A(p,  the  formulas  presented  In  §1.  3 
can  be  used  to  find  the  heat- transfer  coefficient,  making  provision  for 
consideration  of  the  pressure  produced  by  the  boundary  layer.  As  a  resuH 
of  the  Increase  In  pressure,  the  air  parameters  for  the  boundary  layer 
change,  and  this  produces  a  change  In  the  displacement  thickness.  There¬ 
fore,  In  order  to  Increase  accuracy,  second-approximation  calculations 

5 

can  be  carried  out.  However,  in  the  case  of  M„  <  15  and  R  >  10^,  one 

00  0 

approximation  Is  sufficient. 

The  Influence  exerted  by  the  boundary  layer  on  pressure  is  a  strong 
function  of  the  Reynolds  number.  Figure  1. 22  shows  the  graph  for  pres¬ 
sure  at  the  surface  of  a  flat  plate  In  the  case  of  two  Reynolds  numbers 
with  respect  to  numbers.  We  can  see  from  the  graph  that  for  Re^  >  10^ 
the  effect  of  the  boundary  layer  on  pressure  Is  not  great. 

The  boundary  layer  exerts  a  significant  effect  on  the  pressure  In 
the  case  of  low  Reynolds  numbers.  Therefore,  It  Is  particularly  Impor¬ 
tant  to  take  Into  consideration  this  effect  In  the  case  of  flights  at 
great  altitudes,  as  well  as  In  the  determination  of  heat  transfer  In  the 
area  of  the  leading  wing  edge  or  close  to  the  pointed  nose  of  a  cone. 

The  boundary  layer  Is  generally  turbulent  for  Re,^  >  10°,  and  there¬ 
fore  Its  effect  on  pressure  Is  not  great. 

Since  the  effect  exerted  by  the  boundary  layer  diminishes  with  the 
Increase  In  the  Reynolds  number  (Fig.  1.22)  a  negative  pressure  gradient 
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can  be  observed  along  the  surface.  This  will  enhance  a  rise  In  the 
critical  Reynolds  number. 

The  heating  of  the  air  behind  the  shock  wave  Increases  the  Intens¬ 
ity  of  air  radiation.  At  near-cosmic  velocities,  the  Intensity  of  radia¬ 
tion  from  the  air  behind  the  normal  shock 
wave  Is  so  great  that  It  can  sometimes  pro¬ 
duce  a  radiant  energy  flux  that  impinges 
on  the  forward  surface  and  Is  of  the  same 
order  of  magnitude  as  the  heat  flows  which 
occur  as  a  result  of  aerodynamic  heating. 

The  importance  of  this  radiant-energy 
flux  In  the  over-all  transfer  of  heat  re¬ 
sulting  from  the  Increase  In  flight  veloc¬ 
ity  exhibits  an  extremely  rapid  rise.  In¬ 
deed..  If  the  specific  convection  heat  flow 
Qq  at  the  critical  point  of  a  blunt-nose 
body  Is  approximately  proportional  to 

the  heat  flow  of  radiant  energy  emitted  by  the  air  behlind  the  shock 
wave  Is  approximatly  proportional  to  (see  [l6]) 

We  can  also  see  from  the  last  relationship  that  the  radiant  heat¬ 
ing  being  examined  here  will  Increase  with  a  drop  In  altitude,  and  this 
will  take  place  more  Intensly  than  aerodynamic  heating.  It  Is  also  char¬ 
acteristic  that  with  an  Increase  in  the  radius  of  nose  curvature  the 
aerodynamic  heating  is  reduced  and  radiant  heating  Increases. 

The  gases  Included  In  the  con?)osltlon  of  the  air  do  not  radiate  In 
an  Identical  manner,  and  so  the  ln5)ortance  of  the  radiation  from  the 
various  components  of  the  air  changes  with  temperature  and  density.  For 
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Pig.  1. 22.  Relative  pres¬ 
sure  at  surface  of  plate 
at  a  =  0,  said  pressure 
arising  as  a  result  of 
the  bounda3?y  layer. 
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example,  at  a  temperature  of  8,000°  K,  l.e.,  characteristic  of  the  air 
behind  the  shock  wave.  In  the  case  of  near-cosmic  velocities  a  great 
role  Is  played  In  the  total  radiation  by  nitrogen  monoxide  NO  which  Is 
formed  at  high  temperatures  In  the  air,  and  molecular  nitrogen  Ng  Is 
also  extremely  Important.  At  a  temperature  of  12,000°  K,  characteristic 
of  velocities  close  to  the  second  cosmic  velocity,  molecular  nitrogen 
and  nitrogen  monoxide  dissociate  and  the  primary  emitters  are  the  atoms 
and  Ions  of  nitrogen  and  oxygen. 

The  thickness  of  the  emission  layer  which.  In  first  approximation, 
can  be  assumed  to  be  equal  to  the  distance  between  the  shock  wave  and 
the  critical  point,  exerts  great  Influence  on  the  Intensity  of  air  rad¬ 
iation.  The  distance  between  the  shock  wave  and  the  critical  point.  In 
turn.  Is  a  function  of  a  number  of  factors  Including  the  radius  of  nose 
curvature,  the  density  of  the  air,  and  the  degree  of  dissociation. 

In  connection  with  the  great  number  of  factors  which  affect  the 
radiation  of  the  air  behind  the  shock  wave,  a  purely  theoretical  deter¬ 
mination  of  the  emlsslvlty  of  the  air  behind  the  shock  wave  Is  extremely 
complex  and  none  too  reliable.  For  this  reason,  experimentation  plays 
a  more  Important  role.  At  the  present  time,  these  experiments  are  being 
carried  out  In  wind  tunnels. 

It  was  found  In  the  experiments  that  were  carried  out  that  given 
a  nose  radius  of  0. 3  m,  and  M  =  20  and  at  altitudes  of  the  order  of  30 
km,  radiative  heat  transfer  will  amount  approximately  to  10^  of  the  tota 
heat  flow  to  the  critical  point  (see  [13])*  A  similar  result  was  obtalne 
by  R.  Meyerott  [l6]  In  his  theoretical  determination  of  radiation  for 
a  nose  having  a  radius  of  curvature  of  0.  6  m  at  V  =  7- 5  km/ sec  for  an 
altitude  of  37  km. 

However,  at  greater  flight  velocities  the  heat  flow  due  to  the  rad¬ 
iation  of  the  air  may  exceed  the  convection  heat  flow.  For  example.  In 
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accordance  with  calculations  carried  out  by  K.  Gezll  and  D.  D.  Messon 
during  a  vertical  ascent  Into  the  atmosphere  at  a  velocity  of  10.7 
km/sec  of  a  sphere  having  a  radius  of  0.46  m,  at  altitudes  between  30 
and  40  km  the  heat  flow  due  to  radiation  was  greater  approximately  by 
a  factor  of  two  than  the  convection  heat  flow  (see  [l6]). 

A  rough  estimate  of  the  heat  produced  by  the  radiation  of  the  air 
behind  the  shock  wave  and  absorbed  at  the  critical  point  of  a  blunt- 
nosed  body  can  be  obtained  In  accordance  with  the  following  formula; 


(1.79) 


This  formula  was  derived  on  the  basis  of  an  evaluation  of  calcu¬ 
lation  values  from  the  approximate  functional  relationship  presented 
above.  This  formula  can  be  employed  for  the  determination  of  the  limits 
of  applicability  for  Formulas  (1.62),  (1.73)^  and  (1.75).  If  the  value 
Is  of  the  same  order  of  magnitude  as  Qq  according  to  the  cited  formula, 
where  If  the  value  Is  greater,  then  It  Is  evident  that  It  Is  Impossible 
to  determine  heating  by  means  of  the  aerodynamic  heating  formulas  alone. 
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Chapter  2 

RADIATIVE  AND  ABSORPTIVE  PROPERTIES  OP  BODIES 
§2. 1.  THE  BASIC  LAWS  OP  RADIATION  AND  THE  ABSORPTION  OP  RADIANT  ENERGY 

The  significance  of  radiation  In  the  general  thermal  balance  depends 
on  the  surface  temperature  and  the  magnitude  of  the  exterior  heat 
flows.  With  aerodynamic  heating  In  flight  at  velocities  of  M  <  3j 
the  significance  of  surface  radiation  Is  comparatively  slight;  how¬ 
ever,  It  should  be  considered.  In  the  case  of  a  flight  In  which  M  >  5# 
radiation  may  be  considerable;  here,  radiation  may  be  the  most  effec¬ 
tive  means  for  reducing  surface  temperature. 

Por  Instance,  In  the  example  of  aerodynamic  heating  discussed  In 
§3.1  where  M  =  5  at  an  altitude  of  30  km,  the  temperature  of  the  sur¬ 
face  at  a  distance  corresponding  to  10^  of  the  wing  chord.  In  the  ab¬ 
sence  of  radiation,  would  equal  860°C,  while  as  a  result  of  radiation 
It  Is  reduced  to  380°C,  1. e.,  by  480°C.  When  flight  velocity  increases, 
the  radiation  required  to  reduce  surface  temperature  Is  increased. 

The  heat  of  radiation  In  the  case  of  a  steady-state  thermal  regime 
has  the  greatest  Influence  on  the  surface  temperature,  whereas  In  the 
case  of  aerodynamic  heating 

In  the  case  of  a  nonsteady -state  regime  of  aerodynamic  heating,  the 
influence  of  radiation  decreases,  since, 

Por  space  vehicles,  radiation  Is  the  only  practical  method  of 
dissipating  the  heat  of  radiant  solar  emd  planetary  energy  absorbed  by 
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the  surface  of  the  craft.  Using  the  radiation,  we  may  also  dissipate 
the  heat  given  off  by  the  electronic  equipment. 

Also  of  great  significance  for  space  craft  is  the  reflective 
capacity  of  the  wall,  since  the  amount  of  the  absorbed  radiant  energy 
from  the  sun  and  planets  depends  on  the  magnitude  of  this  reflective 
capacity. 

The  radiation  of  an  absolutely  black  body  obeys  Planck's  lawj 


‘  (2.1) 
Where  is  the  Intensity  of  the  monochromatic  radiation;  X  is  the 
wavelength  of  the  radiation;  T  is  the  absolute  temperature  of  the  ra¬ 
diating  body. 

If  the  letter  E  denotes  the  total  energy  radiated  per  unit  body 
surface  per  unit  time,  then 


The  constants  c^  and  Cg  in  Jiq.  (2.1)  are  equal  to  (see  [^])j 
c^  =  8.85  •  10"^°  kcal  •  m^/sec;  Cg  =  1.44  •  10"^  m  •  deg. 

The  integration  of  Eq.  (2.2)  simultaneously  with  Eq.  (2.1)  en¬ 
ables  us  to  find  the  radiation  energy  of  an  absolutely  black  body; 
the  radiant  energy  is 


k<C|7'« 
■  I5cj 


.o7•^ 


where  o  is  the  Stef an -Boltzmann  constant,  equal  to 


(2.3) 


a  =  13  •  6  •  10  kcal/m^’ sec*  deg!^ 

Real  bodies  are  not  entirely  subject  to  Eq.  (2.3).  The  magnitude 
of  their  radiation  is  less  than  the  radiation  of  an  absolutely  black 
body,  so  that  for  a  real  body 


E  =  qi2  = 


eoT 


4 


where  e  is  the  emlsslvity  or  blackness  coefficient  and  e  <  1. 


-  59  - 


The  radiant  energy  falling  on  the  surface  of  the  body  is  com¬ 
pletely  absorbed  only  in  the  case  of  an  absolutely  black  body.  Real 
bodies  either  reflect  or  admit  a  part  of  the  incident  energy,  or 
they  go  through  both  processes  simultaneously.  The  relation  between 
the  absorption  and  radiation  capacities  Is  established  by  the  Klrchhoff 
law. 

If  Is  the  emissivity  for  a  given  spectral  wave,  and  Is  the 
absorption  coefficient  for  the  sajne  wave,  equal  to  a  portion  of  the 
absorbed  energy,  then  In  accordance  with  the  Klrchhoff  law,  at  a 
given  surface  temperature 

Px=*ei. 

The  coefficients  of  Intergral  radiation  and  absorption  at  a  given 
surface  temperature  do  not  equal  one  another  and  depend  on  the  spec¬ 
tral  composition  of  the  radiated  and  absorbed  energy.  The  value  of 
these  coefficients  may  be  defined  as  the  mean  integrals,  using  Eq. (2. 3); 


«» 

^  ^X2 


(2.4) 


(2.5) 


Here,  the  subscript  "1”  is  used  to  designate  the  magnitudes  of  the 
body  being  Investigated,  for  a  surface  temperature  Tj^,  while  the  sub¬ 
script  ”2"  is  used  to  denote  the  magnitudes  of  the  source  of  radiant 
energy  absorbed  by  the  body  at  temperature  Tg.  The  values  of  and 
E^2  may  be  determined  from  Ponnula  (2.1). 

The  intensity  of  the  monochromatic  radiation  has  a  maximum  for 
each  temperature.  As  an  example.  Pig.  2. 1  gives  the  curves  of  the  values 
of  for  6000°K  (temperature  of  sun)  and  1000°K;  these  curves 

were  constructed  according  to  Pormula  (2.1).  It  is  apparent  from  the 
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graph  that  when  the  temperature  of  the  surface  Increases,  the  maximum 
measurement  intensity  is  shifted  to  the  shorter  waves  and  becomes  more 
pronounced. 


Fig.  2. 1.  Relative  intensity 
of  the  radiation  of  an  abso¬ 
lutely  black  body  as  a  func¬ 
tion  of  the  length  of  the 
radiation  wave  and  the  temp¬ 
erature  of  the  body. 


Pig.  2.2.  Coefficients  of  mono¬ 
chromatic  radiation  as  a  function 
of  the  wave-length  for  nickel 
alloy  and  ceramics,  l)  ceramic 
A-417/235,  T  =  760  to  1000°C;  2) 
nickel  alloy,  T  =  760°  to  800°C. 


The  ratio  between  the  surface  temperature  of  an  absolutely  black 


body  and  the  wavelength  at  maximum  radiation  intensity  is  determined 
by  Wien's  law  which  may  be  obtained  from  Eq.  (2.1)  by  differentiating. 
According  to  this  law 


2898 
T  ’ 


(2.6) 


where  is  measured  in  microns.  Wien's  law  enables  us  to  estimate  the 
wavelength  in  the  region  in  the  region  in  which  the  radiation  intensity 
is  greatest. 

Solids  can  be  classified  on  the  basis  of  their  radiant  and  absorp¬ 
tion  properties  as  conductors  of  electric  current,  insulators,  and 
semiconductors.  Metals  are  associated  with  the  conductors,  while  ceram¬ 
ic  material  such  as  plastics,  lacquers,  paints,  etc.,  are  associated 
with  insulators. 

Conductors  and  Insulators  radiate  and  absorb  radiant  energy  in 
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various  ways.  Conductors  exhibit  low  emlsslvlty  at  normal  temperatures, 
while  Insulators  exhibit  high  emissivity.  Conductors  and  Insulators 
vary  in  the  spectral  composition  of  radiant  energy:  in  the  case  of 
conductors  radiative  capacity  is  most  important  in  the  shortwave  region, 
while  in  the  case  of  insulators  radiative  is  most  Important  in  the 
longwave  region.  As  an  illustration  of  this.  Pig.  2.2  gives  the  co¬ 
efficients  of  monochromatic  radiation  for  metal  and  ceramics  as  func¬ 
tion  of  wavelength;  these  coefficients  were  obtained  by  experiment  [40]. 
§2.2  EMISSIVITIES  FOR  VARIOUS  BODIES  AND  SURFACE  CONDITIONS 

The  temperature  of  the  surface  can  have  a  considerable  influence 
on  the  radiative  capacity  of  a  body.  Metals,  as  was  Indicated  in  §2.1, 
have  low  emissivity  in  the  case  of  long  waves.  However,  when  the  wave¬ 
length  decreases,  the  emissivity  of  metals  Increases.  If  we  take  into 
consideration  the  fact  that  according  to  Wien's  law  the  wavelength  of 
maximum- Intensity  radiation  decreases  with  an  Increase  in  the  surface 
temperature,  the  Integral  radiative  capacity  of  metals  should  Increase 
with  temperature.  This  is  confirmed  by  experiments. 

Figure  2.3  (see  [40])  gives  the  curves  of  the  integral  emissivity 
for  heat-resistant  alloys  (lower  curves).  As  can  be  seen,  at  tempera¬ 
tures  higher  than  600° C,  the  Increase  in  emissivity  can  be  considerable. 

The  integral-radiation  coefficients  for  some  metals  are  given 
in  Table  2. 1;  this  table  also  gives  the  temperature  ranges  of  the  ex¬ 
periment  and  the  corresponding  ranges  of  the  change  in  emissivity. 

It  is  also  apparent  from  Table  2. 1  that  as  the  temperature  in¬ 
creases,  emissivity  Increases. 

Emissivity  for  metals  depends  on  their  electric  conductivity.  The 
best  conductors  have  lower  emissivity.  We  can  see  from  Table  2. 1  that 
good  conductors  of  electricity,  such  as  aluminum,  gold,  copper,  and 
silver,  have  low  emissivity. 
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Pig.  2. 3*  The  influence  of 
the  surface  temperature  on 
the  coefficient  of  integral 
radiation  from  heat-resis¬ 
tant  metals,  l)  stainless 
steel;  2)  nickel-cobalt  al¬ 
loy;  5)  stainless  steel;  6) 
nlchrome;  A)  oxidized  sur¬ 
face  with  t  =  1150  C;  B)  un- 
oxldlzed  surface. 

It  is  possible  to  coat  a  metal  with  lacquers  or  paints  which  are 
poor  conductors  of  electricity  and  heat. 

In  this  case,  the  emlsslvity  Increases  sharply  and  approaches 
unity. 

Table  2. 2  gives  the  emlsslvity  for  certain  paint  coatings. 

It  should  be  noted  that  although  certain  paint  coatings  are  white, 
their  radiative  capacity  is  close  to  unity.  This  is  explained  by  the 
fact  that  their  good  reflective  capacity  (and  consequently,  poor  ra¬ 
diative  capacity)  is  manifested  only  in  comparatively  short  wavelengths 
corresponding  to  the  range  of  visible  radiation.  These  coatings  radiate 
infrared  rays  well.  This  same  phenomenon  may  also  be  found  in  the  case 
of  certain  other  dielectrics.  For  example,  at  20° C,  the  emlsslvity 
of  chalk  is  0.81  [9],  for  gypsvun,  it  is  0. 90  [17],  and  for  glazed 
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TABLE  2.1. 

Integral  Emlsslvlty  for  Some  Metals  (see  [17]) 


1 


HaiiMeHOBaHHe 


o 

^  noJiiiposaHHuA 

>Keac3o  no4iipoBaHiioe 
^  >Ke.ie30  JiHToe  Heo6pa6oTaHiioe 
CTa.ibiioe  JiHTbe  noaiipoeaHHoe 
^  Cra.ib  .lucTOBaa  iu.ii((|>OBaHHaii 
Y  HepiKaBeioiuaB  craab  |40) 

3  HitKeab  MiiCTbifl,  noaHpoBaHuuft 
9  XpoM 

HiixpoM  140) 

.11  HitKe.ib'KotSaabTOBbiii  cnaas  (40) 

12  Umhk  (9^,1h),  noaHpOBaiiHUft 

13  OmiiiKooaHHoe  aHcroBoe  Hceaeao 
3oaoTo  noaHpoBaiiHoe 

^5  jsUjib  noaiipoBaHHaa,  saeKipoaHTHaa 
l6Cepe6po  noaMpOBaHHoe,  MticToe 

17  riaaTHiia  noanpoBaHHan,  MHcraa 

18  BpoHaa  ujaM4)OBaiiHaa  (43) 

19  BpoHaa  nopHcraii  [43] 

20Moa«64eH  |I1] 

Moan6jiCH  Ill] 


225^575 

425-J-1C20 

925^1115 

770-J-1040 

940-1100 

480-:-800 

225.-;>375 

100-5-1000 

480-5-800 

450^800 

225-5-325 

28 

2254-635 

80-5-115 

2254-625 

2254-625 

65 

75-M75 

6004-1000 

15004-2200 


0,0394-0.057 
0,1444-0,377 
0,874-0,95  ; 

0,524-0.56  I 

0,554-0,61  I 

0,22-^0,575 
0,074-0,087 
0,08-5-0,26 
0,19-5-0.35 
0,25-s-0.65  I 

0.(MS^.053  I 

0.228 

0,018^0.035 

0.018^-0,023 

0.020^0,032 

0.054-5-0,104 

0,04 

0.57 

0,08^0,13 

0.19-i-0,26 


I)  Itemj  2)  polished  alumlnumj  3)  polished  Iron;  4)  crude 

cast  Iron;  5]  polished  steel  casting;  6)  sheet  ground  steel;  7)  stain¬ 
less  steel  [40];  8)  nickel,  pure:  9)  chrome;  10)  nl chrome  [40]; 

II)  nickel-cobalt  alloy  [4o];  12)  zinc  (99. 156),  polished;  13)  zinc- 
plated  sheet  iron;  14)  polished  gold;  15)  polished  copper,  electrolytic 
16)  polished  silver,  pure;  17)  polished  platinum,  pure;  lo)  ground 
bronze  [43];  19)  porous  bronze  [43];  20)  molybdenum  [11];  21)  molyb¬ 
denum  [11] 

porcelain  it  is  0.92  [17]. 

The  use  of  paint  coatings  is  limited  because  they  are  stable  only 
to  comparatively  low  temperatures,  e. g. ,  heat-resistant  enamel  IOI/19 
and  black  asphalt  lacquer  are  stable  when  subject  to  extensive  heating 
to  250°C.  Heat-resistant  enamel  with  an  aluminum  coating  is  stable 
approximately  to  400°C;  however,  aluminum  paints  have  a  low  emlsslvlty 
e  =0.35  [17]*  The  addition  of  alumlniim  powder  has  a  telling  effect. 

We  can  use  the  oxide  films  of  certain  metals  as  heat-resistant 
coatings  which  have  high  radiation  capacity.  For  example,  iron  oxide 
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TABLE  2.2. 


Emlsslvltles  for  Certain  Paint  Coatings 


Th!1  nOKpUTHfl  ^ 

r  c 

t 

/iMtepaxyp- 
huA  hctqii- 

HHK 

3.  ./laK 

40-r95 

0,80--0,95 

[171 

.  .laK  MepHuft  MaroBuft 

40-^95 

0,96-^0,98 

I«71 

J]aK  NepHbift  CjiecTBiuHA 

25 

0,87 

117] 

BaKeaitTOBuA  aaK 

80 

0,93 

I36J 

7  lUejiiflaK  MepiibiA  MaTOBuA 

74-M45 

0,91 

{171 

0  MacaflHbie  KpacKii  paajiHMHbix  ubctob 

1  100 

0.92-5-0,96 

[17] 

9  SiiiAeBafi  KpacKa 

20 

0,85-^0,95 

(36| 

LOCaika  aaMnoBag  (0,075  mm  h  Ooxbiue) 

!  40-f-370 

0.94 

[17] 

|.lCa)Ka  c  HCMJtKHM  CTeK;ioM 

lOO-J-185 

0,96-5-0.95 

[17] 

L2A;itoMHHKeBaji  KpacKa 

1  150-f-315 

'  0,35 

[17] 

l)  Type  of  coating  2)  literature  source;  3)  white  lacquer; 

4)  dull  black  lacquer;  5)  bright  black  lacquer;  6)  bakellte  lacquer; 
7)  dull  black  shellac;  8)  oil  paints  of  various  colors;  9)  enamel 
pain;  10)  lamp  black  (0.75  mm  and  more;  ll)  carbon  black  with  liquid 
glass;  12)  alumln\im  paint. 

In  the  temperature  range  from  500  to  1200° C  exhibits  an  emlsslvlty 
which  ranges  from  0.85  to  0.95;  the  emlslvlty  of  nickel  oxide  In  the 
temperature  range  from  65O  to  1255° 0  ranges  from  0.59  to  0.86  [17]. 

The  use  of  metal  oxide  films  to  Increase  radiation  capacity  Is 
advantageous  In  that  they  are  easy  to  apply,  while  the  cohesion 
strength  with  metal  In  the  case  of  a  thin  layer  Is  substantial.  The 
simplest  method  of  obtaining  a  film  Involves  the  preliminary  heating 
of  a  metal  and  keeping  It  at  a  high  temperature  In  air  for  several 
minutes. 

Figure  2.4  shows  the  emlsslvlty  curves  for  stainless  steel  at 
various  temperatures  as  functions  of  the  oxidation  temperature.  The 
tested  specimens  were  kept  at  the  oxidation  temperature  for  fifteen 

minutes.  The  graph  given  In  Pig.  2. 4  was  constructed  on  the  basis  of 

\ 

experimental  data  given  In  work  [40], 

It  Is  apparent  from  the  graph  that  the  oxidation  of  stainless 
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steel  at  50°C  for  fifteen  minutes  Increases  the  emlsslvlty  approximate¬ 
ly  to  e  =  0.9*  This  Is  also  true  for  a  number  of  other  alloys  (see 
Fig.  2.3). 


While  the  oxide  films  are  effective 
In  the  case  of  those  metals  Indicated  In 
Pig.  2. 3>  they  are  not  as  effective  with 
others.  For  certain  metals,  even  though 
the  oxide  film  does  Increase  the  radiation, 
emlsslvlty  remains  very  far  from  unity. 

For  example.  In  the  case  of  an  aluminum 
oxide  film,  obtained  by  means  of  soaking 
at  600” C,  the  emlsslvlty  Is  Increased  by 
a  factor  of  approximately  three,  but  the 
value  of  this  coefficient  Is  still  small; 
for  the  temperature  range  from  200  to 
600”  C  this  value  ranges  from  0.11  to  0.19  [17]. 

Table  2.3  gives  the  emlsslvltles  for  a  number  of  metals  with  oxi¬ 
dized  surfaces.  For  comparison,  the  emlsslvltles  of  these  same  metals 
with  nonoxldlzed  surfaces  are  also  given  In  this  table. 

The  emlsslvlty  Increases  as  the  surface  roughness  Increases.  If 
the  height  of  the  roughness  tubercles  Is  several  times  greater  than 
the  radiation  wavelength  the  emlsslvlty  of  the  rough  surface  as  a 
function  of  the  emlsslvlty  e  of  a  smooth  surface,  may  be  expressed  by 
the  following  formula; 

+2,8(1-.)*].  (2.7) 

The  structure  of  this  formula  Is  theoretically  derived  from 
Lambert's  law,  according  to  which  the  quantity  of  radiant  energy  at  an 
2uigle  f  to  the  normal  Is  proportional  to  cos  f. 


Pig.  2. 4  The  Influence  of 
the  oxidation  temperature 
on  the  emlsslvlty  of  stain¬ 
less  steel.  1)  Radiation 
temperature;  2)  oxidation 
temperature. 
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TABLE  2.3. 

The  Effect  of  Oxidation  of  a  Metal  Surface  on  its  Radiation  Capacity 


HaiiMetJOBaHiie 

/*C 

t 

TlHreparyp- 
HblA  HCTOH* 

1  HHK 

aHTan  no;iHpooaiiHaii 

770 -s- 1040 

0.524-0,56 

1171 

4CTa;ib  OKHcaeHHag  npii  600*  C 

200'r’600 

0.80 

1171 

^  HepwaBeiomag  craab 

480^800 

0,224-0.575 

[401 

g  H6p)KaBeiou(ag  craab,  OKHcaeiiiiaa  npH 
tJOO*  C 

450-5-840 

0,434-0,63 

[401 

7  Hep>KaBeK)ma5i  craab,  OKiicaeHHaa  npH 

480-r970 

0,624-0,73 

(401 

1000*  C 

8  Mep^aseiomag  craab,  OKHCaeiiHaa  npM 
1150*0 

480^1150 

0,854-0,96 

[401 

Q  Hiixeab  noanpoBaHHuff 

225-5-375 

0,074-0,087 

*  [171 

^QliiKe.ib,  OKHcaeiiHuri  npii  600®  C 

200-5-600 

0,374-0,48 

[171 

1  iHHXpOM 

4S0^m 

0.194-0.35 

[401 

■*'4lHxpoM,  OKiicacHHuii  npii  1I50®C 

480^1150 

0,754-0,90 

[401 

l^iiKeab-KoOa.ibTOBbiit  cnaae  j 

450-5-800  , 

0,254-0,65 

[401 

^%liiKeab-Ko6aabTOBbift  cn.ias,  OKMcaeBiibill 

4504-1150  1 

0,85-5-0,98 

[401 

npM  1150*  C 

1 

JIO.lIlpOBaflHUft 

2254-575 

0,0394-0.057; 

1171 

OKHcacHHud  npii  600®  C 

2C04;COO 

0,114-0.19 

(171 

XT^^CAb  noaMpOBaHHaa 

804-115  i 

0,0184-0,023 

(171 

L8>Me.ib,  OKifcaeHJiag  npii  600®  C 

2304-600 

0,574-0.87 

[>71 

l^naryBb  MarOBag 

504-350 

0,22 

1171 

^QlatyHb,  OKUCaetmaR  npii  600®  C 

2004-600 

0,6140,59 

(171 

'^Ihhk  noaiipOBaHHbiA 

2254-325 

0,04540,053 

(171 

OKIICaeHilblfi  HpH  400*  C 

400 

0,11 

(171 

1)  Item;  2)  literature  source:  3)  polished  cavSt  steel; 

4)  steel  oxidized  at  600°C;  5)  stainless  steel;  6)  stainless  steel 
oxidized  at  800°C;  7)  stainless  steel  oxidized  at  1000°C;  8)  stainless 
steel  oxidized  at  11502°C;  9)  polished  nickel;  10)  nickel  oxidized  at 
600°C;  11)  nlchrorae;  12)  nlchrome  oxidized  at  1150°C;  13)  nickel- 
cobalt  alloy;  14)  nickel-cobalt  alloy  oxidized  at  1150°C;  15)  polished 
aluminum;  I6)  aluminum  oxidized  at  600°C;  17)  polished  copper;  I8) 
copper  oxidized  at  600°C;  19)  dull  brass;  20)  brass  oxidized  at  600°C; 
21)  polished  zinc;  21)  zinc  oxidized  at  400°C. 

The  coefficient  2. 8  was  determined  on  the  basis  of  experiments. 
Figure  2.5  gives  the  results  of  an  experimental  check  of  Formula  (2.7). 

The  radiation  capacity  of  porous  bodies  is  substantially  increased, 
since  each  pore  may  be  looked  upon  as  an  absolutely  black  body.  For 
'  example,  porous  bronze  in  the  temperature  range  from  75  to  175°  has  an 
emisslvlty  of  e  =  0.  57j  while  ground  bronze  without  pores  at  65°  has 


an  emlssivlty  of  e  =  0. 04  [43]. 


At  the  present  time,  broad  distribution  has  been  given  to  lami¬ 
nated  plates,  which  together  with  good  structural  properties  possess 
superior  heat -Insulation  properties.  With  respect  to  radiation  capacity, 
they  behave  like  all  other  dielectrics,  l.e.,  they  have  high  emlssivlty 
at  temperatures  lower  than  those  which  characterize  their  heat  resis¬ 
tance. 

Laminated  plates  are  based  on  binding  substances,  generally  resins 
and  fillers.  In  laminated  plates  with  high  mechanical  properties,  we 
use  flberglas,  glass  fabric,  asbestos  fabric,  cotton,  or  silk.  Table 
2.4  gives  the  emlssivlty  of  certain  filler  materials.  We  can  determine 
the  radiation  of  resin  from  the  emlssivlty  of  bakellte  lacquer,  which 
is  given  in  Table  2.2.  As  can  be  seen  from  Table  2.4  (and  also  Table 
2.2),  the  emisslvities  of  the  materials  forming  the  laminated  plates 
are  close  to  0.9. 


t  HI  o,!t  f^s  o,t  {oe. 


Fig.  2. 5.  The  effect  of 
roughness  on  the  emissl- 
vlty.  l)  Experimental 
points;  2)  nlchrome; 

3)  nickel-cobalt  alloy; 

4)  stainless  steel;  5) 
brass;  6)  aluminum. 
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TABLE  2.4. 

Emlssivitles  of  Filler  Materials  for  Laminated  Plates ^  etc. 


HaHueHO0aHHe  ^ 

i 

a 

1 

TlMTeparyj^ 

HCTOMHHK 

2^3  Ac6ccTOBaa  <5yMara 

40-5-370 

0,93^0,95  ! 

H71 

Bywara  roHKaii,  HaKJieeHHaa  Ha  Meraa* 

JIH<ieCXHft^aKCT 

95 

0.89 

II7J 

^  jlepeBo 

70 

0,91 

[361 

6  Ctckjio 

20-4- 100 

0.94-5-0.91 

[Ill 

0  CTexao 

1  250-1-1000 

0,87-r0.72 

[111 

CreKao 

1  1100-5-1500 

0.70-5-0.67 

[111 

9  Vroab 

100-5-600 

0,81-5-0,79 

[III 

IcTHnc 

!  20 

'  0.8^0.9 

[111 

1)  Item;  2)  literature  source;  3)  asbestos  paper;  4)  thin 

paper  glued  to  a  metal  sheet;  5)  wool;  6)  glass;  7)  glass;  8)  glass; 

9)  carbon;  lO)  gypsum. 

§2.3.  ABSORPTION  COEFFICIENTS  FOR  VARIOUS  BODIES 

In  the  case  of  flight  at  very  high  altitudes  (higher  than  150  km) 
and  In  the  case  of  cosmic  flights,  a  basic  source  of  heat  is  the  sun. 
In  this  case,  as  will  be  shown  further  In  §3.3j  the  basic  parameter  in 
deteraiining  the  temperature  of  the  wall  Is  the  relation  of  the  absorp¬ 
tion  coefficient  to  the  radiation  coefficient  (emlssivlty) ,  that  is, 
P/e  [see  Formulas  {3>9),  (3*10)^  and  (3.11)].  Let  us  note  that  P  cor¬ 
responds  to  the  solar  radiation  spectmzm  while  e  corresponds  to  the 
spectrum  of  Infrared  rays  at  the  temperature  of  the  wall. 

Because  of  the  selective  nature  of  the  absorption  capacity,  the 
ratio  p/e  will  be  located  in  a  broad  range  of  values.  Figure  2.6  gives 
the  coefficients  P^  of  monochromatic  absorption  for  certain  polished 
metals  depending  on  the  length  of  absorbed  radiant  energy  (see  [15] 
p.  24l).  It  is  apparent  from  the  graph  that  the  absorption  capacity  of 
metals  for  the  energy  of  solar  rays  with  maximum  intensity  (the  wave¬ 
length  is  0.  5  M-)  increases  greatly  in  comparison  with  the  absorption 
capacity,  and  consequently  with  the  radiation  capacity  (the  Klrchhoff 
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law).  For  infrared  rays  corresponding  to  the  possible  temperature  of 
the  wall  (the  wavelength  of  a  ray  with  maximum  Intensity  Is  6  to  8  u)* 

In  the  case  of  metals  with  very  good  electric  conductivity  (silver, 
copper,  gold)  there  is  a  sharp  increase  In  the  absorption  coefficient, 
which  begins  at  wavelengths  of  approximately  1  u-  For  other  metals, 
the  absorption  coefficient  Increases  gradually  as  the  wavelength  de¬ 
creases. 

h 

0,8 

0,0 

O.h 

0,2 


0  /  2  J  «  S  S  TKhk 

Pig.  2. 6.  Coefficients  of  monochromatic 
absorption  of  various  metals,  ll  silver; 

2)  gold;  3)  copper;  4)  steel;  5/  alumi¬ 
num;  6)  nickel;  7)  steel;  8)  chrome;  9) 
nickel. 

The  Integral  absorption  and  radiation  capacities  are  determined 
to  a  significant  degree  by  coefficients  corresponding  to  absorption 
and  radiation  at  the  wavelength  of  maximum  Intensity. 

Table  2.5  gives  values  of  the  ratio  found  from  Pig.  2.6; 

In  addition,  determined  at  =  0.  5  |j,,  while  was  deter¬ 

mined  at  Xjjj  =  7  M-* 

It  is  characteristic  of  metals  that  the  ratio  ^^^2/^x1  greater 
than  one.  Silver  has  the  lowest  value  of  Pj^^^^Xl  gold  has  a 
higher  value  of  Pj^^^^Xl'  region  of  maxirnum  absorption 

intensity  (x^  =  0.5  |jl),  the  absorption  coefficient  for  gold  drops 
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sharply  and  consequently,  the  Integral  coefficient  of  absorption  will 
be  lower  than  the  value  indicated  in  Table  2.  5.  The  same  applies  to 


copper.  For  example,  for  gold  at  X  =  0.  7  ^l,  =  0. 06  and,  consequently. 


TABLE  2.3. 


The  Ratio  of  Absorption  and  Radiation  Coefficients 
Maxlmum-^Intenslty  WaYei§ngtii§, 


iMera^nji 


1 


2 

iipif  X— 0,5  MK 


o  Cepe6po 
?  XpoM 


0,06 

0,31 


2 

IipH  X=7  MK 


0.025 

0,05 


c;  AJlIOMIIHHft 
D  Craab 

I  Mcab 

HMKeab 

9  SoaoTo 


0,35 

0,48 

0,38 

0,40 

0,60 


0,05 

0,05 

0.C25 

0,025 

0,025 


(Emlsslvlty)  for 


2.4 

6,2 

7,0 

9,6 

15,2 

16,0 

24,0 


l)  Metal;  2)  3.^  at  X  =  0. 5  ;  3)  silver;  4)  chrome;  3)  alu¬ 
minum;  6)  steel;  7)  copper,*  8)  nickel;  9)  gold. 

It  should  be  borne  in  mind  that  a  finely  polished  surface  lowers 
the  value  of  the  ratio  3/e  In  comparison  with  its  value  for  ordinary 
polishing.  Actually,  In  the  case  of  fine  polishing,  the  average  height 
of  roughness  tubercles  will  be  0.1  to  0.5  while  In  the  case  of  or¬ 
dinary  polishing  it  will  be  0.  5  to  6  [l.  During  the  transition  from 
ordinary  polishing  to  fine- polishing,  the  radiation  coefficient  (emi- 
sslvity)  hardly  changes,  since  in  relation  to  the  wavelength  of  the 
radiation  (-7  p),  the  surface  in  the  case  of  both  ordinary  and  fine 
polishing  will  be  smooth.  As  concerns  the  coefficient  of  absorption, 
the  surface  with  ordinary  polishing  with  respect  to  the  wavelength  of 
maximum  absorption  intensity  (x^  =  0.5  p,)  will  be  rough,  since  the 
height  of  the  roughness  tubercle  will  be  several  times  greater  than 
the  wavelength  (1  to  12  times  greater).  Therefore,  the  absorption 
coefficient  in  the  case  of  fine  polishing  should  be  lower  than  in  the 
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case  of  ordinary  polishing,  and  consequently,  the  value  of  the  ratio 
P/e  should  be  lower. 

To  shield  against  solar-ray  heating  we  may  use  a  dielectric  sur¬ 
face  coating  which  has  high  emlsslvlty  and  a  small  solar -energy  ab¬ 
sorption  coefficient. 

In  particular,  it  is  possible  to  use  a  paint  coating.  Table  2.2 
shows  that  the  emlsslvlty  e  of  paint  coatings  is  on  the  average  equal 
to  0.9,  while  the  absorption  coefficient  of  white  paint  for  example 
ranges  from  0.12  to  0.26.  Consequently,  the  ratio  p/e  will  range  from 
0. 13  to  0.  3. 

Certain  materials  have  very  low  absorption  coefficients  when 
emlsslvlty  is  of  the  order  of  0.  9*  For  example,  gypsum  has  an  absorp¬ 
tion  coefficient  of  0.05  to  0.10,  while  magnesium  oxide  even  has  an 
absorption  coefficient  of  0. 1  to  0.  2. 

Table  2.6  gives  the  absorption  coefficient  of  certain  materials 


at  room  temperature. 

TABLE  2.6. 

Coefficients  of  the  Integral  Absorption  of  Solar  Energy  of  Certain 
Materials 


JlHreparypHuA  g 
HCrOlHHK 


0,26 

0,38 

0,26 

0,45 

0,74 

0,66 

0,12-5-0,26 

0,97-5-0,99 

0,35-^0.40 

0,25-5-0,35 

0,05-5-0,10 

0,01•^0,02 

0,27 

0, 15-5-0,20 


I17J 

117) 

1I7J 

I17J 

1171 

I17J 

[17J 

121 

(2) 

12] 

12] 

l>71 

(21 


1  MarepHaa 


3  AaioMiiiiHff  noanpoBaHHiiifl 
^  AaioMHHiiA  MaroBbiA 
^  A'eab  no.iiipoBaHHaa 
^  /Keaeao  noaHposaHHoe 
OKiicacHHoe 

^IKeaeao  oUHHKOBanHoe 
^^h'pacKa  Ocjiaa 
npacKa  'lepuafl 
A.iioMHHHCBaii  KpacKa 
l^ric.iaa  4*9P4>0P<>'>ait  »Maa;, 

itOxHCb  HarHHB  (MaTOBaa) 
^ByMara  6eaa* 


l6 


Eynara  BiTMaHCKaii  5eaaa 
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Key  to  Table  2.6.  l)  Materials^  2)  literature  source;  3)  polished  alu¬ 
minum;  4)  dull  aluminum;  5)  polished  cooper;  6)  polished  iron;  7)  oxi¬ 
dized  iron;  8)  iron;  9)  white  paint;  10;  black  paint;  11)  aluminum 
paint:  12)  white  porcelain  enamel;  13)  gypsum;  l4)  magnesium  oxide 
(dull);  15)  white  paper;  I6)  white  Whatman  paper. 


Manu¬ 

script 

Page 

No. 

58 

58 

58 

66 


[List  of  Transliterated  Symbols] 

“  '^iz  ~  ^Izluchayemyy  “  “^radiated 

~  '^a  “  ^aerodinamlcheskyy  ~  ^aerodynamic 

^nar  “  *^nag  ~  ^nagrev  ~  '^heating 

6  =  £  =  £  =  £ 
lu  sh  sherekhovatyy  rough 
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Chapter  3 

THE  TEMPERATURE  OP  THE  WALL  DURING  A  STEADY -STATE  THERMAL  PROCESS 
§3.1.  THE  TEMPERATURE  OP  THE  WALL  WITHOUT  INTERNAL  COOLING 

The  steady-state  thermal  process  is  characterized  by  a  tempera¬ 
ture  which  Is  constant  with  respect  to  time  at  any  point  of  the  flow 
and  of  the  streamlined  body.  This  process  is  of  interest  in  the  case 
of  a  long  flight  with  constant  velocity  at  constant  altitude  or  in  the 
case  of  space  flights  at  sufficiently  great  distances  from  the  earth 
or  from  another  planet. 

In  the  case  of  the  steady -state  thermal  process,  the  equipment 
is  not  heated  and  consequently  q^^^  =  0.  When  there  is  no  cooling  of 
the  inner  surface,  the  thermal-balance  equation  (1.13),  with  consid¬ 
eration  of  Formula  (l. 11),  will  be: 

+  ^M3-^T=0-  (3.1) 

Here  the  temperature  of  the  wall  will  be  constant  with  respect 
to  its  thickness.  This  is  the  temperature  at  which  the  heat  flows 
approaching  and  leaving  the  plating  reach  a  state  of  equilibrium. 
Consequently,  this  temperature  is  called  the  equilibrium  temperature. 

The  heat  loss  q^  along  the  surface  is  produced  by  a  temperature 
gradient  along  the  surface;  for  a  cylindrical  surface  with  the  length 
J.  from  the  generatrix  perpendicular  to  the  flow,  the  quantity  of  heat 
given  off  per  second  will  be: 


where  x  Is  the  distance  from  the  Investigated  point. 

The  quantity  of  heat  entering  the  same  surface  element  will  be: 


-X 


Here  "1"  indicates  the  cross  section  through  which  the  heat  enters, 
while  "2"  indicates  the  cross  section  through  which  the  heat  exits.  If 
the  distance  between  cross  sections  is  dx,  then 


It  is  apparent  from  the  formula  obtained  that  the  heat  loss 
through  the  cross  section  of  the  plating  will  take  place  when  the 
temperature  gradient  along  the  chord  is  not  constant.  Substantial 
losses  may  be  observed  in  the  vicinity  of  the  pointed  leading  edge,  at 
the  nose  of  the  cone,  at  the  point  at  which  the  transition  from  a 
laminar  to  a  turbulent  boundary  layer  takes  place,  and  at  the  places 
where  the  streamline  profile  is  bent. 

Since  the  heat-transfer  coefficient  and  the  recovery  factor  in 

q^,  and  also  the  emisslvity  in  depend  on  the  temperature  of  the 

wall,  we  can  most  conveniently  determine  the  latter  graphically.  For 

this  purpose,  two  or  three  values  of  are  given  and  we  find 

from  Formula  (3.1);  then  we  construct  <l^g^g  from  T^^,  and  taking  into 

consideration  the  fact  that  q^„„  =  0,  we  find  T„4.. 

^nag  s  t 

The  number  of  points  for  the  determination  of  is  taken  de¬ 

pending  on  the  available  possibilities  of  obtaining  the  necessary 
accuracy  for  a  preliminary  evaluation  of  Tg^.  For  example,  if  T^  < 

<  500° K  and  there  is  no  cooling,  the  temperature  of  the  wall  in  the 
majority  of  the  cases  will  be  only  somewhat  less  than  the  recovery 
temperature;  consequently,  we  may  give  two  values  of  T  .  differing 
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from  each  other  by  30  to  50°.  If  the  preliminary  estimation  of  can 
be  made  only  with  an  accuracy  of  several  hundred  degrees,  we  must  cal¬ 
culate  three,  and  possibly  even  four  values  of 

It  should  be  noted  that  the  wall  radiates  from  two  sides;  how¬ 
ever,  the  radiation  of  heat  from  whlthln  the  structure  in  the  case  of 
a  steady-state  thermal  process  will  be  compensated  by  the  absorption 
of  the  heat  radiated  by  the  structure  or  the  opposite  wall.  Consequent¬ 
ly,  Formula  (3. 1)  should  take  into  consideration  only  one-sided  rad¬ 
iation. 

In  order  to  make  more  clear  the  above  method  for  determining  the 
temperature  of  a  streaimllned  surface  in  the  case  of  a  steady -state 
thermal  process,  we  give  below  an  example  of  the  determination  of  the 
temperature  of  the  lower  wing  surface. 

The  data  derived  are  as  follows;  M  =  5,  H  =  30  km,  and  a  =  10°, 

00 

the  profile  is  lenticular, "c  =  0.03,  chord  length  is  8  m,  and  surface 
emlssivity  is  e  =  0. 8.  This  calculation  is  for  the  lower  surface. 

We  will  determine  the  coefficients  c^^,  Cg,  and  c^  from  Pig.  1.5; 
for  =  5)  they  will  be  c^  =  0. 4l,  c^  =  1.22,  c^  =  1.92.  The  angle 
of  surface  inclination  for  the  profile  at  the  leading  edge  in  relation 
to  the  chord  is  =  2c  =  0. 06.  The  angle  between  the  lower  surface 
of  the  leading  edge  and  the  direction  of  the  flow 

?  =•  “  +  ^0=  0, 1 75  +  0,06  =  0,235. 

The  pressure  behind  the  compression  wave  is  determined  from  Formula 
(1.40); 

f “  0.  *  885. 

Hence,  from  Formula  (1.59)  we  find; 

PjP^  ”  1  +  = 4,3. 
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Let  us  determine  the  relative  density  behind  the  compression 
wave  from  Formula  (1.57): 

- 2.6. 

—4-6 

Table  3. 1  gives  the  calculation  formulas  and  the  numerical  cal- 
,  culatlons  for  the  flow  characteristics  along  the  surface  of  the  profile. 

To  determine  the  temperature  of  the  surface,  we  calculate 
for  two  given  temperatures  and  then  graphically  or  analytically  we  find 
T  ,  (from  the  linear  function  between  q  and  T., )  at  which  q„-„  =  0. 

S  u  Hel^  O  t 

For  a  laminar  boundary  layer,  we  assume  surface  temperatures  of  65O 
\  and  700°K  while  for  a  turbulent  boundary  layer  we  assume  temperatures 

of  700  and  850°  K. 

To  determine  the  temperature  of  the  adiabatic  wall,  let  us  find 
the  Prandtl  number.  Prom  r^  =  0.  83  we  determine  the  temperature  of  the 
adiabatic  wall  according  to  Formula  (l. l)  for  a  laminar  boundary  layer 
and  for  the  average  Mg  =  3.66  and  Tg  =  3*53: 

r,  ==  7'^  ( 1  +  0,2r,M?)  =  1 1 30*  K: 

then  from  Formula  (1.6)  where  T„.  =  700°K,  we  find  the  determining 

0  U 

temperature; 

r* = ft  +  0,5  (T„  -  Ti)  4-  0.22  (r,-  T%)  »  697*  K. 

Prom  Pig.  1.1,  let  us  find  the  Pr  value  for  the  obtained  determin¬ 
ing  temperature  (according  to  Van  Drelst); 

Pr* =0,684. 

Similarly,  for  a  turbulent  boundary  layer,  from  r^  =  0. 89,  for 
average  Mg  =  3*87  and  Tg  =  325°K  we  will  have; 

7',=  n90*K,  7^  =  703*K  h  Pr*«0,684. 
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At  other  wall  temperatures,  the  Prandtl  number  changes  very- 
little.  With  =  650° K  in  the  case  of  a  laminar  boundary  layer, 

Pr*  =  0.683;  when  Tg^  =  850°K  in  the  case  of  a  turbulent  boundary  layer, 
Pr*  =  0.688.  This  change  in  the  Prandtl  number,  in  accordance  with 
Formulas  (1.4)  and  (l.  5),  has  very  little  Influence  on  r^^  and  r^.  In 
this  connection,  the  found  values  of  are  used  as  constants  for  fur¬ 
ther  calculations. 

It  should  be  borne  in  mind  that  since  the  braking  temperature  does 
not  change  when  the  local  flow  characteristics  change,  the  Prandtl 
number  may  be  taken  as  constant  for  all  points  of  the  surface;  con¬ 
sequently,  we  may  assume  that  the  temperatures  of  the  adiabatic  wall 
calculated  above  are  the  same  for  all  points  on  the  surface.  The  values 
of  c  *  =  0.257  are  also  constant  for  all  points  on  the  surface  (see 

*r 

Pig.  1.2). 

The  specific  radiation  heat  flow  will  depend  only  on  the  surface 
temperature  and  consequently  according  to  Formula  (l. l4)  it  will  equal: 

at  a  wall  temperature  of  Tg^  =  650‘’K . q^^  =  1.94  kcal/m  sec 

at  a  wall  temperature  of  Tg^  =  700°K . q^^^  =  2.  6l  " 

at  a  wall  temperature  of  Tg^  =  850°K . =  5*68  " 

The  heat  flow  of  solar  radiation  will  have  no  effect  on  the  lower 
surface;  the  heat  flows  from  the  earth  will  be  insignificant  (fraction¬ 
al)  in  comparison  with  the  heat  flows  from  the  aerodynamic  heating 
and  radiation  and  they  may  therefore  be  neglected. 

The  beginning  and  the  end  of  the  transition  of  the  laminar  boun¬ 
dary  layer  into  a  turbulent  boundary  layer  are  roughly  determined  at 
first  on  the  basis  of  a  visual  estimate  of  T„,./T^  »  0.7.  Then,  accord- 
Ing  to  the  data  given  in  §1.5,  the  critical  Reynolds  number  for  the 


I.SIC* 


beginning  of  the  transition 

Re» 


2,14.  lO*. 


TABLE  3.1 


The  Determination  of  Plow  Characteristics  on  the  Profile 


Be4H^Hlla  H  4’opMyaa^ 

M  4)op- 
Myau 

MKcaeKHue  sHaseHHa  ^ 

BejiHHmoti'x  3aAaeMcii4 

0 

:  0,05 

O.IO 

0,20 

0,60 

'  0,75 

1,00 

1 

o 

1 

ac 

— 

0,06 

0,054 

0,048 

0,036 

0 

—0,03 

— 0,C6 

V  *  +  « 

0’235 

1  0,229 

0,223 

0,211 

0,175 

0,145 

0,115 

(1.40) 

0,1885 

0,1807 

0,1734 

0,1586 

0,1195 

0.0900 

0,0661 

=  i  +  o.7MiJ, 

(1.59) 

4.3 

4,16 

4.04 

3,78 

3,11 

2.58 

2,16 

Pj/Pc  =  Pa//’.) 

(1.58) 

1 

0,980 

0,957 

0,911 

0,795 

0.694 

0,611 

P«fP.“(P»/Pc)(Pc/P-) 

2,6 

2,55 

2,49 

2,37 

2,07 

1,81 

1,59 

-  (a//’-)(P«/P«) 

(1.48) 

1,65 

1,63 

1.62  i 

1,595 

1,50 

1,43  1 

1,36 

mJ-5((7-./7-j)(1+0.2M*J-.1] 

(1.49) 

13,2 

13,4 

13.5  ; 

13.7 

15,0 

16,0  : 

1,90 

M, 

— 

:  3,64  i 

3,66 

3,67  1 

3,70  1 

3,87 

4, CO 

4.12 

10^XPl“fP-(Pj/p.)]Xl03 

— 

I  4,81 

4,73 

4.61  1 

4.39 

3,84 

3.35  , 

•2,94 

— 

357 

353 

351 

346 

325 

310 

294 

V,- 

— 

1375  1 

1375 

1375  , 

1375 

1395 

1410  j 

1410 

lOe  X  A  -Il,4910-7rJ'V(7',  +  110)]  X  106 

(1.34) 

—  I 

2.14 

2.12  1 

2,10 

2,01 

1,94  1 

1,83 

10«  X  V,  -  (A/Pj)  X  103 

—  s 

4,53 

4,60 

4,78 

5,24 

5,79 

6,24 

10-6xRe, -(V,4-;v,)x10-6 

(1.20) 

1.21 

1 

2.4 

4,61  i 

10,65 

14,6  1 

18.1 

l)  The  value,  and  formulaj_2)  number  of  formula;  3)  numerical 
values;  4)  the  magnitude  x  is  given. 

the  critical  Reynolds  number  for  the  end  of  the  transition 


Rej 


3.5-10* 

TcrlTr 


5-10*. 


It  is  apparent  from  Table  3* 1  that  the  beginning  of  the  transi¬ 
tion  will  take  place  at  a  distance  approximately  corresponding  to  10^ 
of  the  chord,  while  the  end  of  the  transition  will  take  place  at  a 
distance  corresponding  to  20$^  of  the  chord. 

The  calculated  values  of  are  given  Table  3*2.  Prom  the  known 

values  of  =  0.  If  the  temperatures  of  the  wall,  which  are  as¬ 
sumed  in  the  determination  of  are  designated  as  and  Tg,  and 

the  values  of  q^^  corresponding  to  these  are  designated  as  q^^  and  qg, 

then  from  the  linear  function  between  q„«^  and  T„4.  (for  the  Interval 

^nag  St  ' 
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under  consideration)  we  will  have: 


T  —T  A-a 

The  values  of  determined  from  this  formula  are  also  given 
in  Table  3.  2. 

In  the  determination  of  we  take  q.  =  0,  since  for  most  of 

nag 

the  chord  length  the  temperature  gradient  is  almost  constant.  At  the 
point  of  transition  from  a  laminar  boundary  layer  to  a  turbulent 
boundary  layer  (at  a  distance  corresponding  to  10  to  20^  of  the  chord), 
the  sign  of  the  temperature  gradient  will  change,  as  a  result  of  which 
heat  flows  will  develop  along  the  plating;  consequently,  the  breaks 
in  the  temperature  curve  should  be  rounded  off. 

On  the  basis  of  the  temperatures  obtained  in  Table  3*2,  the  crit¬ 
ical  Reynolds  numbers  can  be  made  more  exact. 

For  the  beginning  of  the  transition 

0,575; 

T,  1130 

for  the  end  of  the  transition 

*--^=0,705. 

Consequently,  the  critical  Reynolds  number  for  the  beginning  of  the 
transition 


for  the  end  of  the  transition 


According  to  Table  3.1,  the  points  which  are  at  distances  corres¬ 
ponding  to  11  and  21$^  of  the  chord  correspond  to  these  Reynolds  numbers. 
The  temperature  values  found  along  the  chord  on  the  lower  surface 
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of  the  wing  profile  are  given  In  Pig.  3. 1. 

TABLE  3.2. 

The  Determination  of  the  Surface  Temperature  of  a  Profile  In  a  Steady- 
state  Thermal  Process 


h  4>opMy4a 

•  '  1 

AT 

(top- 

My;iu 

1  2 

Witc.ieiiHue 

aiiawetiKR 
_ 2 

;taMiiHapHbiii  norpa* 

lIHMHUSl  c.ioh  ^ 

- ^ - „ 

TypOy^ieiiiHUli  norpainmabiM 

caofl 

4 

BeaimHuoji  x  aaAaeMcn  6 

!  ! 

0,05 

0,10  1 

1 

0.20 

1  1 

0. 

1 

£0 

0.75 

uo 

Be;iH>(i(Hofi  Tfi  aasaeMca  7 

650 

700 

650 

7C0 

700 

850 

700 

850 

700 

850 

700 

850 

7-*  -  r,  +  0.5  (T„  -  n)  +  0.22  {Tr  -  r,) 

(1.6) 
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670 
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1 
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i  778 

1 
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103  X  p*  =  (p,7-j/7-.)  X  103 

(1.33) 

:  2,49 

2.40 

2,42 

2,34 

2,19 

1.94 

'  1,78 

1,605 

1,49 

1,345 

1,25 

1,123 
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(1.34) 

3,32 

3.39 

3,32 

3.39 

3,38 
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3,61 
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1 
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1 
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3.70 
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— 
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— 

— 

— 
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1 

7,9 
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(1.7) 

2,74 

2,43 

1,91  ' 

1,70 

8,24 

5,26 

5,89 

3,81  ' 

4,76i 

3,09 

3,87 

2,52 
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(3.1) 

0,80 

-0,18 

0,03 

-0,91 

5,62 ; 
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3,28 
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-2.59 

1,26 

-3,11 

^CT  (7*  —  ^2) 

— 

691 

648 

840 

1  796 

770 
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l)  Value  and  formula;  2)  number  of  formula;  3)  numerical  values;  4)  tur¬ 
bulent  boundary  layer;  5)  lamlneu?  boundary  layer;  6)  the  value  of  x  Is 
assumed;  7)  the  value  of  T  ^  Is  assumed.  "" 


Pig.  3*1*  Temperature  distribution 
along  the  chord  on  the  lower  sur¬ 
face  of  the  wing  profile.  M  ■  5; 

H  -  30  km;  o  -  10®;  e  -  0.8!* 


§3.2.  THE  TEMPERATURE  OP  THE  WALL  WHEN  ITS  INNER  SURFACE  IS  COOLED 

With  a  cooled  wall  and  the  coolant  absorbing  the  specific  the 

heat-balance  equation  for  small  surface  heat  losses  will  be 

+  (3.2) 

With  a  cooled  outer  surface^  the  wall  temperature  will  be  constant  In 
width,  and  Its  determination  will  follow  §3. 1* 

With  an  liuier  cooled  wall  surface  (e.g. ,  by  circulating  a  liquid 
coolant)  the  wall  temperature  will  vary  In  width.  A  temperature  distribu¬ 
tion  diagram  for  the  wall  width  Is  given  In  Pig.  3*2.  For  generality. 

Pig.  3.2  assumes  that  the  wall  consists  of  Insulation  and  plating.  Great 
aerodynamic  heat  flow  cools  the  inner  surface  and  Increased  surface  ra¬ 
diation  properties  do  not  yield  acceptable  wall  temperatures.  Here,  the 
radiant -energy  heat  flow  q^  can  be  neglected,  and  Eq.  (3.2)  Is  simplified. 

Prom  this  we  can  find  the  outer  surface  temperature  '^iz*  To  determine 
the  plating  temperature  we  study  the  Insulation's  thermal  conductivity. 

With  a  steady -state  thermal  process,  neither  the  Insulation  nor  the 
plating  are  heated;  consequently,  the  specific  heat  flow  will  pass  through 
any  surface  parallel  to  the  outer  Insulation  surface 

specific  heat  flow  will  pass  through  any  similar  plating  surface 

If  we  assume  that  within  the  temperature  difference  across  the  width, 
thermal  conductivity  does  not  depend  on  temperature,  we  will  have 

Consequently,  the  insulation  and  plating  temperature  will  change 
linearly  In  width.  The  temperature  differences  across  the  Insulation 
and  plating  will  be 


111 


A7b 
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(3.4) 

(3.5) 


Pig.  3. 2.  Diagram  of  the  temper¬ 
ature  distribution  along  the 
width  of  the  wall  in  a  steady- 
state  thermal  process,  l)  Boun¬ 
dary  layer;  2)  insulation; 

3)  plating 

The  temperature  of  the  outer  plating  surface  will  be 

(3.6) 

Similarly,  using  Formulas  (3*5)  and  (3.6),  we  can  determine  the 
temperature  of  the  inner  plating  surface. 

§3.3.  THE  TEMPERATURE  OP  THE  WALL  DURINO  SOLAR  HEATING 

At  altitudes  higher  than  15O  km  the  aerodynamic  heating  becomes 
insignificant  and  may  be  neglected;  the  equation  of  heat  balance  in 
the  case  of  a  steady-state  thermal  process  will  take  the  following  form: 

(3.7) 

If  we  do  not  take  into  consideration  the  heat  flows  along  the 
surface,  l.e. ,  assvime  that  =  0,  then,  expanding  the  value  of  q^^ 
from  Formula  (l. l4),  we  obtain: 


(3.8) 


For  space  craft  several  thousand  kilometers  away  from  the  earth, 


the  effect  of  terrestrial  radiation  and  the  solan?  rays  reflected  from 


the  earth  Is  insignificant  in  comparison  with  the  effect  of  direct 
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solar  rays,  and  consequently  the  specific  heat  flows  and  q^^  can  be 
neglected.  In  the  case  of  lower  altitudes,  the  effect  of  terrestrial 
radiation  and  solar  rays  reflected  from  the  earth  becomes  greater; 
however,  as  can  easily  be  seen  from  Formulas  (1.9),  (l. lO),  and  (3.8), 
this  effect  Is  still  quite  small  and  these  types  of  radiation  will 
increase  the  temperature  by  less  than  Consequently,  it  is  entirely 
permissible  In  a  number  of  cases  to  neglect  q^  and  q^^. 

If  we  neglect  q^  and  q^^,  and  take  into  consideration 
Formula  (1.8),  Formula  (3.8)  can  be  presented  In  the  following  form; 

7'ct=|/  (3-9) 

Of  particular  Interest  Is  the  temperature  for  f  =  0,  when  the 
temperature  Is  maximum;  In  this  case,  taking  Into  consideration  S  and 
a,  we  will  have: 


(f){T)  ■T- 

If  the  craft  rotates  rapidly  about  an  axis  perpendicular  or  al¬ 
most  perpendicular  to  the  direction  of  the  solar  rays,  and  if  the 
craft  has  a  sufficiently  thick  plating  with  good  heat  conduction,  the 
temperature  of  the  plating  may  be  Identical  over  the  entire  surface. 
The  same  will  be  true  when  there  is  heat  transfer  through  the  Inner 
wall.  This  average  temperature  may  be  found  from  the  condition  that 
the  total  flow  of  solar  radiation  falling  on  a  hemisphere  heats  the 
entire  surface  of  the  sphere. 

The  total  heat  flow  falling  on  the  sphere  and  absorbed  by  it  will 

be; 
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The  total  radiation  heat  flow 


4k/<’5971p, 

where  Is  the  mean  temperature  of  the  entire  surface. 

Prom  the  equality  of  these  heat  flows  we  will  find; 

The  value  of  P/e  in  Formula  (3. ll)  may  be  found  from  the  material 
given  in  the  §§  2. 2  and  2.  3- 
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Chapter  4 

THE  HEATING  OP  A  WALL  DURING  A  NONSTEADY -STATE  THERMAL  PROCESS 
§4. 1.  The  Temperature  of  a  Thin  Wall 


In  the  case  of  a  nonsteady-state  thermal  process,  the  temperature 
at  various  parts  of  the  body  may  change  with  time.  Generally,  to  de¬ 
termine  the  temperature  of  a  body  at  a  given  Instant  In  the  case  of  a 
nonsteady -state  thermal  process,  we  must  solve  simultaneously  Eqs. 
(1.13),  (1.15),  and  (1.16).  However,  this  solution  Is  possible  only 
In  certain  specific  cases  In  which  a  number  of  simplifying  assump¬ 
tions  are  made. 

Taking  Into  consideration  the  fact  that  generally  the  loss  of 
heat  along  the  surface  produced  by  the  temperature  gradient  Is  small 
In  comparison  with  the  heat  which  heats  the  body,  we  will  In  the  fu¬ 
ture  examine  the  one -dimensional  problems  for  a  flat  plate,  l.e. ,  let 
us  assume 


Let  us  also  assume  that  the  thermal-conductivity  coefficient  does 
not  depend  on  the  temperature  and,  consequently.  Is  constant  along  the 
width  of  the  plating.  This  Is  entirely  permissible  for  metals  and  non- 
porous  Insulation  materials  (laminar  plates,  ceramics).  If  we  tedce  the 
mean  value  of  thermal  conductivity. 

In  the  case  of  the  assumptions  made  Eqs.  (1.11),  (1. 13),  (1. 13)> 


and  (1.16)  may  be  represented  in  the  following  way; 
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(4.1) 

(4.2) 


where  a  Is  the  coefficient  of  thermal  diffusivity. 

Structures  with  thin  walls  made  of  metals  possessing  good  con¬ 
ductance  are  in  wide  use.  In  these  cases  the  temperature  difference 
across  the  width  of  the  plating  may  be  small  and  we  can  then  assume 
that  the  plating  is  suddenly  heated  along  the  entire  width.  In  §4.  3 
we  presented  a  criterion  characterizing  the  permissible  constant  temp- 
erture  along  the  width  of  the  plating,  1. e. ,  the  criterion  of  thin 
shielding. 

Assuming  for  a  thin  plating 

^0«=-7'cTt 

let  us  find  for  the  specific  heat  flow  which  heats  the  plating,  the 
expression 

In  this  case  we  can  substitute  the  found  expression  for  for 

the  first  term  of  Eq.  (4.1):  then 

Even  though  the  plating  radiates  from  two  sides,  in  the  case  of 
the  hollow  unfilled  casing  formed  by  the  plating  (wing,  nose  of  the 
body),  as  a  result  of  the  mutual  heat  transfer  between  the  inner  sec¬ 
tions  of  the  plating,  the  heat  will  not  be  able  to  escape  through  the 
inner  surface.  When  there  is  a  substantial  temperature  difference  bet¬ 
ween  the  lower  and  upper  surface,  the  Internal  radiation  of  the  hotter 
plating  section  will  be  greater  than  the  absorption.  This  may  be  taken 
into  consideration  by  increasing  the  radiation  coefficient  (emlsslvlty) 
determined  on  the  basis  of  a  first -approximation  calculation. 

In  the  presence  of  substantial  heat -absorbing  masses  Inside  the 
casing,  the  total  emlsslvlty  may  be  approximately  determined  from  the 
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following  formula 


where  Is  the  Inside  temperature  of  the  structure,  e  is  the  emissivity 
for  one  side  of  the  plating;  is  the  emissivity  of  the  inside  of  the 
structure. 

Since  a  and  in  the  case  of  nonsteady -state  heating  are  generally 
constant,  Eq.  (4.3)  should  be  solved  by  numerical  Integration  with 
finite  temperature  Increments 

It  should  be  borne  in  mind  that  in  the  case  of  a  nonsteady-state 
thermal  process  produced  by  aerodynamic  heating,  generally  can  be 
neglected. 

In  order  better  to  explain  the  method  for  determining  the  temper¬ 
ature  of  the  plating,  we  examine  below  an  example  for  the  determination 
of  the  plating  temperature  of  a  cone  with  a  flare  angle  of  20° ,  at 
a  =  0;  the  flight  altitude  is  15  km.  The  flight  regime  was  as  follows; 
the  flight  was  initially  stabilized  at  M  =  1.  5;  then  the  velocity  was 
uniformly  accelerated  for  20  seconds  up  to  M  =  4;  this  was  followed  by 
a  velocity  decrease  corresponding  to  passive  flight  according  to  the 
law 


1_ 

V 


where  Xp  is  the  passive-flight  time. 

Table  4. 1  gives  the  calculation  formulas  and  numerical  calculation 
of  the  flow  characteristics  at  the  surface  of  the  cone  (beyond  the 
limits  of  the  boundary  layer). 

The  temperature  determination  for  an  aluminum  alloy  plating  with 
a  thickness  of  5  mm  in  a  nonsteady -state  flight  regime  is  given  in 
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Table  4.2.  The  same  table  gives  the  calculation  formulas.  The  initial 
temperature  of  the  plating  (302°)  was  determined  from  the  steady -state 
thermal  process.  The  initial  temperature  in  the  given  case  is  close  to 
the  temperature  of  the  adiabatic  wall.  The  temperature  of  the  plating 
was  determined  for  a  point  at  a  distance  x  =  0. 1  m  from  the  nose  of 
the  cone.  The  emlssivlty  of  the  plating  was  taken  as  =  0.8,  the 
specific  heat  capacity  was  taken  as  c  =  0.224,  and  the  bulk  weight  was 
taken  as  y  =  2800. 


TABLE  4.1. 


Plating  Temperature  in  the  Case  of  Variable  Plight  Velocity 
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1.062 

1.090 

1,118 
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1)  Value  and  formula;  2)  number  of  formula;  3)  flight  time  in  seconds. 
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The  minimum  Reynolds  number  Re^  will  be  evident  at  the  beginning 
of  the  flight  when  the  velocity  Is  lowest.  In  this  case,  on  the  basis 
of  the  data  given  In  Table  4.1,  we  will  have: 


7I-.216, 5  1,052  =  228;  c  =  20K?7-302! 
Vi=302-1,38  -  416  m/sec,i,a.l,52.10-''; 


l,14S-0.01974  -  0.0226; 
0,77.10- 


5,4- 10*. 


TABLE  4.2. 


Determination  of  the  Cone -plating  Temperature 
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— 

P«- 0.01974 

— 

0,0182 

0,0188 

0,0190 

0.0189 

0,0183 

0,0170 

0,0161 

1 

0,0149 

0,0144 

0,0142 

— 
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— 
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— 
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— 
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(4.4) 

0 

6,3 

17,7 

44 

58 

38 

50 

16,5 

0.5 

-7 

1)  Value  and  formula;  2)  number  of  formula  3)  flight  time  In  seconds; 
4;  according  to  Pig.  1. 2. 
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fC 

500 

>iOO 

300 

200 

WO 

0 

Pig.  4. 1.  Plating  temperature  in 
the  case  of  variable  flight  veloc¬ 
ity.  (plating  made  of  aluminiun 
alloy,  6=5  nim).  1)  Temperature 
of  the  adiabatic  wail;  2)  temper¬ 
ature  of  the  plating. 

In  the  case  of  the  obtained  Reynolds  number,  the  boundary  layer 
will  be  turbulent.  Since,  as  the  flight  progresses,  the  Reynolds  number 
increases,  it  follows  that  the  boundary  layer  will  be  turbulent  for  the 
entire  theoretical  flight  time. 

The  results  of  the  calculation,  given  in  Table  4.2,  are  repre¬ 
sented  graphically  in  Pig.  4. 1.  It  is  characteristic  that  the  temper¬ 
ature  rise  for  the  plating  lags  behind  the  temperature  rise  for  the 
adiabatic  wall.  IHiis  property  is  called  thermal  inertia  and  is  deter¬ 
mined  by  the  heat  capacity  of  the  plating.  The  thermal  inertia  enables 
us,  during  short-term  heating,  to  use  materials  less  heat  resistant 
than  those  which  are  required  for  extended  heating. 

In  the  case  of  flight  at  high  altitudes,  when  the  aerodynamic 
heating  is  negligible  or  completely  absent,  Eq.  (4.3)  is  simplified: 


^  «» 
dt 


(4.5) 


If  the  magnitude  of  the  absorbed  radiant  energy  q^^  is  constant 
with  respect  to  time,  this  equation  is  easy  to  integrate.  Given  the 
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variable  q^,  Eq.  (4.5)  should  be  solved  by  numerical  integration. 

§4.2.  TEMPERATURE  DISTRIBUTION  ALONG  THE  THICKNESS  OF  THE  PLATING 

When  the  thermal  conductivity  of  the  plating,  or  at  least  a 
major  part  of  It,  Is  small,  the  method  for  determining  the  temperature 
of  the  plating  indicated  in  §4.1  may  lead  to  very  high  errors.  Ac¬ 
tually,  we  do  sometimes  encounter  conditions  of  aerodynamic  heating 
in  which  the  inner  surface  of  the  plating  and  the  adjoining  layers  of 
materials  cannot  be  heated  substantially  when,  at  the  same  time,  the 
outer  surface  may  already  have  a  high  temperature  and  may  even  be  melt¬ 
ing. 

To  explain  the  heat -shielding,  heat -resistance,  and  strength  prop¬ 
erties  of  the  plating  it  is  necessary  to  know  the  temperature  distri¬ 
bution  along  its  thickness,  while  taking  into  consideration  its  thermal 
conductivity  and  heat  capacity.  To  do  this  we  solve  the  Fourier  equa¬ 
tion  (4.2); 

IL^a  — 

dx  dy*  ■ 

The  following  boundary  conditions  are  assumed; 

l)  there  is  no  heat  transfer  on  the  inner  wall;  consequently,  the 
temperature  gradient  on  the  inner  wall  should  be  zero; 


2)  the  heat  transfer  on  the  outer  wall  obeys  Newton's  law;  the 
heat-transfer  coefficient  and  recovery  temperature  are  constant.  In 
this  case  we  will  have; 


(4.7) 

This  equation  is  a  simplified  version  of  Eq.  (4. 1)  in  which  the 
radiated  and  absorbed  radiant  energy  are  assumed  to  be  zero; 
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3)  the  initial  surface  temperature  of  the  plating  is,  in  general. 


constant  along  the  thickness. 


Pig.  4. 2.  Diagram  of  tempera¬ 
ture  distribution  along  the 
thickness  of  the  boundary- 
plating.  l)  Boundary  layer: 

2)  plating. 


In  the  integration  of  the  Fourier  equation  we  will,  for  the  most 
part,  observe  the  method  Indicated  by  A.  V.  Lykov  [l4].  Let  us  intro¬ 
duce  the  dimensionless  relative  temperature 


T,-T 
r,-ro  ’ 


(4.8) 


where  T  is  the  temperature  in  any  layer  of  the  plating.  Let  us  also 
Introduce  the  relative  coordinate 


(4.9) 


For  symbol  denotation  see  Pig.  4.2.  The  Fourier  equation  (4.2) 
for  the  relative  temperature  and  relative  coordinate  will  be; 

at  a  att 

(4.10) 

Let  us  find  the  particular  integral  of  this  equation  in  the  form 


where  X  is  only  a  function  of  if,  while  Z  is  only  a  function  of  t,  1.  e. , 
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Substituting  the  expression  for  S  in  Eq.  (4.10),  we  get 

dZ  a  d"iX  V 

Since  the  right  and  left  sides  of  this  equation  are  Independent 
(they  are  functions  of  various  Independent  arguments),  evidently,  each 
side  should  equal  the  constant  B,  1. e. , 

dZ  •wo  ^  'Y  D 


Let  us  Integrate  the  first  equation.  Designating  the  value  of  Z 

(4.11) 


as  Zq  for  T  =0,  we  obtain 


Let  us  designate 


Z=Zoe»’. 


Then  the  second  equation  will  be 


d\‘  ^ 


The  Integral  of  this  equation 

■X  s=  Cfg*  ^ 

Condition  (4.6)  gives; 


(4)  “0. 

V  d(  /*-« 


consequently. 


C,  =  f|. 

Let  us  designate  the  value  of  X  as  for  ^  =  1,  then 
Z,=»2f| — - =2tf,coi*, 

whence 


consequently. 


c,-±-X. 

'  2  CM* 


X=»X, 


eoi(t() 

CM* 
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The  value  of  0  may  now  be  expressed  In  the  following  form 

6  I -/?  » 


0 = ( cB  *  COS  (k-)  =  0„o  , 

\cQsJi  /  '  cos  * 

where  It  is  denoted  that 


^1^0  —  ®ciO' 

Here  Is  the  relative  temperature  of  the  outer  surface  for 

T  =  0. 

Let  us  express  the  value  of  B  In  terms  o'f  k;  then 


The  expression  in  this  equation 


at 

®o6 


(4.12) 


is  the  Fourier  similarity  criterion.  This  criterion  is  the  fundamental 
quantity  determining  the  temperature  of  a  body  in  a  nonsteady -state 
thermal  process. 

Expressing  0  in  terms  of  k  and  the  Fourier  criterion,  for  the 
particular  integral  of  Eq.  (4. lO) ,  we  will  obtain 


®*=®cT0  cos  (A!;). 


COS  k 


(4. 13) 


The  value  of  k  in  Eq.  (4.13)  is  determined  from  the  boundary  con¬ 
dition  (4.7).  Using  Expressions  (4.8)  and  (4.9),  we  represent  this 
condition  in  the  form 


'^1  --f 


where 


while  p  is  the  Blot  number 


a  ^r~^CT 


(4. 14) 


^o6 

X  ' 


(4. 15) 


Differentiating  (4.13)  with  respect  to 


di  *"®-  cos* 
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Substituting  this  equation  Into  Expression  (4. l4)  and  taking  into 
consideration  the  fact  that  when  i  =  1,  6  -  we  will  derive  the 

following  equation  In  which  Eq.  (4.13)  will  have  been  taken  into  con¬ 
sideration; 


This  equation  is  useful  for  determining  k.  Since  tan  k  is  a  periodic 
function  with  the  period  it,  k  will  be  a  multivalued  function.  Equation 
(4. l6)  is  transcendental  and  cannot  be  exactly  solved  analytically.  We 
will  deal  with  the  practical  determination  of  k  below. 

The  total  Eq.  (4.10)  integral  can  be  a  linear  function  of  the  par¬ 
ticular  Integrals  given  by  Eq.  (4.13).  If  k^  is  a  value  of  k,  from  Eq. 
(4.16),  the  total  integral  of  Eq.  (4.13)  will  be: 

where  the  constants  ^n  included  in  the  coefficients  of 

the  series  A^. 

The  values  of  are  determined  from  the  initial  conditions  at 
T  =  0.  Let  0Q  =  f(i)  at  the  initial  instant  of  time;  then 


8« 


S  ^,cos(*,?)=/'a). 


Let  us  multiply  both  sides  of  this  equation  by  cos  (kjjj^)d^  and 
let  us  Integrate  the  new  equation  obtained  within  the  limits  of  change 
in  ^  from  0  to  1; 


cos(*^E)<«. 


(4. 18) 


Replacing  the  derivative  of  the  cosines  by  the  sum  of  the  cosines 
and  Integrating,  we  obtain; 

I  I 

/= j  cos  ( cos  (**,?)  ds  “  Y  J 


1  [sin  (*/,-»«)  . 

2  1  *«-*•  *«  +  *•  J' 
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Taking  Into  consideration  the  fact  that 


Where  =  +k^. 

Where  k^^  using  Function  (4.  l6),  we  will  obtain: 

irl'k^  ^>n~Ksink„ cos k„)  =. 

“*  ’iTI**"  *<■  ■“  COS  A,)  ==  0. 

It  in 

Consequently,  In  the  series  of  Expression  (4.18),  of  all  the 
terms  only  two  remain,  for  which  k^  =  k^^^  and  k^  =  Expression 

(4. l8)  will,  in  this  case,  be  written  In  the  following  way: 

0 

Hence 


4sinfe^  kn 
2M/I  -|"  sin  ft/f 


1 

j/(E)cos(A,0^5. 


When  “  ICq  only  one  term  remains  In  the  series  of  Expression 

(4. 18) ! 


whence 


^  4  sin  *0  ^0 


,/(?)cos(AoO<«. 


The  total  Integral  in  (4.1/)  may  now  be  represented  In  the 
following  form  ,  4si„*,  .s,  *  i 

0 

Here  we  take  Into  consideration  the  fact  that  k__^j  and  that  cos 
(kj^O  is  the  even  function.  This  is  the  final  expression  for  the  re¬ 
lative  temperature. 


For  the  sake  of  convenience  we  will  assign  the  subscript  ”1' 
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rather  than  "0"  to  the  first  term 


B  =— 

*  2*,+  sin2*«’ 


then 


(4.20) 

(4.21) 


(4.22) 

The  value  of  the  coefficients  B  should  be  determined  while  tak- 

n 

Ing  Into  consideration  Eq.  (4. 16)  which  determines  the  coefficients 

k^.  The  first  value  of  k^  corresponds  to  the  value  of  the  angle  In 

the  first  quarter  of  the  circumference;  the  second  value  corresponds 

to  the  angle  In  the  third  quarter,  etc.  Thus  the  range  of  the  values 

of  k  will  be; 
n 

0<A,<-^«;  2n<k^C±K 

The  values  of  corresponding  to  these  values  of  k^  will  agree 
with  Formula  (4. 20); 


1  <5, <1,273:  0>i5, >-0,425;  0<fi, <0,255. 

The  values  of  are  determined  for  the  most  part  by  the  functions 
f(0,  and  also  by  the  values  of  the  coefficients  k^.  At  constant  tem¬ 
perature  along  the  width  of  the  plating  at  the  beginning  of  heating, 
the  functions  0  =  f(0  =1  34id  “  1*  If  the  temperature  In  vaj?lous 
layers  of  the  plating  Is  higher  than  the  temperature  of  the  outer  sur¬ 
face  at  the  beginning  of  heating,  f(^)  <1  and,  consequently,  <  1. 
Table  4.3  gives  the  values  of  the  first  tliree  coefficients  k^^, 

while  Table  4.4  gives  the  values  of  the  coefficients  B„  obtained  from 

n 

the  work  of  A.  V.  Lykov  ([14],  pp,  155  and  l6o). 
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Since  ^  =  1  for  the  outside  surface, 


^  cos  *««“*"*.  (4.23) 

n*l 

For  the  Inside  surface  ^  =  0,  f(0  =  1,  and  =  1;  consequently, 

(4.24) 

In  the  majority  of  cases,  for  calculation  of  6  from  Formulas 
(4.22)  (4.23)  and  (4.24),  it  is  entirely  permissible  to  take  three 
terms  of  the  series,  since  the  series  for  6  converges  rapidly. 

Formula  (4.22)  can  be  used  to  calculate  the  temperature  distri¬ 
bution  along  the  thickness  of  the  plating  and  when  a  and  are  vari¬ 
able.  In  this  case  the  given  heating  or  cooling  times  for  the  plating 
'  are  subdivided  into  ranges  such  that  within  these  ranges  the  change 

in  a  and  is  comparatively  small.  Within  these  ranges  it  is  assumed 
that  a  and  are  constants  and  Formula  (4.22)  is  used.  Let  us  note 
that  if  Tj,  is  variable,  6  at  the  end  of  the  preceding  range  does  not 
equal  the  value  of  0  in  the  beginning  of  the  subsequent  range.  This 
is  explained  by  the  step-wise  nature  of  T^. 

Formula  (4.22)  did  not  take  into  consideration  the  radiation  of 
V  the  plating.  It  should  be  noted  that  in  a  nonsteady -state  aerodynamic 

heating  regime  far  from  temperature  equlllbrlvun,  the  effect  of  the 
radiation  heat  on  the  temperature  regime  of  the  plating  is  not  great. 
Consequently,  in  this  case,  the  radiation  heat  may  be  calculated  ap¬ 
proximately  by  introducing  the  equivalent  heat-transfer  coefficient: 

®  ?iu/(  “*  ^t)  • 

Since  vau?lable,  the  calculation  taking  the 

radiation  into  consideration  should  be  ceu?rled  out  in  the  same  way  as 
that  of  the  variable  heat -transfer  coefficient,  l.e. ,  by  dividing  the 
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flight  time  Into  Intervals. 

At  a  wall  temperature  close  to  equilibrium,  the  radiation  effect 
may  be  taken  into  consideration  if  we  represent  the  specific  heat  flow 
involved  in  the  heating  of  the  plating  in  the  form 


^utr  —  *  (  ^ct)  (^plH  T'cr) > 

Where  is  the  equilibrium  temperature  and  determined  from  equa¬ 

tion 

Calculating  the  latter  formula  from  the  preceding  and  factoring 
’^ravn  ' 


For  materials  with  low  conductivity,  for  example,  heat -insulation 
materials,  after  a  short  heating  period,  gets  very  close  to  T^^vn 
and  consequently  in  the  latter  expression  we  may  assume  ^st/^ravn  “ 
then 

We  can  sometimes  be  limited  to  one  term  of  the  series  for  0. 
Actually,  in  this  case  when  p  <  0. 1,  the  second  term  of  the  series 
will  be  less  than  2^  of  the  first,  while  the  third  will  be  less  than 
o.  3^  of  the  first.  For  qp  >  0,  the  relative  values  and  the  second  and 
subsequent  terms  will  be  even  lower. 

For  p  <  0. 1,  we  can  assume  in  approximate  terms 

tg 

then  •»  P  and  k^  «  P. 

The  error  in  the  determination  of  k^^  under  this  assumption  will 
be  less  than  1.6$^.  In  this  case 
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at 

»o«eT 


t 


TABLE  4.3. 

Values  of  Coefficients  k 

n 


p 

0 

an  k„ 
0,001 

=  P) 
0,002  j 

0,00-4  j 

0.006 

0,008 

0,01 

0,02 

0.04 

0.06 

*1 

0.0000 

0,0316 

0,0632 

0,0774 

0,0893 

0,0998 

0,1410 

0,1987 

0,2425 

ItJ 

3»14I6 

3,1419 

mm 

3,1429 

3,1435 

3,1441 

3,1448 

3,1479 

3.1543 

3,1606 

*3 

6.2832 

6,2833 

6,2838 

6,2841 

6,2845 

6,2848 

6,2864 

6,2895 

6,2927 

P 

0,08 

0,1 

m 

■1 

0,5 

0,6 

0,7 

0,8 

0,9 

0,2791 

n 

0,4328 

0,5218 

0,5932 

0.6533 

0,7506 

0,7910 

0,8274 

3,1668 

3,2039 

3,2341 

3.2636 

3,2923 

3,3477 

3,3744 

.  *3 

6,2959 

B 

6,3148 

6,33C5 

6,3461 

6.3616 

6,3770 

6,3923 

6,4074 

6,4224 

p  ! 

1.0 

1,5 

2.0 

IB 

m[ 

5,0  1 

6,0 

jggH 

8.0 

9.0 

0.8603 

0,9882 

1,0769 

1,1925 

1,3138 

1,3496 

HB 

m 

fl2 

3,4256 

3,5422 

3,6436 

3,9352 

4,1116 

■Bl 

^3 

6,4373 

6,5097 

6,5783 

B 

6,8140 

6,9096 

6,9924 

7,1263 

Bfl 

P 

10 

15 

20 

30 

40 

60  ■ 

60 

80 

100 

oo 

*« 

1,4269 

1,4729 

1,4961 

1.5202 

1,5325 

1,5400 

1,5451 

1,5514 

1,5552 

1,5708 

*! 

4,3058 

4,4255 

4,4915 

4,5615 

4,5979 

4,6202 

4.6353 

1,6543 

4,6658 

4,7124 

7,2281 

7,3959 

7,4954 

7,6057 

7,6647 

7,7012* 

7,7259 

7,7573 

7,7764 

7,8540 

Where  we  designate 

(Xg  is  the  time  constant  of  thermal  inertia). 

When  p  <  0. 1,  we  may  also  simplify  Expression  (4.20)  for  B^. 
Assuming  »  sin  k^^,  we  obtain; 

B  — J— — 2_ 

‘  1  +  cos  A|  1  -{•  CO*  yif  ' 


With  the  assxunption  made,  Q  may  be  represented  in  the  form 


i+co*/p  1^.25; 

In  the  case  of  equilibrium  temperature  distribution  along  the 
thickness  of  the  plating,  when  r  «  0,  (;^  <=  1,  and  consequently, 

where  3  is  determined  from  Formula  (4.15). 


(4.26) 
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TABLE  4.4. 


Values  of  Coefficients 


(B/. 


4$in  k„  _ 
2k„  +  sin  24,  ’ 


<g  kn  ^ 


p 

0 

0,001 

0,002 

0,004 

0,006 

0,008 

0.01 

0,02 

0,04 

0,06 

Bi 

1,0000 

■n 

■■■ 

1,0012 

1,0015 

1,0016 

1,0030 

■n 

Bj 

0,0000 

—0,0012 

-0,0016 

—0,0020 

-0,0040 

Ba 

0,0000 

0.0000 

0,0002 

0,0003 

0,0004 

0,0005 

0,0010 

BH 

P 

0,08 

0,10 

0,20 

0,30 

0,40 

0,50 

0,60 

0.70 

0,80 

0,90 

Bi 

1,0159 

1,0312 

1,0581 

1,0701 

1.0813 

BBI 

m 

1,1107 

Bi 

SB 

-0,0197 

-0,0381 

—0,0719 

-0,1025 

—0.1493 

Bj 

■ 

0,0050 

0,0100 

0.0148 

0.0196 

0,0289 

0,0335 

IH 

0.0423 

P  . 

1.0 

.  1.5 

2,0 

3,0 

4,0 

!  5.0 

1 

6,0 

8.0 

9,0 

Bi 

i.n92 

I 

1  1,1537 

1,1784 

1,2102 

1.2287 

BSI 

1 

1.2532 

1.2569 

1,2598 

Bi 

-0,1517 

-0,2013 

-0,2367 

—0,2881 

—0,3215 

—0.3722 

'  -0.3880 

Bi 

0,0466 

0.0667 

0,0848 

0,1154 

0,1396 

0,1588 

0,1861 

0,1959 

0,2039 

P 

10 

15 

20 

30 

40 

50 

60 

80 

100 

GO 

1.2612 

1,2677 

1.2699 

1,2717 

1,2723 

1,2727 

1.2728 

1,2730 

1 

1,2731 

.  1,2732 

Bi 

—0,3934 

-0,4084 

-0,4147 

—0,4198 

-0,4217 

•—0,4227 

-0,4232 

—0,4237 

-0.4239 

-0.4244 

B, 

0,2104 

0,2320 

0,2394 

0,2472 

0,2502 

0,2517 

0,2526 

0,2535 

0,2539 

0.2546 

To  evaluate  the  accuracy  of  Pomrula  (4.26),  Table  4.5  gives  the 
errors  which  result  from  this  formula  In  comparison  with  the  exact 
Formula  (4.22)  with  a  constant  Initial  temperature  along  the  thick¬ 
ness  of  the  plating. 

It  Is  apparent  from  Table  4.5  that  the  accuracy  of  Poimila  (4.26) 

Is  determined  not  only  by  the  value  of  but  also  by  the  relative 

heating  time  t/t„.  At  t/t_  =  0,  Ponmila  (4.26)  produces  the  greatest 

error.  However,  these  errors  are  substantially  reduced  when  T/Xg  =  0.01. 

Formula  (4.26)  may  also  be  used  for  P  >  0.1,  If  In  this  case, 

t/t  Is  sufficiently  large.  For  example.  In  the  case  of  P  *0.2  and 
s 

t/t  =0.1,  there  are  few  errors. 

'  s 

It  is  apparent  that  the  accuracy  of  Formula  (4.25)  will  be  approx- 
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imately  the  same  as  the  accuracy  of  Formula  (4.26). 


TABLE  4.5 

Accuracy  of  Formula  (4.26) 


"~H-eosKf^  ’ 


t 

p 

li.,.  (5=0) 

Bcr  (5=1) 

no 

(4.26) 

no 

(4.22) 

louufdKa 

B 

no 

(4.26) 

no 

(4.22) 

ouinOKa 

B  •/• 

0 

0,05 

1,012 

1.000 

+U2 

0,989 

1,000 

~1,1 

0 

0,10 

1,023 

1,000 

+2,3 

0,975 

1,000 

--2,5 

0 

0,20 

1.050 

1,000 

h5,0 

0,947 

1,000 

-5.3 

0.01 

0,05 

1,002 

0.997 

+0,5 

0,977 

0,976 

+0,1 

0,01 

0,10 

1,013 

0,999 

+1,4 

0,963 

0,965 

-0,2 

0.01 

0.20 

1,040 

1,000 

+3,6 

0,9.38 

0,953 

-1.6 

0,10 

0,05 

0,918 

0,917 

+0,1 

0.895 

0,895 

0 

0,10 

0.10 

_o;'928 

0,921 

+0,8 

0.883 

0,879  i 

+0,5 

0,10 

0,20 

0,950 

0,910 

+J.1 

0,860 

0,854  ' 

+0,7 

1 )  Error,  In  % 

In  the  case  metal  plating,  P  in  many  practical  instances  of  aero¬ 
dynamic  heating  will  be  less  than  0.1  and  consequently  to  reduce  the 
amount  of  calculation  while  determining  Q,  we  may  use  Formula  (4.25) 
and  (4.26)  In  these  cases. 

§4.3.  THE  CRITERION  OF  THIN  PLATING 

§4.1  gave  a  solution  of  the  heat -conduction  equation  on  the 
assumption  that  the  plating  was  so  thin  that  the  temperature  differ¬ 
ence  across  its  thickness  was  insubstantial.  The  question  naturally 
arises,  when  can  plating  be  considered  thin;  In  other  words,  we  must 
establish  a  criterion  for  thin  plating. 

Since  the  temperature  difference  across  the  thickness  of  the 
plating  never  equals  zero,  for  all  practical  purposes  It  is  necessary 
to  establish  the  permissible  difference  at  which  the  plating  may  be 
considered  to  be  heating  uniformly.  We  will  characterize  the  heating 
equilibrium  of  the  plating  by  the  relation  of  the  temperature  incre¬ 
ment  of  the  inner  plating  surface  to  the  temperature  Increment  of  the 
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of  the  outer  surface,  1.  e.  , 


Let  us  Introduce  the  concept  of  a  relative 
across  the  thickness  of  the  plate  and  designate 

,  4r„ 


&Tc 


1- 


'CT 

Let  us  express  the  relative  temperature  of 
temperature  Increment,  using  Formula  (4.8); 

r,-/-,  Tr-To  ’ 


temperature  difference 
it  as  Sp,  in  which  case 

(4.27) 

6  in  terras  of  the 


where  AT  =  T  —  Tq, 
hence 


4r=(i-o)(r,-ro). 


Consequently, 


e  =  1  •—  ^»w)  „  (9,„— -j^c,) 

"  (1-ec,)  (1 -»„)'* 


(4.28) 


Assuming  the  values  of  the  permissible  relative  temperature 
difference  e  ,  from  the  latter  equation  we  can  find  the  plating  thick- 
ness  at  which  this  temperature  difference  may  be  obtained.  Evidently, 
thicknesses  less  or  equal  to  that  found  will  correspond  to  thin  plat¬ 
ing. 


Since  it  is  necessary  to  have  an  approximate  solution  for  the 
criterion  of  thin  plating,  we  can  derive  it  from  approximate  functions. 
For  thin  shielding,  the  Fourier  criterion,  as  can  be  seen  from  what 
follows,  is  qp  >  1.  Consequently,  in  the  expression  for  the  relative 
temperature  (4.22)  It  is  entirely  possible  to  limit  ourselves  to  one 

term  of  the  esqpanslon.  For  low  values  of  S,  as  was  shown  in  §4.2,  we 

2 

may  assume  that  k^  ■■  In  addition,  the  value  of  the  coefficient 
according  to  Formula  (4.20)  may  be  assumed  to  be 
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Then^  with  the  asstunptlon  that  Porraula  (4.23)  and  (4.24) 

may  be  represented  in  the  following  form 


The  relative  temperature  difference  may  now  be  represented  In  the 

form 

e  —  ~cos  i/'p 
"  ]  —  cos  y  p'e~'f* 

At  low  values  of  3  we  can  expand  the  functions  cos  and  e“^^ 

In  series  and  limit  ourselves  to  two  term  of  the  expansion,  1. e. ,  we 
may  assume 


008  1/^=1  — 


Then 


..1 


1  l-py 


Hence 


1  —  en 


1  —  en 


%  i  ti-anP 


At  low  values  of  p  we  may  assume 


then 


9 


2  », 


(4. 29) 


(4. 30) 


Taking  Into  consideration  p  with  respect  to  Ponnula  (4.12),  we 
obtain: 


Consequently,  we  may  assume  thin  plating  to  be  that  plating  at 
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which 


«o«<  |/ 2aT .  (4.31) 

It  should  be  borne  in  mind  that  this  criterion  of  thin  plating 
is  applicable  in  cases  in  which  the  values  of  3  are  low,  i.e. ,  those 
values  at  which  the  above  approximations  for  cos  and  e“^l^ 
are  permissible,  and  this  with  a  sufficient  degree  of  accuracy. 

With  a  view  to  comparing  Formula  (4.30)  with  the  accurate  solu¬ 
tion,  Pig.  4.3  gives  the  graph  constructed  from  the  results  of  exact 
calculation  for  9  with  respect  to  P  where  =  0.1.  It  is  appararent 
from  this  graph  that  for  ^  <  0.1,  the  value  of  qp  from  Formula  (4.30) 
is  sufficiently  close  to  the  exact  value.  At  high  values  of  P,  the 
value  of  qp  drops,  so  that  for  P  >  10  we  can  assume  qp  w  1  for  Cp  =  0.1. 
In  this  case; 

(4.32) 

In  the  general  case 


8o8< 


(4. 33) 


This  formula  may  be  used  together  with  the  graph  given  in  Pig. 
4. 3>  from  which  we  determine  the  value  of  qp  at  Sp  =  0.1. 

When  P  >  10,  qp  may  be  determined  from  the  following  formula 


,_0,9331g(iHi).  (4.3,) 

The  latter  formula  was  found  from  the  condition  p  =  ». 

For  metals,  P  is  generally  of  the  order  of  0.1  and  less;  conse¬ 
quently,  we  may  use  Formula  (4.31)  as  a  criterion  of  thin  metal  plat¬ 
ing.  For  laminated  plates,  the  thermal  conductivity  is  less  by  a 
factor  of  a  hundred  than  in  the  case  of  metals;  consequently,  p  will 
be  substantially  lower  in  the  case  of  metals,  and  so  for  laminated 
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Pig.  4. 3*  Graph  of  the 
Fourier  criterion  for 
thin  plating  with  a  re¬ 
lative  temperature  dif¬ 
ference  of  e  =0. 1.  l) 
P 

Prom  the  following  for¬ 
mula 

<P  =  (1  -  ep)/2epi 

2)  from  exact  calcula¬ 
tion. 


plates  and,  in  general,  for  plastic,  we 
should  use  Formulas  (4,33)  and  (4.3^)* 

If  the  heat-transfer  coefficient  is 
small  (a  <  0.1),  then  in  the  case  of  plas¬ 
tic  the  criterion  of  thin  plating  should  be 
determined  from  Formula  (4, 3l)» 

As  far  as  ceramic  materials  are  concerned, 
their  thermal  conductivity  is  lower  by  a  fac¬ 
tor  of  only  ten  than  in  the  case  of  metals, 
and  consequently  the  value  of  qp  will  be  some¬ 
where  between  the  values  for  metals  and  lam¬ 
inated  plates.  These  values  of  9  should  be 
determined  from  graphs  analogous  to  that  given 
in  Pig.  4. 3^  or,  if  we  want  to  be  especially 


cautious,  we  should  calculate  6^^  from  Formula  (4.3l)» 

Table  4.6  gives  the  maximum  thickness  values  of  a  thin  plating  for 

various  materials  where  =  0. 1.  Metal  and  aluminum-oxide  calculations 

P 

were  carried  out  in  accordance  with  Formula  (4. 3l)>  while  asbestos, 
glass,  and  laminated  plastic  calculations  were  carried  out  in  accord¬ 
ance  with  Formula  (4.32). 

The  thermal  conductivity  of  metals  in  Table  4.6  was  determined  at 


200"  C. 

§4.4.  THE  DETERMINATION  OP  THE  THICKNESS  FOR  HEAT -INSULATION  PLATING 
A  nvimber  of  materials  with  low  thermal  conductivity  are  used  for 
heat  insulation.  The  required  heat-insulation  thickness  in  a  steady- 
state  thermal  process  is  easy  to  find  from  Formula  (3.4)  in  the  pre¬ 
sence  of  internal  cooling; 


= 


9oi 


(^•35) 
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TABLE  4.6. 

Maximum  Thickness  Values  of  Thin  Plating  for  Various  Materials  at 

ep  »  0. 1. 


MarepHa;!  ^ 

TcMnepatypo- 
np0B04H0CTb2 
mVcck  (251 

npeaeabHaa  BeaHHHHa 

TOHKOfl  oOuiHBXH  B  MM 

_ \ _ _ 1 _ _ 

npH 

TsslO  CCK. 

cfr- 

npH 

t as  100  ceK. 

^  A4iomhhmA 

8,7. 10-5 

14 

44 

MarHMfl 

5epH44Hfi 

7,2. 10-5 

13- 

41 

5.9- 10-5 

11 

36 

8  >Kejie30 

1,5. 10-5 

5.8 

18 

9  Cia^b  IX18H9T 

0,42. 1 0-5 

3 

9.6 

lOTuraH 

0,6- 10-5 

3,6 

11 

HM04H64eH 

5,  MO-5 

11 

34 

L2  HHKejib 

1,4.10-5 

5,6 

18 

-3  OKHCb  AloOa  npH  1000®  C 

6,210-’ 

1,2 

3.7 

Ik 

Tckcto.iht  npM  20®  C 

1,5.10-’ 

1.2 

3,8 

^^CicKao  npH  200®  C 

5,9.10-’ 

2.4 

7,6 

L6  CTeKJlOTeKCTOJIHT  (6) 

2. 10-’ 

1,4 

4,4 

ITAcCccTp  7=100  KtjM^ 

7.10-’ 

2.6 

8,4 

l)  Material;  2)  thermal  diffusivity;  3)  maximum  size  of  thin  plating, 
in  mm;  4)  where;  5)  aluminum;  6)  magnesium;  7)  beryllium;  8)  iron; 

9)  steel  lKhl8N9T;  lO)  titanium;  11;  molybdenum;  12)  nickel;  13)  alu¬ 
minum  oxide  AlgOg  at  1000°  C  [51] i  14)  textollte  at  20°  C;  15)  glass 

at  200°  C;  16)  glass -textollte  [6];  37)  asbestcc,  -y  =  100  kg/m^. 

The  required  insulation  thickness  is  proportional  to  the  thermal 
conductivity  and  does  not  depend  on  its  heat  capacity. 

Porous  and  loose  fibrous  materials  exhibit  the  lowest  thermal 
conductivity.  Moreover,  as  the  porosity  and  brittleness  Increase  and 
the  weight  of  the  material  correspondingly  decreases,  the  thennal  con¬ 
ductivity  of  the  material  approaches  the  thermal  conductivity  of  air. 
The  thermal  conductivity  of  asbestos,  depending  on  the  degree  of  its 
looseness,  will  be  (see  [25]); 

1  iipii  T=  500 .  x=2,5- 10-5  01+0,0018  0 

.  T=^00  . X  =  l, 8- 10-5  (1+0,0022  /) 

.  I^^IOO  . X=  1,5. 10-5  (1+0,0027  0 

1)  Where. 

The  thennal  conductivity  of  air 
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X  0.56  •  1 0-f-(  1-1-0,0030/). 


The  porous  and  loose  fibrous  materials  have  very  low  mechanical 
properties  and  consequently  can  be  used  only  for  the  Inner  facing  of 
the  plating. 

Among  the  heat -Insulation  materials  of  great  strength  are  certain 
types  of  plastics  and.  In  particular,  laminated  plastics.  However, 
their  thermal  conductivity  Is  greater  than  In  the  case  of  porous  ma¬ 
terials.  For  example,  textollte  has  a  thermal  conductivity  (see  [25]) 

\  «  7  •  10"“^  kcal/m  •  sec  •  deg. 

The  ceramic  materials  are  heat  resistant  and  sufficiently  strong 
to  be  used  for  the  outside  coating;  however,  they  have  even  greater 
thermal  conductivity.  For  example,  aluminum  oxide  AlgO^  has,  at  1000®  C, 
a  thermal  conductivity  (see  [51]) 

X.  =*  6  •  10  kcal/m  •  sec  •  deg. 

In  addition,  ceramic  Insulation  Is  very  brittle  and  cracks  easily 
under  Intense  heating. 

In  a  nonsteady-state  thermal  process,  the  required  heat-lnsulatlon 
thickness  should  be  determined  from  Eq.  (4.24)  In  the  general  case; 
moreover.  It  Is  necessary  to  give  the  permissible  temperature  values 
of  the  Inner  surface.  However,  the  solution  of  this  equation  for  tp, 
and  later  also  for  6^^,  In  the  general  case  can  be  obtained  only  graph¬ 
ically  or  by  Interpolation,  If  the  values  of  q)  are  given. 

The  temperature  of  the  Inner  surface  can  be  assumed  on  the  basis 
of  the  relative  temperature  difference  [see  Formula  (4.27)].  For  ex¬ 
ample,  for  practical  purposes  In  a  niunber  of  cases  the  following  would 
be  entirely  acceptable. 


then 
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we  can  de- 


Glven  the  value  of  and  knowing  the  value  of  6 

p  ^  I 

termlne  the  relative  temperature  of  the  Inner  surface  from  Formula 
(4.28): 

0«  =  *n  +  9ct(l— .). 

For  the  given  equation  we  can  limit  ourselves  to  one  term  of  the 
series  In  Eq.  (4.23)  and  (4.24). 

Then  at  =  1  (the  constant  temperature  along  the  thickness  of 
the  plating  at  the  beginning  of  Its  heating) 

®ct/9,k  =  COS  a,. 

The  calculation  of  exact  values  of  and  subsequent  deter¬ 

mination  of  the  corresponding  values  of  qp  by  successive  approximations 

Indicates  that  In  the  case  of  the  given  -  cos  k, ,  the  errors 

St'  vn 

will  be  at  high  values  of  e  .  In  reality^  when  e  Increases,  qp  decreases 
and  where  e  =  1,  qp  =  0.  This  Is  evident  from  the  structure  of  Eqs. 
(4.23)  and  4.24),  and  also  from  Table  4.7. 

Assiunlng  that  for  all  practical  purposes  we  scarcely  need  a  value 
of  Ep  greater  than  0. 95>  we  can  calculate  the  errors  for  »  0.95 
at  various  values  of  p.  The  relative  errors  In  the  determination  of 
the  Insulation  thickness,  assiuning  ^1  follows: 


1  Stia^rciiiiR  ^ 

0,1 

1,0 

10 

100 

Oltlll^KII  B  Ollpe;ie.1CtllfM  TO;ilUHHU  H30- 

4 

5 

0.7 

0.4 

.iimnii  B 

l)  The  value  of  P;  2)  errors  In  the  determination  of  the  Insulation 
thickness.  In 

Where  <  0.95,  there  will  be  even  fewer  errors  In  the  required 
Insulation  thickness. 

Assuming  ^g^/^yn  *  °°®  ^1*  above-given  value  for  6^, 
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we  find: 


(1— (1— s„)cos*il 


(4.36) 


Giving  the  value  of  we  can  determine  for  the  given  value 
of  p  from  this  formula. 

The  Insulation  thickness  required  to  ensure  the  given  relative 
temperature  difference  can  be  determined.  If  we  know  the  value  of  the 
Fourier  criterion; 
then 


■-/f- 


f  (4.37) 

The  value  of  q)  can  be  calculated  In  the  following  way.  Assuming 
certain  values  of  9,  let  us  find  the  corresponding  values  of  6^  from 
Eq.  (4.24).  Having  determined  the  required  value  of  0^^  from  Formula 

(4.36) ,  we  can  find  the  sought  value  of  9  graphically  or  by  linear 
Interpolation. 

Table  4. 7  gives  the  determination  of  values  of  9  for  various 

values  of  and  P.  We  were  limited  to  two  terms  of  the  series  In  the 
P 

first  three  rows  when  we  determined  the  values  of  0^^  from  Eq.  (4.24). 

In  the  subsequent  rows,  the  values  of  0^  were  determined  from  Formula 

(4.36) . 

The  value  of  9  was  determined  by  linear  Interpolation;  for  P  =  0. 01 
and  0.1,  the  value  of  9  was  determined  graphically. 

The  results  of  the  calculations  are  presented  graphically  In  Fig. 
4.4.  By  using  this  graph,  we  can  determine  rather  rapidly  the  required 
Insulation  thickness. 

It  should  be  noted  that  the  Indicated  method  for  the  determination 
of  the  Insulation  thickness  was  used  formally  under  those  conditions 
In  which  Eq.  (4.24)  was  found.  If  the  values  of  T^,  and  a  are  variable 
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In  flight.  It  Is  possible  approximately  to  determine  the  Insulation 

thickness  by  averaging  T^,  and  a  with  respect  to  time. 

TABLE  4.7. 

The  Determination  of  the  Fourier  Criterion  (qp)  at  Various  Relative 
Temperature  Differences 


p 

0»0l 

O.l 

1.0 

10 

100 

oa 

®BK  ¥^0,1 

0,99987 

0,9989 

0,9919 

0,9690 

0.9496 

0,9495 

Abh  £tt/  y=0i2 

0.99935 

0,9939 

0.9623 

0,8295 

0,7795 

0,7720 

^Bil  ^-=^0^3 

0,99851 

0,9850 

0,8915 

0,6745 

0,6144 

0,6065 

cos  kx 

0,995 

0.952 

0,652 

0,1413 

0.0158 

0 

s„=0.80  ( 

0,99880 
^  1 

0,9877 

1 

0,9199 

0,8230 

0,8030 

0,8000 

0,272  j 

0.275  , 

0,253 

0,204 

0,186 

0,184. 

s„=0,85  ( 

0,99912 

0,99125 

0,9423  ! 

0.8685 

0,8522 

0.8500 

1  9  1 

0.225 

0,237  1 

0.216 

0.172 

0.157 

0,156 

^.=0,90  (  ““ 

0,99950 

0,99470 

0.9628 

0,9130 

0,9015 

0,9000 

1  9 

0,175 

0,184 

0,174 

0,140 

0,1283 

0,1280 

e„=0,95  ( 

0,99974 

0,99748 

0,9820 

0,9567 

0.9506 

0,9500 

\  ? 

0,130 

0.135 

0,125 

0,109 

0,099 

0.097 

To  Illustrate  the  method  for  determining  the  required  Insulation 
thickness,  as  an  example,  let  us  determine  the  required  glass -textollte 
thickness  for  the  plating  of  the  Instrument  bay  of  a  cone  subject  to 
Intense  aerodynamic  heating  for  sixty  seconds.  Let  us  determine  this 
thickness  for  the  Initial  data  corresponding  to  the  example  considered 
In  §4.1. 

It  Is  necessary  to  assure  a  temperature  for  the  Inner  plating 
surface  that  Is  not  higher  than  80®  C. 

Let  us  Initially  determine  the  temperature  of  the  outer  surface 
which,  due  to  the  low  thermal  conductivity  of  glass -textollte,  can 
approximately  be  assiuned  to  be  equal  to  the  Instantaneous  equlllbrliun 
t«nperature,  understanding  by  the  latter  the  temperature  at  which 


Let  us  find  this  temperature  by  successive  approximations,  using 
Table  4. 2.  For  t  »  20  sec. : 
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1  33.aacMca  Tct  i 

1 

750 

780 

779 

778 

1 

2  Ha.NOallM 

6,71 

4,37 

1 

4.45 

4p52 

3  HaxoAHM  qn9 

3.85 

4,50 

4,48 

4,46 

l)  Given  Tg^;  2)  let  us  find  3)  let  us  find 

Hence,  It  Is  apparent  that  at  t  =  20  sec. ,  =  779°  K.  Similarly, 

let  us  find:  for  t  »  60  sec,,  =  490®  K. 

Let  us  average  the  wall  temperature  with  respect  to  time,  assum¬ 
ing  that  the  linear  change  of  Is  In  the  time  Intervals  from  0  to  20 
sec.  and  from  20  to  60  sec. : 

(7;r)cp=— ^•-65-  +  — y— K. 

Consequently,  the  temperature  Increase  for  the  outer  surface 

Ar„ = 604*  -  302’ = 302®; 

the  temperature  Increase  for  the  Inner  surface  according  to  the  condi¬ 
tion 

a4,==353‘’-302“  =  51®. 

The  relative  temperature  gradient 

s«l_®L= 0,831. 

"  302 

Averaging  ot  with  respect  to  the  temperature: 

_ 0y0323  0i0780  20  •  0,0/80  -I-  0,0465  40  _ q 

«cp=  ^  ^  60  ~  ' 

To  determine  the  value  of  p,  a  rough  estimate  of  the  plating  thick 
ness  =  5  mm  Is  assumed;  then 

0t0598*0,0Q5  »  2^ 

Prom  Pig.  4. 4  we  will  find  that  qp  ■  0. 202 
Consequently,  according  to  Pormula  (4.37) 
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0,0077  M. 


60 

0.202 


Let  US  make  the  second  approximation,  making  the  value  of  p  more  pre¬ 
cise: 


0,0598-0.0077  ,,  „ 

^ - - 6-57. 


Prom  Pig.  4.4  we  find  that  9  =  0.192,  and  consequently. 


8o.  =  j/  2.10-'.-^=.0,0079. 

It  is  apparent  from  the  calculation  that  the  effect  of  P  on  the 
result  Is  Insignificant.  It  should  be  noted  that  the  scale  of  the  axis 
of  abscisses  in  Pig.  4.4  Is  logarithmic  and  encompasses  a  range  of 
change  In  P  from  0.01  to  1.00  (by  factor  of  10,000).  Because  P  has 
little  Influence  on  the  result  of  the  calculation,  we  need  not  deter¬ 
mine  the  value  of  a  In  the  value  of  P  with  great  accuracy. 

The  heat -Insulation  calculation  given  Is  based  on  the  assumption 
that  the  temperature  of  the  outer  surface  Is  equal  to  the  Instantaneous 
equilibrium  temperature. 

In  fact,  the  temperature  of  the  outer  surface  and,  consequently, 
the  obtained  Insulation  thickness  Is  somewhat  Increased.  On  the  basis 
of  the  obtained  data,  we  can  determine  the  temperature  of  the  outer 
surface  and  make  more  precise  the  required  Insulation  thickness. 

Por  p  =  6.57,  from  Tables  4.3  and  4.4  we  find:  k^  »  1.36, 
kg  =  4.  l4,  =  1.25,  and  Bg  =  —0.366;  consequently  where  9  ■  0. 192 

from  Eq.  (4.23)  we  find: 

=  1.25*0.2108  +  0.366*0.5446  e“3*28  ^  0.1850  +  0.0075  = 

=  0. 1925. 


To  calculate  the  radiation  we  assume  that  the  temperature  of  the 


adiabatic  wall  Is  equal  to  the  equilibrium  temperature;  then  we  obtain: 


Tr-Tc 


604  —  7’, 


Tr-To 


604  —  302 
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H=0,1925; 


hence  Tg^  »  546°  K.  Consequently,  ATg^  »  546  -  302  =  244  and 


e„  =  l--^=l  =  0,799. 


Afc 


248 


Prom  Pig.  4.4  we  find  qp  =  213;  consequently 


3o6-j/ 2-10-‘. -^  =  0,0075  .«-7,5 


MM. 


Consequently,  the  temperature  taken  equal  to  the  equilibrium 
temperature  produced  an  error  of  +55^  In  the  determination  of  the  heat- 
insulatlon  thickness. 


Pig,  4,4,  Graph  of  the  Pourler 
criterion  at  various  relative 
temperature  differences. 

We  can  approximately  estimate  the  value  of  the  Pourler  criterion 
from  the  following  formula 

?=I.7lglii.  (4.38) 

According  to  this  fonmila,  the  error  In  the  determination  of  the 
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thickness  for  e  >0.8,  will  be  less  than  125^. 

P 

According  to  this  formula,  for  the  above-considered  example 


Consequently 


?  =  1.71g 


1.1 

0,798 


0,238, 


S„i5  =  0,0071  M, 


1. e. ,  the  error  Is  5^. 

When  selecting  the  heat-lnsulatlon  material,  given  satisfactory 
mechanical  strength  properties  and  physical  (adhesion,  breakdown  tem¬ 
perature)  properties  assuring  the  possibility  of  Its  use  under  the 
given  conditions.  It  Is  necessary  to  evaluate  the  feasibility  of 
using  a  given  material  with  respect  to  weight. 

Formula  (4.37)  makes  it  easy  for  us  to  compare  the  various  ma¬ 
terials  with  respect  to  weight  in  a  nonsteady -state  thermal  process. 

2 

In  reality,  the  weight  of  1  m  of  the  sheet  material  for  the  required 
heat-lnsulatlon  thickness  will  be: 

The  value  of  x/q)  Is  In  the  main  determined  by  external  conditions, 
and  this  Is  apparent  from  Formula  (4.38).  Consequently,  the  weight 
criterion,  depending  on  the  physical  properties  of  the  material  and 
proportional  to  the  weight  of  the  required  Insulation,  will  be 


/t- 

Table  4.8  gives  the  value  of  this  criterion  for  certain  heat- 
lnsulatlon  materials. 

Table  4.8  gives  three  groups  of  material:  ceramic  Insulation, 
laminated  plastics,  and  porous  Insulation  (including  loose  fibrous 
Insulation).  Within  each  group,  the  heat-lnsulatlon  materials  have 
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weight  criteria  which  are  comparatively  close.  However,  the  weight 
criteria  of  materials  from  various  groups  differ  greatly. 

Porous  materials  possess  the  lowest  weight  criteria,  but  because 
of  the  low  mechanical  properties,  they  can  be  used  only  for  internal 
insulation. 

It  should  be  noted  that  the  thermal  conductivity  of  porous  mater¬ 
ials  Increases  substantially  as  the  temperature  increases;  consequently 
the  insulation  should  be  calculated  at  an  average  temperature.  Table 
4.8  gives  the  thermal  conductivity  of  porous  materials  at  200°  C. 

TABLE  4.8. 

Ideal-Weight  Criterion  for  Certain  Heat-insulation  Materials 


1 

Marepiia^ 

X 

KKaA 

T 

KlfM^ 

c 

KfCaA  ' 

npc^eab- 
iiaa  TCM- 
neparypa 

npHMe* 

HeHiia 

2 

Jlurepa- 

rypiibrfl 

HC'^'lHHX 

M'CeK'Zpad 

Ki*zpad 

OKHCb  aJIlOMHHIfJI 
^  OKHCb  mipKOHMH 

3200 

0,30 

2.5 

2000 

1511 

4400 

0,18 

2.3 

2600 

151] 

TeKCTOJIHT 

1350 

0.35 

0,52 

250 

(251 

^  CtCK40TCKCT0JIHT 

0,7.10-4 

1600 

0,24 

0,69 

300 

(6| 

Q  Ac6ecT08biA  xap- 

^  TOfl 

0,44.10-4 

1000 

0,20 

0,47 

450 

[25] 

p  <9>ToponjiacT 

0,59.10-4 

2200 

1 

0,25 

0,72 

400 

(32| 

Ac(5acr  pacnyuieH- 

HUil 

0,22- 10-4 

100 

0.20 

1 

0.105 

1 

!  600 

(251 

11  MiiHcpaabHaa  aara 

0,18. 10-4 

150 

0,22 

0.110 

600 

(251 

n  ^  CTexaoBoaoKito 
(Maru) 

0,26.10-4 

120 

0,20 

0,125 

450 

(251 

^3  neHOcrexao 

0,2810-4 

200 

0,20 

0,160 

500 

(251 

I)  Material;  2)  temperature  limit  of  application.  In  °C;  3)  literature 
source;  4)  aluminum  oxide;  5)  zirconium  oxide;  6)  textolite;  7)  glass- 
textolite;  8)  asbestos  cardboard;  9)  fluorine  sheet;  10)  loose  asbestos; 

II)  mineral  cotton;  12)  fiberglas  (mat  finish);  13)  foam  glass. 

Ceramic  insulation  possesses  a  weight  criterion  several  times 
greater  than  the  criterion  for  laminated  plastics,  but  it  has  a  higher 
breakdown  temperature  (melting).  It  is  true  that  the  temperature  ad- 
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vantages  of  ceramic  insulation  cannot  always  be  used,  since  the  cer¬ 
amic  insulation  may  crack  during  intensive  heating  (thermal  shock). 

Like  outer  heat  insulation,  laminated  plastics  possess  good 
weight  qualities  and  have  a  comparatively  low  breakdown  temperature. 

It  should  be  noted  that  in  Table  4.8,  the  practical  temperature  limits 
of  laminated  plastics  correspond  to  the  long-term  effect  (many  hours) 
of  high  temperature.  In  the  case  of  a  short-term  effect  (seconds  and 
single  minutes),  the  permissible  operating  temperatures  for  laminated 
plastics  may  be  substantially  Increased. 

In  the  case  of  a  single  short-term  effect,  the  laminated  plastics 
may  even  be  used  at  air  temperatures  of  several  thousand  degrees;  it 
is  true,  of  course,  that  in  this  case  the  upper  layer  will  be  subject 
to  erosion  and  sublimation,  so  that  with  time  the  thickness  of  the 
plate  will  decrease.  In  addition,  it  should  be  borne  in  mind  that  after 
this  removal  of  metal  the  remaining  thickness  of  the  plate  will  have 
a  carbonized  upper  layer,  while  the  layer  beneath  will  be  broken  down 
(cracks,  folllatlon). 

For  example,  in  tests  of  glass -textollte  in  a  plasmatron  stream 
forSied  at  normal  atmospheric  pressure,  the  loss  rate  depended  on  the 
temperature  of  the  stream  and  ranged  from  0.033  to  0.8  mm/sec.;  in  this 
case  the  temperature  range  was  from  2000  to  13000®  K  [45].  For  the 
temperature  range  2000®  K  <  T  <  6000®  K,  the  loss  rate  was  almost 
linearly  dependent  on  the  temperature  and  consequently  the  thickness 
of  the  loss  can  be  represented  in  the  form 

8„==(1.175-10-<r-0,205)t  MM.  (^.39) 

The  thickness  of  the  carbonized  layer  and  damage  to  the  glass- 
textollte  structure  depend  on  the  loss  rate  (see  [45])  and  for  3300®  K 
<  T  <  6000®  K  may  be  expressed  by  the  formula 
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MH. 


(4.40) 


X 


For  T  <  3300°  K,  the  value  of  begins  to  drop. 


Manu¬ 

script 

[List  of  Transliterated  Symbols' 

Page 

No. 

86 

"^CT  ^st 

=  T  =  T 

stena  wall 

86 

Si  =  ^1  = 

^luchlstaya  ~  "^radiant 

87 

'^0(5  "  '^ob 

=  T  =  T 

obshivka  plating 

87 

Siar  *^nag  ^nagrevaniye  ~  "^heating 

92 

BH  =  vn  = 

vnutrennyy  =  inner 

100 

pasH  =  ravn  =  ravnovesiye  =  equilibrium 

100 

e$  =  ef  = 

effecktlvnyy  =  effective 

107 

M3  =  Iz  = 

izluchayemyy  =  radiated 

107 

ox  =  okh  = 

=  okhlazhdenlye  =  cooling 

113 

cp  =  sr  = 

srednyy  =  average 

118 

yr  =  ug  = 

ugar  =  loss 
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